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ABSTRACT 



Nucleate pool boiling heat transfer is an integral part 
of any vapor-compression refrigeration cycle. With a view to 
improving overall cycle efficiency, the heat transfer 
performance in the evaporator can be improved by using 
enhanced boiling surfaces. This thesis looks at the pool 
boiling characteristics of R-114 (presently used in large 
shipboard AC systems) from 10 enhanced single copper tubes and 
compares performance with a smooth copper tube. Since small 
amounts of oil escape into the refrigerant as it passes 
through the compressor of a refrigeration system, tests have 
also been conducted with up to 10% (by weight) of a miscible 
oil to see what effect this may have on overall evaporator 
performance . 
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I . INTRODUCTION 



A. BACKGROUND 

The United States Navy currently uses R-114 in 
approximately 950 centrifugal chilled-water air-conditioning 
plants. Ongoing efforts to reduce the cost and increase the 
efficiency of these plants have been directed towards the use 
of enhanced heat-transfer tube surfaces and alternative 
refrigerants. The choice of alternative refrigerant must be 
"ozone friendly." A possible "drop in" replacement for R-114 
has been identified: R-124. Studies of the boiling performance 
of R-124 would be complemented well by a baseline of data on 
R-114 . 

B. ADVANTAGES OF USING R-114 

The use of R-114 in Naval applications has been driven by 
concerns for both cost and safety. R-114 operates under 
moderate pressures and has one of the lowest toxicity ratings 
among all refrigerants. System components are significantly 
smaller and more lightweight which is desirable for shipboard 
use and maintenance. While R-ll is widely used in industry, 
its operating requirements make it unsuitable onboard ship. R- 
11 operates under vacuum and has a tendency to absorb ambient 
moisture. This acidic combination of R-ll and moisture attacks 
tubes and causes corrosion. R-114 is considerably stable to 
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heat and may be used with all standard ferrous and non-ferrous 
metals except zinc, magnesium and aluminum alloys containing 
appreciable amounts of zinc and magnesium. The thermophysical 
properties of R-114 are provided in Appendix A. 

C. THESIS OBJECTIVES 

Safety requirements for future studies of R-124 
necessitated the disassembly and relocation of the single-tube 
apparatus described in this thesis. Past studies using the 
apparatus have varied considerably in operating procedure, 
thus it became desirable to establish a consistent operating 
procedure in the acquisition of all data on R-114. With the 
foregoing discussion in mind, the objectives of this thesis 
were to: 

1. Relocate and calibrate the single tube apparatus. 

2. Operate the apparatus to produce repeatable data and 
confirm earlier results. 

3. Obtain boiling data for R-114 and R-114/oil mixtures on 
ten tube surfaces to be used as a database for future 
studies on enhanced surfaces and alternative "ozone 
safe" refrigerants. 
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II. REVIEW OF REFRIGERANT AND REFRIGERANT/OIL POOL 

BOILING BEHAVIOR 

A. NUCLEATE POOL BOILING PERFORMANCE OF PURE R-114 

The heat-transfer performance of pure R-114 on a smooth 
tube in the natural convection region of its boiling curve can 
be reasonably predicted using well known correlations, 
Churchill and Chu [Ref. 1] and more recently Churchill and 
Usagi [Ref. 2] . A lack of complete understanding of the 
nucleate boiling mechanism makes it more difficult the to 
predict the performance in the nucleate boiling region. Well 
known correlations, such as that by Rohsenow [Ref. 1], and 
more recently, Stephan-Adelsalam [Ref. 3] are approximate and 
should be compared with experimental data whenever possible. 

Many pool-boiling experiments have been carried out on 
finned tubes (typically 19 to 40 fpi) in pure refrigerants. 
For various refrigerants, tube material, fin density and 
system pressure, typical enhancements in the heat transfer 
coefficient (over smooth tubes values) of up to 3 have been 
reported, the refrigerants with better wetting characteristics 
giving the greater enhancement. Interestingly enough, the 
boiling performance of finned tubes becomes very similar to 
that of smooth tubes, if the heat flux is recalculated on the 
basis of total outside wetted surface area (rather than root 
diameter, which is commonly used) . 
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Re-entrant cavity surfaces differ from finned surfaces in 
that they provide a higher density of stable nucleation sites 
with minimum wall superheat (Yilmaz and Westwater [Ref. 4] 
Carnavos [Ref. 5] , Marto and Lepere [Ref. 6] and Wanniarachchi 
et al. [Refs. 7, 8]). Significant improvements in heat 
transfer over both finned and smooth tubes have been shown for 
re-entrant cavity surfaces. 

Boiling enhancements were as large as 14.6 and 6.4 for 
High Flux and Turbo-B tubes, respectively, compared to smooth 
tube data. 

Webb [Ref. 9] stated that the key to the high performance 
of the re-entrant surfaces can be attributed to three factors: 
(1) re-entrant cavity within a critical size range, (2) 
interconnected cavities, and (3) nucleation sites of a re- 
entrant shape. If the cavities are interconnected, adjacent 
cavities can activate each other. Re-entrant cavities provide 
a stable vapor trap, which can remain active at low values of 
superheat . 

During nucleate boiling of R-114 at one atmosphere, 
Wanniarachchi et al. [Refs. 7 and 8] reported data for four 
enhanced tubes: a porous coated High Flux, Thermoexcel-E, 
Thermoexcel-HE and a 26 fpi GEWA-T. Enhancements of 9.1, 8.2, 
6,8 and 4.4 respectively were obtained at a heat flux of 30 
kW/m 2 . The enhancement provided by the High Flux tube reduced 
to 6 at higher fluxes of about 100 kW/m 2 , but it still 
outperformed the other three tubes. At heat fluxes below 18 
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kW/m 2 , both the Thermoexcel tubes performed the best. However, 
they noted that the uncertainty in the data increased as both 
superheat and heat flux decreased. 

Whereas the majority of pool-boiling experiments use 
electrically heat tubes, McManus et al. [Ref. 10] used hot 
water to heat both a High Flux and Turbo-B tube in a pool of 
R-114. The range of heat flux was very limited. However, at a 
heat flux of 40 Kw/m 2 , enhancements of 14.6 and 6.4 
respectively, were obtained over smooth tube values. 

B. NUCLEATE POOL BOILING PERFORMANCE OF 

REFRIGERANT-OIL MIXTURES 

Most refrigeration systems use hermetically sealed 
compressors which allow small quantities of oil to escape into 
the working fluid. This oil tends to collect in the flooded 
evaporator, giving up to 10% oil in some cases. The nucleate 
pool boiling mechanism is further complicated when a less 
volatile component, such as miscible refrigerant oil, is mixed 
with the refrigerant. Design and modelling rely heavily on 
accurate experimental data covering a wide range of operating 
conditions . 

Compared with pure refrigerants, there is little data for 
boiling of refrigerant/oil mixtures from finned surfaces. 
Chaddock [Ref. 11] gives a thorough summary of such mixtures 
through 1975. Sauer et al. [Ref. 12] report data from 19 fpi 
copper tube using R-ll with two different oils. For oil 
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concentrations less than 3%, their results indicated no 
difference in boiling heat transfer performance from that of 
the pure refrigerant. For oil concentrations greater than 5%, 
significant degradation in performance was shown. They 
concluded that finned tubes performed best at low wall 
superheats and that this performance was not impaired by up to 
3% oil addition. 

Murphy [Ref. 13] carried out experiments on a 26 and 40 
fpi tube using R-114/oil mixtures. Similar to Sauer [Ref. 11], 
no difference in performance was observed for oil 
concentrations up to 3%. Indeed in some cases, at a higher 
heat flux of about 100 kW/m 2 , a 3% oil concentration showed a 
small increase in performance over the pure refrigerant case. 
The mechanism for this improvement is not well understood, but 
has been found by a number of investigators for both flow and 
pool boiling. At a 10% oil concentration, both tubes exhibited 
a 30% degradation in performance versus the pure refrigerant 
case, which agreed with Sauer's findings. 

Wanniarachi et al. [Refs. 7, 8] reported the performance 
of High Flux, Thermoexcel-E, Thermoexcel-HE and GEWA-T tubes 
in R-114/oil mixtures at one atmosphere. Oil concentrations up 
to 10% (by weight) were studied. A marked increase in the wall 
superheat values required for the incipience of nucleate 
boiling was observed in the presence of oil. Enhancements 
mentioned previously for pure R-114 were reduced at a 3% oil 
concentrat ion, to 7.5, 5.5, 5.5 and 3.8 respectively. With 10% 
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oil concentrations, further reductions to 5.5, 3.8, 3,8 and 
2.9 were observed. At the highest heat fluxes (> 45 kW/m ? ) and 
for oil concentrations greater than 6%, the performance of the 
High Flux tube degraded at a rate that left the GEWA-T tube as 
the best performer. 
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III. DESCRIPTION OF EXPERIMENTAL APPARATUS 



A. OVERVIEW OF THE SYSTEM 

Wanniarachchi et al . [Ref. 7], Murphy [Ref. 13] and 
Karasabun [Ref. 14] provide complete details of the original 
apparatus used in previous studies. Disassembly of the 
original apparatus to facilitate relocation and subsequent 
modifications to meet present site requirements justify a full 
description of the apparatus. The original R-12 cooling 
subsystem was supplanted by a more efficient R-502 system. The 
auxiliary heaters and plastic shield were removed, however 
power extensions were kept available for their use in future 
studies. The Variac power supply to the tube heater was 
calibrated and emplaced into a panel enclosure. See Appendix 
D, Apparatus Calibration for more details. A more modern 
Hewlet-Packard Computer/Data Acquisition system was installed 
and the data reduction program DRP8 was updated to account for 
variac calibration and program language changes. 

An overall schematic representation of the experimental 
apparatus is shown in Figure 3.1, and a photograph is shown in 
Figure 3.2. The apparatus essentially consists of 8 
(components: 1) a Pyrex-glass tee for the pool boiling of the 
R-114 liquid, "the evaporator"/ (2) a Pyrex-glass tee for the 
condensing of the R-114 vapor, "the condenser"; (3) a R-114 
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liquid reservoir; (4) a refrigerant /oil subsystem composed of 
an oil reservoir and a graduated oil cylinder; (5) a cooling 
subsystem composed of a 1/2 ton R-502 refrigeration plant, 2 
positive displacement pumps and a 30 gallon water/ethylene 
glycol sump; (6) a vacuum pump; (7) a data acquisition and 
instrumentation system; (8) an aluminum/Plexi-glass framework 
within which components (1) through (4) were housed. 

The apparatus was operated as follows: liquid R-114 was 
boiled in the evaporator (1), passed as vapor through the 
aluminum L-shaped connection and condensed in the condenser 
(2) . R-114 condensate was returned by gravity via the 
distribution tube to the condenser. The R-502 refrigeration 
plant maintained the water/ethylene glycol mixture in the sump 
between -12 and -18 °C. Cooling of the sump was accomplished 
by a countercurrent flow heat exchanger through which an 8-gpm 
turbine-type pump recirculated the water/ethylene glycol 
mixture. An identical positive displacement pump circulated 
the cooled mixture through the copper condensing coils in the 
condenser via the control valve VC. 

Prior to dismantling the evaporator in order to change out 
boiling tubes, the liquid R-114 was boiled out of the 
evaporator, condensed on the coils of the condenser then 
recovered to the R-114 liquid reservoir rather than returned 
to the evaporator. Upon reassembly, the evaporator was 
recharged with liquid R-114 by gravity from the reservoir. The 
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level in the reservoir could be topped off via the condenser 
from an external supply cylinder through valve V10. 

Exact quantities of oil could be added to the liquid R-114 
in the lower glass tee from the graduated cylinder (5) by 
gravity. The graduated cylinder was filled with oil by gravity 
from the oil reservoir through valves V14 and V2 . Refilling 
the oil reservoir was done through valve V3 located at the top 
of the framework. All oil left over from the recovery of the 
liquid R-114 was drained out of the lower glass tee and 
disposed of. 

The boiling tube was mounted horizontally in the 
evaporator and held in place by Teflon bushings with O-ring 
inserts. The whole system was connected to the vacuum pump (6) 
through valve V8 in order to remove noncondensible gases. A 
relief valve mounted on the aluminum L connection was set to 
a gage pressure limit of 20 psi, which is 50% less than the 
manufacturer-recommended working pressure limit of the Pyrex 
glass tees. All connections were assembled with Swageiok 
fittings sealed with Teflon ferrules to ensure leak tightness. 
Copper tubing of 3/8 inch diameter was used for all piping, 
except in the oil section, where 1/4 inch copper tubing was 
used . 
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B. BOILING TEST SECTION 



1 . Evaporator 

The evaporator is a T-shaped container made of Corning 
Pyrex glass. Figure 3.3 shows a schematic of the evaporator 
with dimensions. Pyrex glass has several advantages: it is 
corrosion-resistent , transparent, has a smooth interior 
surface (this minimized nucleate boiling at the inner surface 
of the container) and is stronger with temperature and 
pressure variations compared to ordinary glass. The 
manufacturer recommended operating pressure limit of the Pyrex 
glass tee is 30 psi gage. Each end of the evaporator was 
fitted with a cast-iron flange and a gasket. A detailed sketch 
of the cast-iron flange is shown in Figure 3.3 also. To reduce 
weight upon the glass tees, two aluminum flanges, 210 mm in 
diameter and 12.7 mm in thickness, were bolted to the cast- 
iron flanges. All fittings were connected through the aluminum 
flanges to the Pyrex-glass evaporator. 

2. Test Tube 

A schematic drawing of the test tube is shown in 
Figure 3.4. The boiling tube was held in place by two Teflon 
bushings, which were attached to the aluminum flanges at both 
ends of the T-shaped evaporator. Four studs were used in order 
to attach the Teflon bushing on the aluminum flange. The 
Teflon bushing and both ends of the test tube were sealed by 
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means of two viton O-rings that were placed between the 
aluminum flange and the Teflon plug at each end. 

3 . Boiling Tubes 

A smooth, hard-copper tube, 15.9 mm (5/8 in.) in outer 
diameter, 12.7 mm (1/2 in) in inside diameter and 431.3 mm (17 
in.) in length, was used to provide comparison data for the 
nucleate pool boiling heat-transfer performance of the 
enhanced tubes tested. All of the tubes were heated and 
instrumented identically. Dimensions of each tube are provided 
in Table V of Appendix E, Sample Calculations. Wanniarachchi 
et al. [Refs. 8, 9] and Murphy [Ref. 14] provide complete 

details of the tube surfaces. 





The ten tubes 


tested 


were : 


1 . 


Smooth 








2. 


GEWA-K 26 


fpi 


(Figure 


3.5) 


3. 


GEWA-K 40 


fpi 


(Figure 


3.5) 


4. 


GEWA-T 19 


fpi 


(Figure 


3.5) 


5. 


GEWA-T 26 


fpi 


(Figure 


3.5) 


6. 


GEWA-YX 26 


fpi 


(Figure 


3.6) 


7. 


High Flux 


(Figure 3.7) 




8. 


Thermoexcel-HE 


(Figure 


3.8) 


9. 


Thermoexcel-E 


(Figure 


3.9) 


10. 


Turbo-B (Figure 3.9) 
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The heater was a 1000-Watt, 240-Volt stainless-steel 
cartridge, 6.35 mm (1/4 in.) in outer diameter and 203.2 mm (8 
in) in length. It was inserted into a copper sleeve, which was 
1/4 inch in inside diameter, 1/2 inch in outer diameter and 8 
inches in length. In order to provide a uniform heat flux, the 
cartridge heater and the copper sleeve were soldered together. 
To minimize the thermal contact resistence, the copper sleeve 
was then tinned (the cartridge heater was used to maintain the 
solder in molten state) and inserted as a unit into the middle 
portion of the test tube. The active boiling length of the 
test tube was therefore eight inches in the middle portion of 
the boiling tube. In order to compute the actual average heat 
flux in the heated portion of the tube, a suitable correction 
was applied for the heat lost by natural convection out both 
ends. See Appendix E, Sample Calculation for details. 

To measure the wall surface temperature of the boiling 
tube, 8 thermocouples (Type-T Teflon coated copper-constantan; 
0.01 inch in diameter) were inserted into eight grooves 
machined on the outside of the copper sleeve. As shown in 
Figure 3.10, these thermocouples were located at different 
axial and circumferential locations. The exact locations of 
the thermocouples and dimensions of the thermocouple grooves 
are given in Figure 3.10. All the thermocouple grooves were 
axially machined from the location of the thermocouple hot 
junctions to the nearest end of the copper sleeve. 
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C . CONDENSER SECTION 



The condenser was also a T-shaped container made of 
Corning Pyrex glass. It was identical to the evaporator. The 
position of the condenser can be seen in Figure 3.1. R-114 was 
condensed on a helical copper coil, which was inserted in the 
Pyrex-glass condenser. Copper tubing (3/8 inch outer diameter) 
was fabricated into a three inch diameter coil. The active 
condensation length was estimated to be 15 ft. 

The top portion of the condenser was connected to a 
portable, mechanical vacuum pump to remove noncondensible 
gases from the apparatus. The bottom of the condenser was 
connected to the evaporator via valve V5 in order to return 
the condensed R-114 liquid to either the evaporator or the 
reservoir. The coolant, i.e., water/ethylene glycol mixture, 
entered from the top portion of the condenser. The condenser 
was placed vertically and connected to the vapor outlet of the 
evaporator using a fabricated, L-shaped aluminum tube, two 
inches in diameter. A Bourdon gage with a range up to a gage 
pressure of 150 psi and a relief valve which was set to 20 psi 
gage were mounted on the L-tube. 

D. OIL ADDING SECTION 

To study the boiling performance of R-114/oil mixtures, a 
cylindrical aluminum oil reservoir, 6 inches in diameter and 
6 inches in height, and a graduated glass oil cylinder were 
installed above the evaporator as shown in Figure 3.1. The 
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graduated cylinder was 355 mm in length and had a diameter of 
25.4 mm. A resolution of 0.5 ml was achieved with the 
graduated cylinder. The oil reservoir was connected to the oil 
cylinder through valves V14 and V2 . Refilling the oil 
reservoir was accomplished through valve V3 located on top of 
the reservoir. The addition of oil into the evaporator was 
achieved through VI by gravity. The oil cylinder was refilled 
by gravity from the oil reservoir located above. 

E . COOLING SECTION 

1 . Water-Ethylene Glycol Mixture Tank 

In order to store the water-ethylene glycol mixture, 
a special tank was manufactured. The total volume of the tank 
was 9243 in 3 (18.84 in x 18.84 in x 26.04 in) and it was made 
of 0.5 inch thick Plexiglas sheet. All sides of the tank were 
glued together with methylene chloride solution. The joints 
were held together with small screws for extra strength. The 
low thermal conductivity of Plexiglas was especially suited to 
minimize heat transfer (from room to water-ethylene glycol 
mixture) through the tank walls. The tank was placed on top of 
square wood platform on the floor and all sides were insulated 
with 7/8 in thick rubber sheets. The cooling mixture contained 
13 gallons of ethylene-glycol and 25 gallons of distilled 
water. The freezing point of this mixture was about -25 deg C. 
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2 . R-502 Refrigeration Plant 

A 1/2 ton R-502 refrigeration plant was installed to 
cool the water/ethylene glycol mixture. Figure 3.12 shows a 
schematic of the R-502 refrigeration plant. It consists of a 
compact-type, air-cooled condenser, a compressor, a receiver, 
a filter-drier unit, a pressure regulator, a temperature 
control switch and a thermostatic expansion valve. The 
evaporator of the R-502 refrigeration plant was constructed 
using a 3.8 inch copper tube which was run through a 
countercurrent heat exchanger opposite the water-ethylene 
glycol mixture. The temperature of the water /ethylene glycol 
mixture was controlled by both a thermostatic expansion valve 
and a temperature control switch. The R-502 refrigeration 
plant was adjusted to keep the temperature of the cooling 
liquid at about -17 deg C. 

3 . Pump and Control Valve 

An 8 gpm, 115 VAS Burks turbine-type, positive- 
displacement pump was installed on the floor and the 1 inch 
diameter suction side of the pump was directly connected to 
the water/ethylene glycol tank. Cooling liquid was pumped from 
the tank to the condenser through the control valve VC. Also, 
a by-pass valve V9 was placed before the control valve VC on 
the discharge line. The use of the by-pass line served two 
important purposes: (1) it avoided overloading the pump in the 
event valve VC is completely closed, and (2) it provided 
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proper mixing for the "warm" stream returning from the 
condenser. The by-pass valve required adjustment only at the 
highest heat fluxes to 'satisfy' proper cooling of the 
condenser . 

F. R-114 RESERVOIR 

An aluminum cylindrical reservoir, 9 inches in diameter 
and 10 inches in height, was placed vertically between the 
evaporator and condenser in order to store R-114 as a liquid. 
The liquid level of the R-114 could be observed by means of a 
sight glass attached on the reservoir. The R-114 reservoir was 
connected to the vapor line through valve V7 and to the liquid 
line through valve V6. See Figure 3.1 for the arrangement of 
the reservoir. 

6 . CHAMBER 

An aluminum frame (42.13 in x 20.07 in x 24.02 in) was 
constructed to locate all the parts of the apparatus, except 
the cooling section. All four vertical sides of the frame were 
covered with 1/2 inch thick Plexiglas sheets and both left and 
right sides were provided with hinges to enable easy access to 
the components of the apparatus. Aluminum and plywood plates 
were used to cover the bottom and top sides of the frame, 
respectively. The valve bodies of VI through V8 were placed 
inside of the front Plexiglass box. The whole frame was placed 
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above the water-ethylene glycol tank with aluminum supports so 
that the system was very compact. 

One of the main advantages of this chamber was that the 
temperature surrounding the evaporator could be reduced 
relative to the ambient temperature. Also, in case of 
emergency, the thick Plexiglas chamber provides a safety 

barrier to personnel and data acquisition equipment. 

H. INSTRUMENTATION 

1 . Power Measurement 

A 240 volt AC source was used as the power supply, and 
it was adjusted by a variac in the range of 0-220 volts and 
0-5 amperes according to the desired heat flux at the surface 
of the boiling tube. Power input to the boiling tube was 

measured with an AC current inductive sensor (output in volts) 
and a voltage sensor directly coupled to the variac output. 
The voltage sensor was run through a AC-DC true R.M.S. 
converter which output a proportional signal in volts. Figure 
3.12 shows a schematic representation of the power-measurement 
setup. Both the AC current sensor and the AC-DC true R.M.S 
converter were connected to the data acquisition/control unit. 

2 . Temperature Measurement 

Various temperatures were monitored throughout the 

system to include: 

1. Boiling tube wall (8 thermocouples at various 

longitudinal and circumferential positions) 
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2. Liquid temperature (one thermocouple) 

3. Vapor temperature (two thermocouples), and 

4. Water/ethylene glycol mixture temperature (one 

thermocouple) 

The locations of the wall thermocouples in the sleeve of the 
boiling tube are shown in Figure 3.5. The liquid and vapor 
thermocouples were inserted into the three specially 
manufactured thermocouple housings. Figure 3.14 shows a 
schematic of these thermocouple housings. The stainless-steel 
portion minimizes (owing to low thermal conductivity) errors 
resulting from the axial conduction of heat from the 
surroundings. The copper tip of the housing helps to minimize 
the temperature drop from the area being measured to the 
thermocouple location (owing to the high thermal conductivity 
of copper) . 

All the temperature measurements were accomplished by 
30 gage copper-constantan thermocouples. Each thermocouple 
measurement was read directly by a Hewlett-Packard 3497A data 
acquisition system, which was controlled by a Hewlett-Packard 
9826 computer. Each thermocouple was scanned twenty times over 
5 seconds and the readings averaged to obtain a more accurate 
measurement . 
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Figure 3.1 Schematic of Single Tube Apparatus 



20 




Figure 3.2 Photograph of Experimental Apparatus 
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a) Corning Pyrex Glass Evaporator (D x Dj. = 402x51 mm, 



L = 178 mm, L^= 127 mm) 
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b) Cast Iron Flange and Gasket (d^ 
L = 14mm, A = 21°) 



190mm, d 0 = 210mm, 



Figure 3.3 Schematic of the Pyrex Glass Evaporator 
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Figure 3.4 Schematic of the Boiling Test Tube 





(a) Standard fin (GEWA-K) and GEWA-T fin 




(b) GEWA-T fin 



Figure 3.5 



Schematic and Photograph of Fins 
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(a) GEWA-YX tube 




(b) GEWA-YX fin 

Figure 3.6 Photographs of GEWA-YX Tube and Fin 
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VAPOR BUBBLE TRAPPED 
IN ACTIVE SITE 




Figure 3.7 Schematic and Photograph of High Flux Surface 
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Figure 



3.8 Schematic and Photograph of 
Thermoexcel-HE Surface 
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(a) Thermoexcel-E 




(b) Turbo-B 



Figure 3.9 Photographs of Re-Entrant Cavity Surfaces 
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SECTION A-A (DIMENSIONS IN MM ) 
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Figure 3 11 Schematic of R-502 Ref r i ger at i on System 
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POWER SUPPLY 93DD 



THERMOCOUPLE 
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Figure 3.13 Sketch of a Thermocouple Well 



IV. DATA ACQUISITION AND REDUCTION 



A. DATA ACQUISITION AND STORAGE 

A Hewlett-Packard 3852A automatic data acquisition unit 
was used to read temperatures from the thermocouples and to 
read current and voltage values of the cartridge heater from 
the AC current sensor and the AC-DC true RMS converter. A 
Hewlett-Packard 9000 Series computer was used to control the 
measurements done by the 3852A, analyze the collected data and 
store it. 

The iterative data collection/reduction program DRP8 was 
loaded and run. Information was entered through the keyboard 
to prompt the data acquisition unit when to take data. Channel 
assignments are listed in Table I. The raw data was processed 
and transferred to a file on the computer's hard drive. A 
printout of the data reduction was also provided. 

B . DATA REDUCTION 

Following data acquisition for each data point, results 
were computed according to the stepwise procedure outlined in 
the next section, and then printed out using a Hewlett-Packard 
Inkjet printer. Graphs were plotted using commercial software. 
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C. STEPWISE DATA-COLLECTION AND SOLUTION PROCEDURE 



1. Input name of user-specified file to be stored on 
computer hard drive. 

2. Select tube type (all dimensions of the boiling test 
tube are included) . 

3. Input desired heat flux (W/m 2 ) and saturation 

temperature (°C) of the boiling liquid. 

4. Set desired heat flux by adjusting VARIAC rheostat. 

5. Set saturation temperature by adjusting flow of coolant 
through condenser coils with control valve VC. 

6. Once saturation temperature is achieved, wait for 
steady-state conditions (at least 5 minutes) prior to 
taking data readings. 

7. Prompt data acquisition unit to scan all channels 

listed in Table I (thermocouples, AC current sensor and 
AC-DC true RMS voltage converter) . All channel readings 
are made in volts and stored in user specified file. 

8 . Convert these voltage readings to corresponding units 
(temperature in °C, current in amperes, voltage in 
volts) . Karasabun [Ref. 15] describes in detail 
conversion equation for temperature. See Appendix 13- 
Data Acquisition Apparatus Calibration for details of 
current and voltage conversion. 

9. Compute the heat-transfer rate from the cartridge 

heater. See Appendix D-Sample Calculations. 

10. Compute the average wall temperature of the boiling 
tube and calculate the wall superheat of the current 
data set. See Appendix D-Sample Calculations. 

11. Compute the physical properties of R-114 using given 
correlations at film temperature. See Appendix D-Sample 
Calculations . 

12. Compute the natural-convection heat-transfer 

coefficient of R-114 from the smooth non-boiling ends 
of the test tube. See Appendix D-Sample Calculations. 

13. Compute heat losses from the non-boiling ends of the 
boiling tube. See Appendix D-Sample Calculations. 
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14. Calculate the heat flux from tube to the refrigerant. 
See Appendix D-Sample Calculations. 

15. Calculate the heat transfer coefficient of the R-114 
from the boiling tube. See Appendix D-Sample 
Calculations . 

16. Store the heat flux versus wall superheat values for 
each data set in user-specified file. 

17. Use commercial software. Harvard Graphics, to plot 
data . 



TABLE I. IIP 3497A CHANNEL ASSIGNMENTS 



Channel 


Assignment 


500-507 


Boiling tube wan temperature 


508-509 


Boiling liquid temperature 


510 


Vapor temperature 


511 


Sump temperature 


512 


RMS voltage Tube heaters 


513 


Tube heater current sensor 


514 


Auxiliary heater current 
sensor 
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V. EXPERIMENTAL PROCEDURE 



A . PREPARATION 

1 . Vacuum test of the Apparatus 

Upon changing the tube and reassembling the apparatus, 
the system was evacuated down to approximately 25 in. Hg. by 
a portable mechanical vacuum pump through valves VI 1 and V8 . 
The apparatus was allowed to stand for 30 minutes under vacuum 
to observe any possible leaks. During this period, 
thermocouples were connected to the data acquisition system, 
the electrical integrity of the tube heater was checked and 
the tube heater connected to the variac panel. If a leak was 
detected (i.e. an observed rise in pressure), a pressure test 
was conducted to isolate the source. (See V.A.2.) 

2 . Pressure Test of the Apparatus 

In order to detect sources of leaks, the apparatus was 
pressurized with air up to 15 psi gage and a soap/water bubble 
test was carried out upon suspect joints. All detected leaks 
were successfully and systematically fixed. 

3. Charging the Apparatus with R-114 

Upon completing the pressure tests (if any) , the 
evaporator of the apparatus was filled with R-114 liquid up to 
a marked position (20 mm above the top of the boiling tube) 
using the following procedure. 
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1. Valve V7 was opened to allow the pressures within the 
reservoir and the evacuated apparatus to equalize. It 
was noted that the liquid R-114 within the reservoir 
boiled since the equilibrium pressure was less than the 
R-114 vapor pressure at room temperature. 

2. Valve V6 was opened to allow the liquid R-114 to drain 
into the evaporator by gravity. 

3. At the designated evaporator level, valves V6 and V7 
were closed to isolate the reservoir. 



4. Degassing and Data Acquisition Channel Check. 



1 . With the condenser operating, the boiling tube was run 
at a high flux of about 100 kW/m 2 for ten minutes in 
order to de-gas both the tube surface and refrigerant. 
Care was taken not to exceed a pressure of 15 psi gage 
(well below manufacturer recommendation) within the 
apparatus . 

2. Any noncondensible gases present collected in the top 
of the condenser. These were then removed using the 
vacuum pump through valves, V8 and VI 1. 

3. The program, SETUP8, was run to check output on all 
thermocouple channels as well as power output. If an 
erroneous reading was detected, then appropriate action 
was taken. In the case of a faulty boiling tube wall 
thermocouple, the faulty thermocouple was tagged for 
future exclusion in data acquisition procedure. (Each 
thermocouple had X thermocouples embedded in its tube 
wall and the loss of one made no difference to the 
average wall temperature) . 

4. The apparatus was shut down and the boiling tube 
allowed to soak in the refrigerant overnight. This 
allowed for the surface to become wetted. 



B. NORMAL OPERATION 

The following procedure was followed to obtain the heat- 
transfer coefficient of pure R-114 and R-114/oil mixtures from 
the boiling tube being tested: 
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1. The R-502 refrigeration unit was operated 1-2 hours in 
advance in order to reduce the temperature of the 
water-ethylene glycol sump to a minimum of -8 °C. 

2. The data acquisition/control unit, computer and variac 
panel were switched on. 

3. The computer program Setup8 was loaded and run: 

a. All the data acquisition channels were rechecked. 

b. A power output of approximately 5W was input to 
the boiling tube. This has been measured to 
produce a heat flux of 500 to 600 W/m 2 in all tubes 
tested which facilitated reaching the first data 
point on the increasing flux run. 

c. The temperature of the refrigerant was slowly 

reduced down to 2.2 °C by circulating a small 
amount of coolant through the condenser, regulated 
by control valve, VC. This reduction was monitored 
to coincide with the reduction in sump 

temperature . 

4. The iterative computer program DRP8 was loaded and 
run . 

5. Two data runs were made: an increasing heat flux run 
and a decreasing heat flux run wherein the heat flux 
was incrementally increased or decreased. Ten different 
heat fluxes (0.6, 1, 2, 3.5, 6, 10, 20, 35, 60, and 100 
kw/m 2 ) were selected for each increasing or decreasing 
heat flux data run. 

6. The desired heat flux and saturation temperature were 
input to the program as reference points for each data 
point . 

7. The variac was adjusted to set the actual heat flux 
which was continuously compared to the desired heat 
flux by the program until they agreed within 2 percent. 

8. The control valve VC was adjusted to regulate the flow 
of cooling liquid through the condenser to maintain 
nearly constant saturation temperature at a given heat 
flux. Desired versus actual saturation temperatures 
were monitored continuously by the program until they 
agreed ± 0.1 °C. 

9. For each data point, conditions in the evaporator were 
allowed to stabilize for at least 5 minutes prior to 
measuring the raw data. The following raw data were 
measured and stored in a user specified file: boiling 
tube wall thermocouple readings, liquid bulk 
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thermocouple readings, vapor thermocouple readings, 
sump thermocouple readings, current sensor readings, 
and voltage sensor readings. 

10. Two data points were taken at a given heat flux and 
saturation temperature to comprise a data set. The 
following processed data were recorded as a printout: 
wall temperatures of the boiling section, liquid bulk 
temperatures, vapor temperature, sump temperature, wall 
superheat, heat transfer coefficient of the R-114/R- 
114-oil mixture and the actual heat flux. 

11. For each data set, the above procedure beginning with 
step 5 was repeated. 

12. Table II provides a listing of all data runs. 



TABLE II. LISTING OF DATA RONS 



Data 
Run # 


Data File 


Tube 


Heat 

Flux 


Oil 

(%) 


Purpose 


1 


DATA630I5 


HF 


I 


0 


CAL 


2 


DATA630D5 


HF 


D 


0 


CAL 


3 


DATA702I5 


HF 


I 


0 


CAL 


4 


DATA702D5 


HF 


D 


0 


CAL 


5 


DATA0705I5 


HF 


i 


0 


CAL \ DATA 


8 


DATA0705D5 


HF 


D 


0 


CAL \ DATA 


7 


DATA07 0 914 


SM 


I 


0 


DATA 


8 


DATA070 9D4 


SM 


D 


0 


DATA 


9 


DAT0710I6 


GK-4 0 


I 


0 


DATA 


10 


DAT0710D6 


GK-4 0 


D 


0 


DATA 


11 


DAT0711I7 


GK-26 


I 


0 


DATA 


12 


DAT0711D7 


GK-26 


D 


0 


DATA 


13 


DAT 071518 


GT-19 


I 


0 


DATA 


14 


DAT0715D8 


GT-19 


D 


0 


DATA 


15 


DAT0716I9 


GT-26 


I 


0 


DATA 
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Data 
Run # 


Data File 


Tube 


Heat 

Flux 


Oil 

(%) 


Purpose 


16 


DAT0716D9 


GT-26 


D 


8 


DATA 


17 


DAT0717I11 


T-HE 


i 


% 


DATA\DG 


18 


DAT0717D11 


T-HE 


D 


0 


DATA 


19 


DAT0722I12 


TB 


I 


0 


DATA 


20 


DAT0722D12 


TB 


D 


0 


DATA 


2 % 


DAT0723I09 


GY-26 


I 


0 


DATA 


22 


DAT0723D09 


GY-26 


D 


0 


DATA 


23 


DAT0724 I 1 1 


T-HE 


I 


0 


DATAXDG 


24 


DAT0724D11 


T-HE 


D 


0 


DATA \ REP 


25 


DAT0803I4 


SM 


I 


0 


DATA\REP 


26 


DAT0803D4 


SM 


D 


$ 


DATA\REP 


27 


DATO 8 0414 


SM 


I 


• 


DATA 


26 


DAT0804D4 


SM 


D 


• 


DATA 


29 


DAT080514 


SM 


I 


10 


DATA 


30 


DAT0805D4 


SM 


D 


10 


DATA 


31 


DATO 808 17 


GK-2 6 


I 


0 


DATA\REP 


32 


DAT0806D7 


GK-2 6 


D 


0 


DATA\REP 


3* 


DAT0807I7 


GK-26 


I 


3 


DATA 


34 


DAT080ID7 


GK-26 


D 


• 


DATA 


35 


DAT0808I7 


GK-26 


I 


10 


DATA 


36 


DAT0808D7 


GK-26 


D 


10 


DATA 


37 


DAT0809I6 


GK-40 


I 


8 


DATA\REP 


38 


DATO 8 H ©D6 


GK-4 0 


D 


0 


DATA\REP 


36 


DAT0810I6 


GK-4 0 


I 


3 


DATA 


40 


DAT0810D6 


GK-40 


D 


3 


DATA 


41 


DAT 081216 


GK-4 0 


I 


10 


DATA 


42 


DAT0812D6 


GK-4 0 


D 


10 


DATA 


43 


DATO 8 151 8 


GT-19 


I 


0 


DATA\REP 
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Data 
Run # 


Data File 


Tube 


Heat 

Flux 


Oil 

(%) 


Purpose 


44 


DAT0815D8 


GT-19 


D 


5 


DATA\REP 


45 


DAT0816I8 


GT-29 


I 


3 


DATA 


46 


DAT081608 


GT-19 


D 


3 


DATA 


4 1 


DATO 82018 


GT-19 


I 


10 


DATA 


48 


DAT0820D8 


GT-19 


D 


10 


DATA 


49 


DAT0821I9 


GT-26 


I 


3 


DATA 


50 


DATO 8 2 ID 9 


GT-26 


D 


3 


DATA 


51 


DATO 822 I 9 


GT-26 


i 


10 


DATA 


52 


DAT0822D9 


GT-26 


D 


10 


DATA 


53 


DAT 0823110 


T-E 


I 


3 


DATA 


50 


DAT0823010 


T-E 


D 


3 


DATA 


55 


DATO 82 411 0 


T-E 


I 


10 


DATA 


56 


DAT0824D10 


T-E 


D 


10 


DATA 


57 


DAT082519 


GY-26 


I 


3 


DATA 


56 


DAT0825D9 


GY-26 


D 


3 


DATA 


56 


DATO 82 61 9 


GY-26 


i 


10 


DATA 


60 


DAT0828D5 


GY-26 


D 


10 


DATA 


61 


DATO 821 15 


HF 


i 


3 


DATA 


62 


DAT0827D5 


HF 


D 


3 


DATA 


63 


DATO 82 8 15 


HF 


i 


10 


DATA 


64 


DAT0828D5 


HF 


D 


10 


DATA 


61 


DATO 82 911 2 


TB 


I 


3 


DATA 


66 


DAT0829D12 


TB 


D 


3 


DATA 


61 


DAT 0830112 


TB 


I 


10 


DATA 


68 


DAT 0 8 3 0 D 1 2 


TB 


D 


10 


DATA 


69 


DATO 8 31 111 


T-HE 


I 


3 


DATA 


70 


DAT0831D11 


T-HE 


D 


3 


DATA 


71 


DATO 90 1 111 


T-HE 


I 


10 


DATA 
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Data 
Run # 


Data File 


Tube 


Heat 

Flux 


Oil 

(%) 


Purpose 


72 


DAT0901D11 


T-HE 


D 


10 


DATA 


73 


DAT0906I10 


T-E 


i 


0 


DATA 


74 


DAT0906D10 


T-E 


D 


0 


DATA 



SM = Smooth GK-26 = GEWA-K 26 fpi GK-40 = GEWA-K 40 fpi 

HF = High Flux GT-19 = GEWA-T 19 fpi GT-26 = GEWA-T 26 fpi 

TB = Turbo-B GY-26 = GEWA-YX 26 fpi T-E = Thermoexcel-E 
T-HE = Thermoexcel-HE 

I = Increasing Flux D = Decreasing Flux 

CAL = CALIBRATION DATA DATA = BASELINE DATA 

REP = REPEATABILITY DATA DG = DEGASSING DATA 
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VI . RESULTS AND DISCUSSION 



A. REPRODUCIBILITY 

To test the reproducibility of the experimental data, 9 
data runs were repeated on different days. All data runs were 
conducted at a saturation temperature of 2.2 °C. Figure 6.1 
shows a comparison of two pairs of runs, 13, 14, and 43, 44 
for increasing and decreasing heat flux. It can be seen that 
there is excellent agreement between both the increasing and 
decreasing heat flux runs. 

B. COMPARISON OF DATA WITH PREDICTION 

To check the general validity of the data obtained, four 
well-known correlations, two for natural convection and two 
for nucleate pool boiling, were plotted and compared with the 
smooth tube data for pure R-114. The two natural convection 
correlations used were Churchill and Chu [Ref. 1] and a more 
recent update Churchill and Usagi [Ref. 2] . The two nucleate 
pool boiling correlations used were: Rohsenow [Ref. 1] and 
Stephan-Abdelsalam [Ref. 3]. One of the first correlations 
developed, the Rohsenow correlation is based upon the physics 
of pool boiling heat transfer. It applies only for clean 
surfaces and when used to estimate heat flux, errors can 
amount to ± 100%. The Stephan-Abdelsalam correlation for 
refrigerants was developed from application of regression 
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analysis to a large amount of existing data. Application to a 
large number of pure substances of widely differing thermal 
properties is remarkably simple. Figure 6.2 shows the plotted 
data and correlations. It can be seen from the plot that there 
is excellent agreement for natural convection and not bad 
agreement in the nucleate boiling region. The slightly worse 
agreement for nucleate boiling can be attributed to the more 
complex mechanisms at work in this region. 

C. EFFECTS OF DEGASSING 

Early experiments (Runs 1-6) conducted with the High Flux 
tube pointed out the significant effect of degassing both the 
tube surface and refrigerant. It was observed that tubes that 
were not vigorously boiled after placement into the apparatus 
and not allowed to settle overnight, tended to nucleate at 
very low heat fluxes the following day when conducting an 
increasing heat flux run. Similar effects were observed if the 
apparatus was allowed to stand undisturbed for several days. 
The most extreme effects were seen for the re-entrant cavity 
tubes. Figure 6.3 shows the change in behavior (i.e. for 
gassed and degassed) for a Thermoexcel-E tube for an 
increasing heat flux run. For the degassed case, an extreme 
temperature overshoot was observed, which indicates a large 
number of nucleation sites were deactivated. Once all 
nucleation sites were activated, the performance of both cases 
was identical. It can be surmised that the entrapment of vapor 
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bubbles in surface cavities as well as the outgassing from the 
Teflon O-rings of the rig over time, enhanced the nucleation 
on the tube surfaces. Vigorous boiling removed those gases 
from both surface and liquid. 

D. BOILING PERFORMANCE OF ENHANCED TUBES IN PURE R-114 

AND R-114 /OIL MIXTURES 

The results of the 74 data runs are discussed with 
references to the figures (Figures 6.1 to 6.80) that are 
provided at the end of this chapter. Performance of each tube 
will be presented and discussed as follows: 



1. For pure R-114 

2. For R-114/oil mixtures (3% and 10% by weight) 



3. 


Comparison 

increasing 


of pure R-114 and 
heat flux runs. 


R-114/oil 


mixtures 


for 


4 . 


Comparison 

decreasing 


of pure R-114 and 
heat flux runs. 


R-114/oil 


mixtures 


for 


5. 


Comparison 

applicable 


to similar tube types (i.e. finned ■ 
for each R-114/oil mixture for 


- as 
both 



increasing and decreasing heat flux runs) 



1 . Boiling Performance of Smooth Tube 

Figures 6.4, 6.5 and 6.6 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from a smooth copper tube. 
Figures 6.7 and 6.8 show the comparative performance for the 
three ref rigerant /oil mixtures on increasing and decreasing 
heat flux respectively. It can be seen that in Figure 6.7 in 
the natural convection region there is slight degradation with 
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increases in oil. However, in the nucleate boiling region 
there is significant degradation with increasing oil 
concentration. Wanniarachi et al. [Ref. 6] reported similar 
behavior for smooth tubes. This behavior was attributed to two 
competing effects: (1) the 'choking' of the tube surface by an 
oil rich layer which develops as the more volatile refrigerant 
is boiled off; () the diffusion of the miscible oil rich layer 
back into the bulk fluid. In the natural convection region it 
can be surmised that the diffusion effect predominates while 
in the nucleate boiling region the 'choking' effect eventually 
overwhelms the diffusion. This degradation is seen again in 
Figure 6.8. It can be seen that in Figure 6.8 that as 
conditions approach natural convection along the decreasing 
heat flux curve, the three mixtures perform similarly as 
nucleation sites die out. 

2 . Boiling Performance of Finned/Modified Finned Tubes 

Figures 6.9, 6.10 and 6.11 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from a standard 26 fpi tube, 
(GEWA-K) . At a heat flux of 35 kW/m 2 , heat transfer 
enhancement compared to the smooth tube for pure R-114, 
R-114/3% oil and R-114/10% oil was 2.4, 2.9 and 3.2 
respectively. Figures 6.12 and 6.13 show the comparative 
performance of the three ref rigerant /oil mixtures on 
increasing heat flux and decreasing heat flux runs 



46 



respectively. On increasing heat. flux and decreasing heat flux 
runs, the finned tube only shows a small degradation in both 
the natural convection and nucleate boiling regions. The 
smooth tube only shows small degradation in the natural 
convection region. Also, at the higher heat fluxes, the 
R-114/3% oil mixture shows a slightly better performance than 
the pure R-114 and R-114/10% oil mixture. This phenomenon was 
observed for all finned tubes, however, the reasons for such 
behavior are not fully understood; it may be due to bubble 
'scouring' between the fins. This small increase in 
performance at low oil concentrations agrees with results 
reported by Stephan and Mitrovic [Ref. 15] . 

Figures 6.14, 6.15 and 6.16 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from a GEWA-K 40 fpi tube. 
Figures 6.17 and 6.18 show the comparative performance of the 
three ref rigerant /oil mixtures on increasing heat flux and 
decreasing heat flux respectively. Similar behavior to the 
GEWA-K 26 fpi tube was observed including the improved 
performance at higher fluxes at low oil concentrations (3%) . 

Figures 6.19, 6.21, and 6.23 show the comparative 
performance of the three R-114/oil mixtures boiling from GEWA- 
K 26 and 40 fpi tubes on increasing heat flux. For all three 
oil concentrations, the GEWA-K 40 fpi tube outperforms the 
GEWA-K 26 fpi and smooth tubes consistently in the nucleate 
boiling region. At a heat flux of 35 kW/m 2 the GEWA-K 40 fpi 
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and the GEWA-K 26 fpi had heat transfer enhancements of 4.1 
and 3.1 respectively over the smooth tube for pure R-114. 
Similar enhancements were seen at R-114/3% oil and R-114/10% 
oil mixtures. See Table III at end of this chapter. In the 
natural convection region, the results were mixed. Neither 
tube predominated consistently for pure R-114 and R-114/oil 
mixtures. Figures 6.20, 6.22 and 6.24 show similar comparisons 
on decreasing heat flux. It can be seen that the GEWA-K 40 fpi 
performs best in all three cases except at very low fluxes 
where results were less certain. However, the calculation for 
the heat flux did not take into account the actual increases 
in surface area for the finned tubes over the smooth tube. A 
fin root diameter was used to calculate surface area. The 
actual area increase is of the same magnitude as the heat 
transfer enhancements seen, indicating that both smooth and 
finned tubes give similar heat transfer performance. To 
demonstrate this, Figure 6.9, the GEWA-K 26 fpi tube in pure 
R-114, is replotted in Figure 6.81 using the actual fin 
surface area. The ratio of actual fin surface area to root 
diameter area was a factor of 4. Actual surface heat flux was 
calculated by dividing the root diameter surface heat flux by 
four. It can be seen in Figure 6.81 that in the natural 
convection region the GEWA-K 26 fpi tube performs comparably 
to the smooth tube. However, in the mixed boiling and nucleate 
boiling regions the GEWA-K 26 fpi tube shows an enhanced 
performance over the smooth tube up to three, which indicates 
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there are other mechanisms involved besides increased surface 



area . 

Figures 6.25, 6.26 and 6.27 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from a modified fin tube, 
GEWA-T 19 fpi . At a heat flux of 35 kW/m 2 , the heat transfer 
enhancement compared to the smooth tube for pure R-114, 
R-114/3% oil mixture and R-114/10% oil mixture was 3.1, 4.3 
and 5.2 respectively. Figures 6.28 and 6.29 show the 
comparative performance of the three refrigerant /oil mixtures 
on increasing and decreasing heat fluxes respectively. Similar 
to the GEWA-K type tube, only a slight degradation is observed 
in both the natural convection and nucleate boiling regions. 
Again at higher heat fluxes, the R-114/3% oil mixture 
performed best. 

Figures 6.30, 6.31 and 6.32 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from a GEWA-T 26 fpi tube. 
At a heat flux of 35 kW/m 2 , heat transfer enhancement compared 
to the smooth tube for pure R-114, R-114/3% oil mixture and 
R-114/10% oil mixture was 3.7, 4.6 and 4.8 respectively. 
Figures 6.33 and 6.34 show the comparative performance of the 
three refrigerant /oil mixtures on increasing and decreasing 
oil mixtures respectively. It can be seen that at higher heat 
fluxes, the R-114/3% oil mixture performed best and there was 
significant degradation in performance for the R-114/10% oil 



49 



mixture compared to the smooth tube. At lower heat fluxes, the 
pure R-114 outperformed the ref rigerant /oil mixtures. It can 
be surmised that the 'choking' effect predominates for the 
ref rigerant /oil mixtures. 

Figures 6.35, 6.37 and 6.39 show the comparative heat 
transfer performance of each ref rigerant /oil mixture boiling 
from the GEWA-T 19 and 26 fpi tubes for increasing heat 
fluxes. Figures 6.36, 6.38 and 6.40 show similar comparisons 
for decreasing heat fluxes. It can be seen that in the 
nucleate boiling region at the higher heat fluxes that the 
GEWA-T 26 fpi tube performed best. At the lower heat fluxes, 
in the natural convection region for an increasing heat flux 
run the tubes performed equally. For a decreasing heat flux 
run, the GEWA-T 19 fpi performed best. Recalling the channel 
structure of the GEWA-T tube suggests more complex phenomena 
than simple increase in surface area. Indeed Marto et al. 
[Ref. 16] have found that the gap width with GEWA-T tubes has 
a significant effect on boiling performance. Down to some 
limiting value, decreasing the gap width generally increases 
heat transfer coefficient. The enhancement can be attributed 
to strong interaction between heat transfer and hydrodynamic 
effects within GEWA-T channels. 

Figures 6.41, 6.42 and 6.43 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from another modified fin 
tube, GEWA-YX 26 fpi. At a heat flux of 35 kW/m 2 , heat 
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transfer enhancement compared to the smooth tube for pure 
R-114, R-114/3% oil mixture and R-114/10% oil mixture was 3.4, 
4.3 and 4.4 respectively. Figures 6,44 and 6.45 show 
comparative performance for the three refrigerant/oil mixtures 
on increasing and decreasing heat fluxes respectively. It can 
be seen that in the natural convection region slight 
degradation in performance occurs. At the higher fluxes in 
the nucleate boiling region, the R-114/3% oil mixture 
performed best as seen with other types of finned tubes. 

Figures 6.4 6, 6.4 8 and 6.50 show the comparative 
performance for each ref rigerant /oil mixture boiling from the 
GEWA-K, GEWA-T and GEWA-YX 26 fpi tubes for increasing heat 
fluxes. Figures 6.47, 6.49 and 6.51 show similar comparisons 
for decreasing heat fluxes. Both the modified finned tubes, 
GEWA-T and GEWA-YX, show marked increases in performance over 
the standard finned tube, GEWA-K, in the nucleate boiling 
region. Results were not so clear however in the natural 
convection region, possibly due to variation in surface area. 

3 . Boiling Performance of High Flux Tube 

Figures 6.52, 6.53 and 6.54 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from the High Flux tube. At 
a heat flux of 35 kW/m 2 , heat transfer enhancement compared to 
the smooth tube for pure R-114, R-114/3% oil mixture and 
R-114/10% oil mixture was 6.4, 6.2 and 3.8. The High Flux tube 
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displays a much larger relative temperature overshoot compared 
to the smooth tube and its point of incipience occurs at a 
much lower heat flux. Accuracy of the data was greatest at 
high heat flux and large superheat values (± 3.8%) and most 
uncertain at low heat flux and small superheat values 
(± 17.15%). See Uncertainty Analysis-Appendix E. It can be 
seen that at higher heat fluxes that the High Flux tube in the 
R-114/10 oil mixture performed almost equally to the smooth 
tube. At a heat flux of 100 kW/m 2 the heat transfer 
enhancement compared to the smooth tube was 1.1. Wanniarachchi 
et al. [Ref. 7] reported similar results. The marked decrease 
in boiling performance was attributed to the creation of an 
oil-rich layer within the boiling pores. The superior boiling 
performance for R-114/10% oil mixture compared to the smooth 
tube is completely recovered when the heat flux is decreased. 
It can be seen that in the natural convection region the 
performance of the High Flux tube is comparable at all oil 
concentrat ions to the smooth tube which indicates there is no 
effect of either tube type or oil concentration in this 
region. Figures 6.55 and 6.56 show the comparative performance 
of the three refrigerant /oil mixtures on increasing and 
decreasing heat fluxes respectively. The extreme rate of 
degradation of performance of the High Flux surface at 10% oil 
concentration with increasing heat flux is clearly seen as 
well as its complete recovery as the heat flux is decreased. 
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4 . Boiling Performance of the Thermoexcel Tubes 

Figures 6.57, 6.58 and 6.59 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from the Thermoexcel-E Tube. 
At a heat flux of 35 kW/m 2 , heat transfer enhancement compared 
to the smooth tube for pure R-114, R-114/3% oil mixture and 
R-114/10% oil mixture was 5.7, 5.3 and 5.2 respectively. 
Figures 6.60 and 6.61 show the comparative performance of the 
three refrigerant /oil mixtures on increasing and decreasing 
heat fluxes respectively. As seen with the High Flux tube, the 
Thermoexcel-E tube performed equally to the smooth tube in the 
natural convection region for all refrigerant /oil mixtures. 

Figures 6.62, 6.63 and 6.64 show the heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from the Thermoexcel-HE 
tube. At a heat flux of 35 kW/m 2 , heat transfer enhancement 
compared to the smooth tube for pure R-114, R-114/3% oil 
mixture and R-114/10% oil mixture was 4.8, 5.1 and 4.8. 
Similar performance to the Thermoexcel-E tube was observed. 
Figures 6.65 and 6.66 show the comparative performance of the 
three refrigerant /oil mixtures on increasing and decreasing 
heat flux respectively. At a higher heat flux of 100 kW/m 2 the 
heat transfer enhancements compared to the smooth tube for the 
R-114/3% oil mixture and R-114/10% oil mixture was reduced to 
2 . 2 and 1.9. 
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Figures 6.67, 6.69 and 6.71 show the comparative 
performance of each refrigerant /oil mixture boiling from both 
the Thermoexcel-E and HE tubes for increasing heat fluxes. 
Figures 6.68, 6.70 and 6.72 show similar comparisons on 
decreasing heat fluxes. As discussed earlier both tubes showed 
improved performance over the smooth tube in the nucleate 
boiling region. It can be seen that the Thermoexcel-HE tube 
performed marginally better than the Thermoexcel-E at high 
heat fluxes. In the natural-convection region both tubes 
performed comparably to the smooth tube. Wanniarachchi et al 
[Ref. 8] reported similar behavior for Thermoexcel-E and 
Thermoexcel-HE tubes. It can be surmised that in the natural 
convection region the small re-entrant cavities of both tubes 
do not increase the total boiling surface area like fins. The 
Thermoexcel-E and HE tubes 'appear' smooth and would be 
expected to perform comparably to the smooth tube which was 
indeed seen. 

5. Boiling Performance of the Turbo-B Tube 

Figures 6.73, 6.74 and 6.75 show heat transfer 
performance of pure R-114, R-114/3% oil mixture and R-114/10% 
oil mixture respectively, boiling from the Turbo-B tube. At a 
heat flux of 35 kW/m 2 , heat transfer enhancement compared to 
the smooth tube for pure R-114, R-114/oil mixture and 
R-114/10% oil mixture was 6.2, 6.7 and 5.6 respectively. As 
seen with all previous re-entrant tubes tested, comparable 
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performance with the smooth tube is observed in the natural 
convection region and marked improved performance over the 
smooth tube is observed in the nucleate boiling region. 
Figures 6.7 6 and 6.77 show the comparative performance of each 
refrigerant /oil mixtures boiling from the Turbo-B tube on 
increasing and decreasing heat fluxes respectively. Similar to 
the other re-entrant cavity surfaces tested, the R-114/10% oil 
mixture had the greatest rate of degradation for increasing 
heat flux in the nucleate boiling region. 

6 . Overall Performance of Smooth and Enhanced Surfaces 
Figures 6.78, 6.7 9 and 6.80 show the comparative 
performance of 8 of the 10 tube surfaces tested. The other 
two, Thermoexcel-E and GEWA-T 19 fpi, performed similarly to 
the Thermoexcel-HE and GEWA-T 26 fpi respectively, and so are 
not included. Three distinct groups of tubes can be identified 
based upon boiling heat transfer performance at all oil 
concentrations. These are: (1) smooth, (2) finned/modified 
finned and (3) re-entrant cavity. Within the finned/modified 
finned group, it can be seen that the GEWA-T is the best 
performer. Its enhanced performance was attributed to its 
channel geometry. The High Flux tube in R-114/10% oil mixture 
performed comparably to the smooth tube (heat transfer 
enhancement was 1.1) at a heat flux of 100 kW/m 2 . The severe 
degradation in performance for the High Flux tube can be 
attributed to the 'choking' of the porous surface. Table III 
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shows numerically the enhancement in performance over the 
smooth tube for all three ref rigerant /oil mixtures boiling 
from each of the tubes for three heat fluxes of 10, 35 and 100 
kw/m 2 . It can be seen at the moderate heat flux of 35 kW/m 2 , 
the finned/modified finned tubes show enhanced performances of 
2.4 to 5.2 over the smooth tube while the re-entrant cavity 
tubes show enhancements of 3.8 to 6.7 over the smooth tube. 
The effect of oil concentration is significant at high heat 
fluxes. At a heat flux of 100 kw/m 2 for pure R-114, the finned 
and re-entrant cavity tubes show comparable enhancement 
compared to the smooth tube. For R-114/10% oil mixtures at a 
heat flux of 100 kW/m 2 , the finned tubes performed best. 
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Figure 6.1 Repeatabi I ity Comparison For 



Q/CW/ m~2} 
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Figure 6.2 Performance Comparison foi 
Pure R-114 Boiling from a Smooth Tub' 
With Known Correlations 
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Figure 6,3 Performance Comparison F< 
Pure R-114 Boiling From Thermoexce I 
Effect of De-gassing Sur f ace/ F I u i d 
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Smooth Tube Surface 
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Figure 6.7 Performance Comparison for 
i ! ing R- 114/ 0%, 3% & 10% Oi I Mixtures 
Increasing Flux From Smooth Tube 
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Figure 6.20 Performance Comparison For 
Boiling Pure R-114 From 
GEWA-K 26/40 f p i Tubes 
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Figure 6.29 Performance Comparison For 
ot I i ng R- 1 14/ 0% J 3% & 10% 01 i Mixtures 
From GEWA-T 19 f p i Tube 
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Figure 6.45 Performance Comparison For 
oiling R- 1 14/ 0%, 3 % & 10% Oi I Mixtures 
From GEWA-YX 26 fpi Tube 
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Figure 6.53 Performance Comparison For 
Bo i I i ng R- 114/ 396 O i I Mixture From 
High Flux T ube 
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Figure 6.56 Performance Comparison For 
io i I i ng R- 1 14/ 0%, 3% & 10% Oil Mixtures 
Decreasing Flux From High Flux Tube 
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Figure 6 64 Performance Comparison For 
Boiling R- 1 14/ 10% 0 i I M i xt ure From 
Thermoexce I -HE Tube 
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Figure 6.80 Performance Comparison For 
Bo i I i ng R- 114/ 10% Oi I Mixture From 
Smooth arid Enhanced Tubes 
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TABLE III. HEAT TRANSFER COEFFICIENT ENHANCEMENT RATIOS FOR 
SMOOTH VERSUS ENHANCED SURFACES 



VII. CONCLUSIONS 



1. Accurate and repeatable boiling heat-transfer data 
using R-114 and R-114/oil mixtures was obtained on ten 
enhanced tube surfaces. 

2. Three distinct groups of boiling tubes can be 
identified by structure and boiling performance: 
smooth, finned and reentrant cavity. 

3. Significant increase in performance was observed for 
reentrant cavity surfaces versus either the finned or 
smooth surfaces within the nucleating regime at all 
heat fluxes tested. 

4. Degradation in boiling performance accompanied 
increases in oil concentration within the nucleating 
regime. Most notably, the high flux surface performed 
comparably to a smooth surface at the highest flux 
tested. 

5. Within the convection regime all tube surfaces tested 
performed equally well. No significant effect of oil 
concentration can be observed up to incipient 
nucleation . 

6. Distinct temperature overshoot within the convection 
regime occurs repeatably in all tube surfaces tested. 
The reentrant cavity surfaces exhibit the most extreme 
overshoot versus either the finned or smooth surfaces. 
Oil concentration has no significant effect. 
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VIII. RECOMMENDATIONS 



1. Due to the way a real flooded evaporator operates, it 
would be preferable to heat the tube within the single 
tube apparatus with hot water. 

2. As a complement to present single tube and bundle tube 
studies, the present single tube apparatus should be 
modified to accommodate two vertically aligned boiling 
tubes in order to study the effects of multiple tubes 
and variation of tube pitch to diameter ratios. 

3. The physical properties of the R-114/oil mixtures 
tested should be measured for comparison to the 
predictions used in the data reduction program and for 
better understanding of the mechanisms driving the 
degradation in heat transfer performance. 

4 . Neutrally buoyant particles should be placed in the 
pool to facilitate study of the circulation patterns 
within the flooded evaporator in more detail. 

5. A high speed camera should be used to study the 
nucleation process and circulation patterns in more 
detail . 
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APPENDIX A: THERMOPHYSICAL PROPERTIES OF R-114 



The following thermophysical properties of saturated R-114 
are plotted versus temperature °C: 



1. 


Liquid 


Density (kg/m 3 ) 


2. 


Liquid 


Dynamic Viscosity (Pa*s) 


3. 


Liquid 


Thermal Conductivity (W/m*K) 


4 . 


Liquid 


Specific Heat (J/kg*K) 


5. 


Latent 


Heat (J/kg) 
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Fluid Temperat ure/ C CD 
Figure A.1 Density of R-114 
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Fluid Temperature/CC} 

Figure A. 2 Dynamic Viscosity of R-114 



Thermal Conduct i v i ty/CW/m*IQ 
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Fluid Temperature/ C C} 

Figure A. 3 Thermal Conductivity of R-114 



Specific Heat/ CJ/ kg*K} 
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Fluid Temperature/CQ 
Figure A. 4 Specific Heat of R-114 



Latent Fleat/ C J/ kg} 
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Fluid Temperat ure/ C C} 
Figure A. 5 Latent Fleat of R- 1 14 



APPENDIX B: DATA ACQUISITION APPARATUS CALIBRATION 



A. BACKGROUND 

Early data runs revealed a disparity in calculated values 
between previously reported theses and those presently 
obtained for heat flux and heat-transfer coefficients. The 
thermocouple channels were checked for validity of temperature 
readings using other programs and data acquisition equipment. 
Their readings were found to be as accurate as possible 
(± 3%) . Voltage readings originating from the RMS converter 
were also checked out as accurate. The RMS converter converts 
true voltage measured into a proportional voltage output (0-10 
volts) to the data acquisition unit. The proportionality 
constant to convert RMS voltage readings. Vs, to real tube 
heater voltage, V, was equal to 25. This was verified using a 
standard voltmeter across the power leads of the cartridge 
heater measured over several power settings. The discrepancy 
was identified to be the measurement of the current supplied 
to the tube cartridge heater by the VARIAC unit as measured by 
the inductive pickup (Figure 3.13). The inductive measurement 
of current (in proportional volts) is no longer run through 
the RMS converter. Instead, it is run directly to the Data 
Acquisition Unit. A procedure to measure the true current and 
calculate the proportionality constant for converting the 
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inductive current. Is, to the true current, I, was established 



as follows: 



B . EQUIPMENT USED 



1 . Voltmeter 

2. 2 ohm resistor 

3. HP-3852A Data Acquisition Unit 

4. Computer - HP-9300 Series 

5. Current Sensor (Inductive Pickup) 

6. Voltage Sensor/RMS Converter 

7. Cartridge Heater - lkW, 220V rating 

8. VAR I AC - power supply 220V AC 

9. SETUP8 program 



C . CALIBRATION PROCEDURE 



1. Ensure all power is turned off to the VARIAC. 

2. Insert the standard 2 ohm resistor in series with the 
cartridge heater as shown in Figure C.l. 

3. Connect leads of voltmeter across the 2 ohm resistor to 
measure the drop in voltage. At any given power 
setting, the current through the resistor will equal 
the current through the cartridge heater. This assumes 
the resistances of both heater and resistor are 
constant. This was verified as shown in Figure C.2. The 
total resistance of the heater/resistor rapidly and 
asymptotically approached a constant value. 

4. Connect power to the VARIAC. 

5. Load and run SETUP8 program to provide the voltage and 
current readings made by the data acquisition unit. 
Note the SETUP8 program used in calibration differed 
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from that listed in Appendix G. A proportionality 
constant of 2.0 was assumed for calibration then 
modified once the true value was calculated. 

6. Starting with the VARIAC set at zero volts, 
incrementally increase the voltage up to the maximum of 
220V. 

7. At each voltage, measure the drop in voltage across the 
resistor. The current (amps) is calculated: 

I = Drop in voltage (volts) / 2 ohms 

8. Note and compare the current provided by the program. 

9. At the maximum 220V, incrementally decrease voltage 
down to 0, taking comparative current readings as 
before, in order to detect any hysteresis effect. None 
was observed. 

10. Plot current (amps) vs. voltage (volts) . See Figure 
C.3. Note the linearity of both data sets. The true 
meter reading is some factor smaller than the data 
acquisition readings. 

11. The calibration constant, C, is calculated by: 

C = (Z Meter current) * 2 

(I Data acquisition current) 

C = 1.9182 



D. DISCUSSION 

A disparity between the previous reported theses values 
and those presently obtained still exists. Calculated values 
reported in previous theses for the heat flux and the heat 
transfer coefficient are consistently higher. No previous 
thesis reports such a calibration as carried out above and it 
is thought that this is the cause of the discrepancy. 
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Figure B.2 Cartridge Heater Reslstence 
For High Flux Tube 
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Figure 83 Var lac Current Calibration 
Using High Flux Tube 



APPENDIX 



C: AN EXAMPLE OF REPRESENTATIVE DATA RUN 



Honth r deta and tine : 5 Jul 1991 

NDTE : Propran nane : DRP8 
Disk nunber ■ 01 
New file nana: DATA07051S 
TC ib defective at location 
Tube Nunber: 5 



Deta Set Nunber 
TC No : 1 

Tenp : 5.11 

Twa TliQd 
5.08 2.3G 

Thatab Htube 



1 Bulk Dll X - 
2 3 4 

0.00 5.15 5.01 

Tliod2 Tvapr Paet 
2.2E 3.91 - 1 . E5 



Qdo 



2.770 I . 590E+02 4.40EE+02 



Data Sat Nunbar 
TC No: 1 

Tanp : 4.94 

Twa TliQd 
4.95 2.23 

Thatab Htube 



• 2 Bulk Oil X - 

2 3 4 

0.00 5.04 4.87 

TliQd2 Tvapr Peat 
2.18 3.29 -1.7E 

Qdp 



2.740 1 . E03E+02 4.391E+02 



Data Sat Nunbar ■ 3 8ulk Oi 1 X 



TC No: I 

Tanp : E.08 

Twa Tlipd 
E . 07 2.38 

Thatab Htuba 



2 3 4 

0.00 G . 1 9 5.9G 

TliQd2 Tvapr Paat 
2. IE 3.27 -1.E9 

Qdp 



3.80G 1 . E94E+02 G.44GE+02 



Data Sat Nunbar - 4 8u)k Oil X - 

TC No: 1 2 3 4 

Tanp : G . 12 0.00 E.23 E.02 

Twa TliQd Tliod2 Tvapr Paat 

E. 1 2 2.39 2.19 3.29 -1.E7 

Thatab Htuba Qdp 

3.828 1 . E82E+02 E.439E+02 



Data Sat Nunbar 
TC No : 1 

Tanp : 7.13 

Twa TliQd 
7.13 2.24 

Thatab Htuba 



5 Bulk Oil X • 

2 3 4 

0.00 7.27 E.97 

Tlipd2 Tvapr Paat 
2. IE 2.92 -1.77 

Qdp 



4.934 1.791E+02 8.834E+02 



Data Sat Nunbar - 8 Bulk Oil 2 



TC No: 1 

Tanp : 7.09 

Twe Tliod 
7.11 2.20 

Thatab Htube 



2 3 4 

0.00 7.25 G . 95 

Tliqd2 Tvapr Paat 
2.12 2.91 -1.81 

Qdp 



4.945 I.78GE+02 8.834E+02 



Data Sat Nunbar - 7 Bulk Oil X - 

TC No: 1 2 3 4 

Tanp ; 9.5E 0.00 9.74 9.3E 

Twa TliQd TliQd2 Tvapr Paat 

9.57 2.47 2.22 2.94 -1 .El 



Thatab Htuba 



Qdp 



7.223 1 . 8E3E+02 1.34GE+03 



17:21 : 5 1 



0.0 2.11545538805 

5 E 7 8 

5.12 5.08 5.09 5.11 

Tauno 
-13.8 



0.0 2.11545539033 

5 G 7 8 

4.99 4.95 4.9E 4.98 

T aunp 
-13.7 



0.0 2.11 54553975E 

5 E 7 8 

G . 1 4 G .09 E.ll G . 15 

Taunp 
-13. E 



0.0 2.11545539780 

5 E 7 8 

G . 1 8 E.14 E . 1 5 E . 1 8 

Taunp 
-13. E 



0.0 2 . 1 1 545540G28 

5 E 7 8 

7.24 7.15 7.17 7.25 

Taunp 
-13. E 



0.0 2.1 1545540E55 

5 6 7 8 

7.21 7.13 7.15 7.20 

Taunp 
-13.5 



0.0 2.11545541023 

5 E 7 8 

9.71 9.E0 9.E2 9.75 

Taunp 
-13.5 
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Date Sat 


Number 


- 8 8ulk Oil X - 


0.0 2. 


11545541165 


TC No: 


1 


2 3 4 


5 6 


7 


8 


Temp : 


9.33 


0.00 9.52 9.11 


9.49 9.37 


9.39 


9.50 


Two 


Tliqd 


Tliod2 Tvapr Paat 


Taump 






9.33 


2.21 


2.13 2.82 -1.80 


-13.5 






Thatab 


Htube 


Qdo 








7. 1 65 


1 . 874E+02 1.3436+03 








Oata Set 


Number 


- 9 Bulk Oil X • 


0.0 2. 


1 1 54554 1 546 


TC No: 


1 


2 3 4 


5 6 


7 


8 


T amp : 


10.99 


0.00 11.53 11.09 


11.50 11.38 


1 1 .40 


11.40 


Two 


T1 lod 


Tliod2 Tveor Paat 


Taump 






H.2C 


2.30 


2.16 2.66 -1.73 


-13.5 






The tab 


Htube 


Odp 








S.02B 


1 . 9636+02 1 . 7726+03 








Data Set 


Number 


• 10 Sulk Oil X - 


0.0 2. 


11545541589 


TC No: 


1 


2 3 4 


5 6 


7 


8 


Tamp : 


11.00 


0.00 11.52 11.09 


11.50 11.37 


11.40 


1 1.41 


Twa 


T1 iqd 


Tliod2 Tvaor Paet 


Taump 






>1.26 


2.32 


2.22 2.67 -1.69 


-13.5 






That ab 


Htube 


Odp 








6.988 


1.977E+02 1 .7776+03 








Data Set 


Number 


- 11 8ulk Oil X - 


0.0 2. 


1 1545541889 


TC No: 


1 


2 3 4 


5 6 


7 


8 


Temp : 


11.45 


0.00 12.32 11.93 


12.30 12.18 


12.21 


12. 10 


Twa 


T 1 iqd 


Tliod2 Tvaor Paat 


Taump 






11.99 


2.24 


2.17 2.67 -1.76 


-13.5 






That ab 


Htube 


Odp 








9.783 


2.074E+02 2.029E+03 








Data Set 


Number 


- 12 8ul k Oil X • 


0.0 2. 


11545541913 


TC No: 


1 


2 3 4 


5 6 


7 


8 


Tamp : 


1 1.50 


0.00 12.38 11.97 


12.35 12.24 


12.25 


12.14 


Twa 


Tliod 


Tliod2 Tvapr Paat 


Taump 






12.04 


2.26 


2.19 2.69 -1.74 


-13.5 






Thatab 


Htube 


Odp 








9.814 


2.0GSE+02 2.027E+03 








Oata Set 


Number 


• 13 Sulk Oil X - 


0.0 2. 


1)545542290 


TC No: 


1 


2 3 4 


5 6 


7 


8 


Tamp : 


12.08 


0.00 13.11 12.71 


13.09 12.98 


12.99 


12.82 


Twa 


Tliqd 


Tliqd2 Tvapr Paat 


Taump 






12.74 


2.22 


2.20 2.96 -1.75 


-13.5 






Thetab 


Htube 


Odp 








10.527 


2.1126+02 2.223E+03 








Data Sat 


Number 


» 14 Sulk Oil X - 


0.0 2. 


11545542315 


TC No: 


1 


2 3 4 


5 6 


7 


8 


Tamp : 


12.02 


0.00 13.05 12.67 


13.03 12.92 


12.94 


12.75 


Twa 


Tliqd 


Tliqd2 Tvapr Paat 


Taump 






12.68 


2.20 


2.17 2.96 -1.78 


-13.5 






Thatab 


Htube 


Odp 








10.494 


2.11 0E+02 2. 2146+03 








Data Sat 


Number 


- 15 Bulk Oil X - 


0.0 2. 


11545542511 


TC No: 


1 


2 3 4 


5 6 


7 


8 


Temp : 


12.55 


0.00 13.83 13.53 


13.82 13.72 


13.73 


13.43 


Twa 


Tliqd 


Tliqd2 Tvapr Paat 


Taump 






13.42 


2.19 


2.16 2.92 -1.79 


-13.5 






Thatab 


Htube 


Odp 








1 1 .239 


2 . 1 82E+02 2.452E+03 
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Oat a Set 


Number 


- 18 Bui ► Oil X - 


0.0 


2. 


1 1S4SS42 


534 


TC No: 


1 


2 3 4 


S 


8 


7 


8 


Temp : 


12. E2 


0.00 13.85 13.81 


13.83 ! 


13. 7S 


13. 7S 


13.47 


Twa 


T1 iQd 


Tliod2 Tveor Peat 


Tftumo 








13. 4G 


2. 28 


2.22 2.91 -1.72 


-13. S 








Thet ab 


Htube 


Odo 










1 1.217 


2.I83E+02 2 . 449E+03 










Oata Set 


Number 


- 17 Bulk Oil X - 


0.0 


2. 


1 1 545543 


206 


TC No: 


1 


2 3 4 


S 


6 


7 


8 


Temo : 


2.89 


0.00 2.74 2.89 


2.69 


2.S9 


2 .57 


2.68 


Twa 


Tljpd 


Tl jod2 Tvapr Peat 


T aump 








2.S8 


2.44 


2.19 2.01 -1.64 


-13. S 








The tab 


Htube 


Odo 










.270 


1.128E+04 3. 036E+03 










Data Set 


Number 


- 18 Bulk Oil X - 


0.0 


2. 


11545543229 


TC No: 


1 


2 3 4 


S 


6 


7 


8 


Tamo : 


2.88 


0.00 2.71 2.88 


2.67 


2. 56 


2.55 


2.65 


Twa 


Tlipd 


Tlipd2 Tvapr Paat 


Taump 








2. SB 


2.39 


2. IS 1.98 -1.B9 


-13. S 








The tab 


Htuba 


Odp 










.288 


l . 0S2E+04 3.029E+03 










Data Sat 


Number 


- 19 Bulk Oil X - 


0.0 


2. 


11545543435 


TC No: 


1 


2 3 4 


S 


6 


7 


8 


Tamp : 


3. IS 


0.00 2.97 2.92 


2.92 


2.79 


2.75 


2.91 


Twa 


Tliod 


Tlipd2 Tvapr Paat 


T aumo 








2.7B 


2.41 


2.17 2.09 -1.67 


-13. 5 








Thetab 


Htuba 


Odp 










.472 


9. 576E+03 4. S23E+03 










Data Set 


Number 


- 20 Bulk Oil X - 


0.0 


2. 


11545543508 


TC No: 


1 


2 3 4 


S 


6 


7 


8 


Tamp : 


3.12 


0.00 2.93 2.89 


2.89 


2 . 75 


2.72 


2.88 


Twa 


T] iqd 


Tliod2 Tvapr Paat 


Taump 








2.73 


2.47 


2.19 2.08 -1.62 


-13. S 








Thetab 


Htube 


Odp 










.394 


1 . 1 44E+04 4 .506E+03 










Data Sat 


Number 


- 21 Bulk Oil X - 


0.0 


2. 


11545543579 


TC No: 


1 


2 3 4 


S 


6 


7 


8 


Temp : 


3.02 


0.00 2.83 2.79 


2.79 


2.64 


2.62 


2.77 


Twa 


Tljpd 


Tliod2 Tvapr Paat 


Taump 








2.B2 


2.37 


2.09 1.98 -1.74 


-13. S 








Thetab 


Htuba 


Odp 










.394 


1.149E+04 4 . S2SE+03 










Data Set 


Number 


- 22 Bulk Oil X - 


0.0 


2. 


11545543829 


TC No: 


1 


2 3 4 


S 


6 


7 


8 


Temp : 


3 . 45 


0.00 3.20 3.16 


3. IS 


2.96 


2.91 


3.14 


Twe 


Tl iqd 


Tlipd2 Tvapr Paat 


T aump 








2.B9 


2.41 


2.16 2.76 -1.67 


-13. S 








Thetab 


Htube 


Odp 










. B03 


1 . 1 78E+04 7. 099E+03 










Data Sat 


Number 


- 23 Bulk Oil X - 


0.0 


2. 


11545543855 


TC No: 


1 


2 3 4 


S 


6 


7 


8 


Temp : 


3.44 


0.00 3.19 3. IS 


3.14 


2.96 


2.90 


3.13 


Twa 


Tliod 


Tllod2 Tvapr Paat 


Taumo 








2.88 


2.39 


2. IS 2.87 -1.69 


-13. S 








Thetab 


Htube 


Odp 










.812 


1.1S8E+04 7. 092E+03 
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Data Set Number 
TC No: 1 

Temp : 3.95 

Twa Tliod 
3.J2 2.3E 



- 24 Bull. Oil X • 

2 ? 4 

0.00 3.56 3.55 

Tliod2 Tvaor Peat 
2.17 2.71 -1.69 



Thatab Htuba Ode 

.857 1 . 279E+04 1 .096E+04 



Data Set Number » 25 Bull- Dll 1 • 

TC No: 1 2 3 4 

Tame : 3.95 0.00 3.56 3.55 

Twa TJiod TliodC Tvaor Paat 

3.12 2.36 2.16 2.75 -1.70 

Thetab Htuba Ddo 

.859 1.279E+04 1.098E+04 



Data Set 
TC No: 
Temp : 
Twa 
3.36 
Thatab 
1 .073 



Number • 26 Bulk DiJ 1 - 

12 3 4 

4.48 0.00 3.92 3.95 

Tliod TJiod2 Tvapr Peat 

2.38 2.19 2.07 -1.67 

Htube Ddo 

I . 4 1 1 E +04 1 . 5 1 4E +04 



Data Set Number - 27 Bulk Dil 2 • 

TC No: 1 2 3 4 

Tamp : 4.45 0.00 3.B9 3.93 

Twa Tliod Tliod2 Tvapr Peat 

3.33 2.33 2.15 2.04 -1.72 

Thetab Htuba Ddo 

1.092 1 . 385E 4-04 J.512E+04 



Data Set Number • 28 Bulk Ol 1 X - 

12 3 4 

5.12 0.00 4.37 4.43 

Tliod Tliqd2 Tvaor Peat 

2.35 2.17 2.08 -1.70 

Htuba Ddo 



TC No.- 
Temp : 
Twa 
3.65 
Thatab 
1 .390 



1.437E+04 1 . 996E+04 



Data Sat Number 
TC No : 1 

Tamp : 5.12 

Twa Tliod 
3.64 2.34 

Thetab Htuba 



29 Bulk Oil 1 - 
2 3 4 

0.00 4.36 4.42 

Tliqd2 Tvapr Peat 
2.16 2.07 -1.71 

Ddp 



1.393 J.433E+04 1 . 996E+04 



Data Set Number 
TC No : 1 

Tamp : 6.47 

Twa Tliod 
4.23 2.31 

Thatab Htuba 



30 Bulk Oil 1 - 
2 3 4 

0.00 5.25 5.39 

Tliod2 Tvaor Peat 
2.18 2.10 -1.72 

Ddp 



1.989 J.S0SE+04 2.993E+04 



Data Sat Number 
TC No : 1 

Tamo : 6.45 

Twa Tliod 
4.22 2.27 



- 31 Sulk Oil 1 - 

2 3 4 

0.00 5.23 5.37 

Tliod2 Tvapr Peat 
2.15 2.08 -1.75 



Thetab Htube Ddp 

2.003 1.497E+04 2.998E+04 



0.0 

5 

3.50 

T eump 

-13.6 



0.0 

5 

3.49 
T eump 
-13.6 



0.0 

5 

3.85 

Teump 

-13.6 



0.0 

5 

3.83 
T eump 
-13.6 



0.0 

5 

4.27 
T eumo 
-13.6 



0.0 

5 

4.26 
T eump 
-13.6 



0.0 

5 

S.J 1 
T eump 
-13.6 



0.0 

5 

5.09 
T eumo 
-13.6 



2.11545544709 
6 7 8 

3.29 3.20 3.52 



2.1 1545544732 
6 7 8 

3.28 3.19 3.52 



2.1 1545545094 
6 7 8 

3.62 3.49 3.92 



2.1 1545545119 
6 7 8 

3.60 3.46 3.89 



2. 1 1545545460 
6 7 8 

4.04 3.86 4.40 



2.11545545484 
6 7 8 

4.03 3.84 4.39 



2. 1 )545546220 
6 7 8 

4.90 4.58 5.33 



2. 11545546248 
6 7 8 

4.88 4.57 5.31 
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Data Sat 


Nunbar 


- 32 Bulk Oil X - 


0.0 


2. 


11545546659 


TC No: 


1 


2 3 4 


5 


6 


7 8 


T ano : 


8.99 


0.00 G.GG 6.92 


6.39 


6.17 


5.57 6.63 


Twa 


Tliqd 


Tliod2 Tvapr Peat 


T eunp 






5. 16 


2.33 


2.18 2.13 -1.71 


-13.4 






Thatab 


Htuba 


Odp 








2.910 


1 .SS9E+04 4 . S3SE+04 








Data Sat 


Nunbar 


- 33 8ul k Oil X - 


0.0 


2. 


11545546682 


TC No: 


l 


2 3 4 


5 


6 


7 8 


T amp : 


8.97 


0.00 6. 64 6.90 


6.37 


6.)4 


5.56 6.62 


Twa 


Ti lod 


T 1 iod2 Tvapr Paat 


Taunp 






5. IS 


2.30 


2.17 2.10 -1.73 


-13.3 






Thatab 


Htuba 


Odo 








2.911 


1 . SS6E+04 4.S31E+04 








Oata Sat 


Nunbar 


- 34 Bulk Oil X - 


0.0 


2. 


11545546951 


TC No: 


1 


2 3 4 


5 


6 


7 8 


Tanp : 


12. SS 


0.00 8.98 9.39 


8.54 


8.06 


6.96 8.56 


Twa 


Tliod 


Tliod2 Tvapr Paat 


T eunp 






6. S3 


2.29 


2.16 2.15 -1.74 


-13.2 






Thatab 


Htuba 


Odp 








4.305 


1 . G32E+04 7.028E+04 








Data Sat 


Nunbar 


- 35 Bulk Oil X - 


0.0 


2. 


11545546986 


TC No: 


1 


2 3 4 


5 


6 


7 8 


Tanp : 


12.59 


0.00 9.01 9.42 


8.57 


8.09 


6.99 8.60 


Twa 


T 1 iod 


Tliod2 Tvapr Paat 


T aunp 






6.56 


2.31 


2.19 2.18 -1.71 


-13.2 






Thatab 


Htuba 


Odp 








4.315 


1 .627E+04 7.019E+04 








Data Sat 


Nunbar 


- 36 Bulk Oil X - 


0.0 


2. 


11545547257 


TC No: 


l 


2 3 4 


5 


6 


7 8 


Tanp : 


16 .06 


0.00 11.35 11.79 


10.57 1 


10.04 


8.26 10.43 


Twa 


Tliqd 


Tliqd2 Tvapr Paat 


Taunp 






7.84 


2.23 


2.15 2.14 -1.77 


-12.9 






Thatab 


Htuba 


Odp 








S.849 


1.G94E+04 9. 5G9E+04 








Data Sat 


Nunbar 


- 37 Bulk Oil X - 


0.0 


2. 


11545547281 


TC No: 


1 


2 3 4 


5 


6 


7 8 


Tanp : 


1 G . 1 1 


0.00 11.39 11.81 


10.59 1 


10.07 


8.28 10.46 


Twa 


Tliqd 


Tliqd2 Tvapr Paat 


Taunp 






7.87 


2 . 26 


2.18 2.17 -1.75 


-12.8 






Thatab 


Htuba 


Odo 








5. 6 56 


1.G92E+04 9.S72E+04 








NOTE: 37 


data runs wara atorad In fila OATA070SIS 


NOTE: 37 


X-Y paira wara atorad in plot data 


fila 


PLOT0705I5 
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APPENDIX D: SAMPLE CALCULATION 



Data run number 9 (saturation temperature was 2.22 °C and 
heat flux was 20.14 kW/m 2 ) was chosen for the sample 
calculation . 



A. TEST-SECTION DIMENSIONS 



Do = 


0.01588 


m 


Di = 


0.01270 


m 


D1 = 


0.01143 


m 


D2 = 


0.01588 


m 


L 


0.20320 


m 


Lu = 


0.07620 


m 



B. 



MEASURED PARAMETERS 



V 


= 112.08 


volts 


I 


= 1.94 


volts 


T1 


= 16.51 


°C 


T2 


= 15.80 


•c 


T3 


= 15.88 


*c 


T4 


= 15.93 


•C 


T5 


not read due to defective thermocouple 


T6 


= 15.94 


*C 


T7 


= 15.63 


*c 


T8 


= 16.16 


•c 


Tsat 


2.22 


*c 


kc 


= 344.00 


(W/m- K) 


OUTER 


WALL TEMPERATURE OF THE BOILING TUBE 


P 


ti-Do = 3.1416 • 0.01588 = 0.0499 (m) 


Ac = 


7t (Do 2 — Di 2 ) / 4 = 3.1416 • [ (0 . 01588) 2 - (0 . 0127) 2 ] 


Ac = 


7.14E-5 (m 2 ) 


Qh = 


VI = 112.08 • 1.94 = 217.44 (W) 


Tavg 


= £ Tn /n 


. = 15.97 CC) 


Two 


= Tavg - 


Qh- [ In (D2/D1) / (2- 7i- L- kc) ] 


Two 


= 15.97 - 


217.44- 




[In (15. 


88/11.43) / (2- 3.1416- 0.203- 344) ] 



] / 4 
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Two 

e 



= 15.80 CC) 

= Two - Tsat = 15.80 - 2.22 = 13.58 °C 



D. PROPERTIES OF R-114 AT FILM TEMPERATURE 

Tf = (Two+Tsat) /2 = (15 .80 + 2 .22) /2 = 9.01 *C 
p = exp [ -4 .4636 + (1011 . 47/Tf ) ] • IE-3 

p = exp(-4.4636+(1011.47/9.01) ] • IE-3 = 4.16E-4 

(N* s/m 2 ) 

Tc = Critical Temperature (R) = 753.95 (R) 

Tf = 9.01 CC) = 507.89 (R) 

j = 1 - Tf(R)/Tc(R) = 1 - 507.89/753.95 = 0.33 

p = 581.77 + 984 .15- j (1/3) + 263.02*j + 279.99*j 1/2 + 

17.94* j 2 (kg/m 3 ) 
p = 1504.22 (kg/m 3 ) 

v = p/p = 4 . 16E-4/1504 .22 = 2.77E-7 
k = 7.10E-2 - (2.61E-4 • Tf(*C)) 

k = 7.10E-2 - (2.61E-4 • 9.01) = 6.86E-2 (W/m*K) 

Tf = 9.01 CC) = 282.16 (K) 

Cp = 400 + (1.65 * Tf) + (1.51E-3 * Tf 2 ) - (6.68E-7 * 

Tf 3 ) 

Cp = 970.43 ( J/kg* K) 

a = k / p • Cp = 6. 86E-2/ (1504.22* 970.43) = 4.7E-8 (m 2 /s) 

P = - (Ap/AT) /p = 1.92E-3 (1/K) 

Pr = v /a = 5.88 



E. HEAT -FLUX CALCULATION 

Average natural-convection heat transfer coefficient at 
non-boiling ends: 



h = k_ • 0.60+0.387' 

Do 



( g* P* (Do 3 ) * 0* 
{ v* a* ] 



Tanh (m* Lu) 



\ 1/6 



Lu* m 



[l+(0.559/Pr) 9/16 ] 8/27 
m = [h* p/ (kc* Ac) ] 1/2 



h = 235.4 (W/m 2 * K) 



The value for h was computed iteratively. A value for h of 190 
(W/m 2 *K) was assumed at start of iteration. 

Heat-transfer rate through non-boiling ends: 

m = 21.87 (1/m) 

Qf = (h* p* kc* Ac) 1/2 * 0* Tanh (m* Lu) = 6.79 (W) 
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Heat flux through active boiling surface 

Q = Qh - 2 • Qf = 217.44 - 2 • 6.79 = 203.85 (W) 

Ab = it-Do-L = 3.1416 • 0.015875 * 0.2032 = 1.01E-2 (m 2 ) 
q" = Q/Ab = 203 . 85/1 . 02E-2 = 20, 108 (W/m 2 ) 

h = q/6 = 20,108/13.59 = 1480 (W/m 2 - K) 

The following are the results obtained from the computer 
by running the data reduction program DRP8 (See Appendix H) . 

q" = 20, 140 (W/m 2 ) 

0 = 13.64 CC) 

h = 1477 (W/m 2 - K) 
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TABLE IV. DIMENSIONS OF BOILING TUBES 





Dl 

(mm) 


D2 

(mm) 


Di 

(mm) 


Do 

(mm) 


L 

(mm) 


Lu 

(mm) 


Kcu 
(W) 
m 2 - K 


Smooth 


12.44 


15.88 


12.70 


15.88 


203.2 


76.20 


344 


GEWA-K 
26 fpi 


10.1 


12.7 


11.11 


12.7 


203.2 


76.2 


398 


GEWA-K 
40 fpi 


10.1 


12.7 


11 . 11 


12.7 


203.2 


76.2 


398 


GEWA-T 
19 fpi 


11.57 


13.8 


11.8 


13.31 


203.2 


76.2 


398 


GEWA-T 
2 6 fpi 


11.57 


13.8 


11.8 


13.31 


203.2 


76.2 


398 


GEWA-YX 
2 6 fpi 


11.57 


13.8 


11.8 


13.31 


203.2 


76.2 


398 


Thermoexcel- 

E 


11.57 


13.8 


11.8 


13.31 


203.2 


76.2 


398 


Thermoexcel- 

HE 


11.57 


13.8 


11.8 


13.31 


203.2 


76.2 


398 


Turbo-B 


11.57 


13.8 


11.8 


15.8 


203.2 


76.2 


398 


High Flux 
(Copper/ 
Nickel 
95/5) 


12.95 


15.82 


13.2 


15.82 


203.2 


76.2 


45 
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APPENDIX E: UNCERTAINTY ANALYSES 



Four data points from two data runs were chosen for the 
uncertainty analysis: data run number 8 (smooth tube, 0% oil, 
100, 000 & 3500 W/m 2 ) and data run number 74 (Thermoexcel-E, 0% 
oil, 6000 & 100,000 W/m 2 . The dimensions of the boiling tubes 
can be found in Appendix D, Table IV. The measured and 
calculated parameters used in this analysis were obtained as 
shown in Appendix D and are listed in Table V of this 
appendix. All uncertainties are presented as a percentage of 
the calculated parameter. What follows is a sample calculation 
for data run number 8: smooth tube, 0% oil at a heat flux of 
93,130 W/m 2 using the methods outlined by Kline and McClintock 
[Ref. 17], All other data point uncertainties were calculated 

similarly, the results of which are listed in Table V of this 

appendix . 

A. UNCERTAINTY IN SOURCE HEAT-TRANSFER RATE 

Qh = VI W 

Is = 2.122 volts 61 = ± 0.025 amp 

Vs = 9.45 volts 6V = ± 0.05 volts 

I = 1.91 82 -Is = 4.07 amp 

V = 25-Vs = 236.26 volts 

where : 

6 = uncertainty in measurement and calculation 

6Qh/Qh = ((6V/Vs) 2 + (6l/Is) 2 ) 1/2 

SQh/Qh = ( (0.05/9. 45) 2 + (0 . 025/2 . 122) 2 ) K2 

6Qh/Qh = 1.29 percent 
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B. UNCERTAINTY IN SURFACE AREA 



Ab = %• Do* L 

Do - 15.88 (mm) ADo = 0.1 (mm) 

L = 203.20 (mm) AL = 0.1 (mm) 

AAb/Ab = ( (ADo/Do) 2 + (AL/L) 2 ) 1/2 

AAb/Ab = ( (0 .1/15.88) 2 + (0 . 1/203 .2) 2 ) 1/2 

AAb/Ab = 0.63 percent 



C. UNCERTAINTY IN WALL SUPERHEAT 
AT = Two ~ Tsat 

Tsat = 2.22 *C AT = 0.01 ‘C 

Two = Tavg - Qh [ (In (D2/D1) ) / (2* it* L* kc) ] 

Tavg = (ETn/8) where n = 1 to 8 

Tn = thermocouple readings 

T1 = 23.86 ‘C T2 = 19.62 °C T3= 20.38 *C 

T4 = 19.87 *C T5 = Defective thermocouple 

T6 = 20.18 *C T7 = 20.16 *C T8 = 20.20 *C 

Tavg = 20.18 *C 

S.D. = ( (E (Tn - Tavg) 2 ) /7) 1/2 = 0.19 *C 

where S.D. = standard deviation 



The logarithmic term in equation of Two is very small 
compared to the standard deviation, this term can be neglected 
for the uncertainty analysis. 



Two = Tavg = 18.72 *C ATwo = S.D. = 0.19 ’C 

AT = 17.96 *C 



AAT/AT = ((ATwo/AT) 2 + (-ATsat /AT) 2 ) 1/2 
AAT/AT = ( (0.93/16. 5) 2 + (0. 01/16. 5) 2 ) 1/2 

AAT/AT = 1.06 percent 



D. UNCERTAINTY IN HEAT FLUX 

q = (Qh-2Qf ) /Ab 

Qh = 961.58 W AQh = 12.42 W 
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Assuming the same proportion in the uncertainty for Qf : 



Qf = 8.99 W 6Qf = 0.12 W 

Qh - 2Qf = 943.6 W 



6q/q 


= [ (6Qh/ (Qh-2Qf) ) 2 + (26Qf/ (Qh-2Qf) ) 2 ) + 

(6Ab/Ab) 2 ) ] 1/2 


6q/q 

6q/q 


= [ (12. 42/943. 6) 2 + (0 . 24/943 . 6) 2 + (0 . 0063) 2 ] 1/2 

= 1.46 percent 



E. UNCERTAINTY IN BOILING HEAT-TRANSFER COEFFICIENT 



h = 


q/AT 


6h/h 

6h/h 

6h/h 


= [ (6q/q) 2 + (6AT/AT) 2 ] 1/2 

= [ (0 . 014 6) 2 + (0 . 01 . 06) 2 ] 1/2 

= 1.8 percent 
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TABLE V. UNCERTAINTY ANALYSIS OF FOUR DATA POINTS 



PARAMETERS 


SMOOTH 
100 kW/m 2 


SMOOTH 
3500 W/m 2 


THERMO-E 
6000 W/m 2 


THERMO-E 
100 kW/m 2 


Is (amp) 


2 . 122 


0.454 


0.494 


2.06 


Vs (volts) 


0.05 


2.09 


2.20 


8.88 


I (amp) 


4.07 


0 . 90 


0.98 


3.95 


6 I (amp) 


0.025 


0.025 


0.025 


0 . 025 


V (volts) 


236.26 


52.29 


55.28 


222.05 


8V (volts) 


0.05 


0.05 


0.05 


0.05 


Qh (W) 


961.58 


47.06 


54.17 


877.1 


5Qh (W) 


12.42 


2.83 


3.00 


7.43 


5Qh/Qh (%) 


1.29 


6.00 


5.55 


0.95 


Do (mm) 


15.88 


15.88 


13.31 


13.31 


o 

Q 


0.1 


0.1 


0.1 


0.1 


L (mm) 


203.2 


203.2 


203.2 


203.2 


8l (mm) 


0 . 1 


0 . 1 


0.1 


0.1 


Ab (m 2 ) 


1 . 014E-2 


1 . 014E-2 


8 . 50E-3 


8 . 50E-3 


6Ab/Ab (%) 


0.63 


0.63 


0.75 


0.75 


6Tsat CC) 


0.01 


6.0# 


0.01 


2.91 


Two CO 


20.18 


12.92 


2.77 


11.13 


6Two CO 


0.19 


0.28 


8.53E-2 


0.36 


AT CO 


17.96 


10.70 


0.53 


9.80 


6AT/AT (%) 


1.06 


2.62 


16.20 


3.67 


Qf (W) 


8.99 


5.02 


5.69 


2.91 


6Qf 


0.12 


0.03 


0.15 


2 . 4 8E-2 


Qh-2Qf (W) 


943.60 


37.02 


54.06 


871.3 


q (W/m 2 ) 


93, 081 


3664 


6109 


102,542 


5q/q (%) 


1.46 


7 . 67 


5.63 


1.14 


h (W/m 2 «K) 


5182 


342.4 


11,200 


10, 463 


8h/h (%) 


1.80 


8.10 


17.15 


3.84 
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APPENDIX F: SETUP PROGRAM 



The following is a listing of the complete setup program 
used in the preparation of the apparatus for testing. It 
consists of five sections: 

1 . Monitor the Sump temperature 

2. Monitor the evaporator liquid temperature 

3. Measurement and readout of all thermocouple channels. 

4 . Measurement and readout of the power supplied to the 
tube cartridge heater. 

5. Measurement and readout of the power supplied to the 
auxilliary heaters if used. 

The program was written in Hewlett-Packard Basic 5.0 for both 
the Hewlett-Packard 9300 series computer and 9852A series data 
acquisition/control unit. 
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1 I 

2 I 

3 I 
13 
Z 3 
25 

27 

28 
3® 

31 

32 

33 

34 

35 

36 

37 

38 
4® 
41 
4Z 

43 

44 

45 

46 

47 

48 
491 
5® 

51 

53 

54 

55 

56 
6 ® 
9® 
ie® 

1 1® 
12® 
14® 

141 

142 
15® 
151 
1521 

155 

156 

158 

159 
16® 
162 
163 

165 

166 

167 

168 
169 
17 ® 

171 

172 
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PROGRAM: SETUP PROGRAM FOR W 3852 A 
DATE: AUGUST 13.1991 

PROGRAMMER: LT DEAN SUGIYAMA 



REAL T( ®: 12) ,Sun( ®: 1 2 > ,Tenp(0: 12 ) ,E( 1:3) 

ON KEY 0,15 GOTO 27 

PRINTER IS 1 

PRINT 

PRINT 

PRINT USING "4X, ""SELECT OPTION 

PRINT USING "6X , " "0-MONI TOR SUMP 

PRINT USING "6X.""1 -MONITOR LIQUID 

PRINT USING "6X, HM Z-CHECK THERMOCOUPLES"" 

PRINT USIN6 "6X, ""3-CHECK MAIN HEATER 

PRINT USIN6 "6X. ""4-CHECK AUX HEATERS 

PRINT USIN6 "6X "5-EXIT PROGRAM 

PRINT USING "4X, *' "NOTE: KEY ® - ESCAPE 

BEEP 

INPUT Ido 

IF Ido>5 THEN Ido-5 
IF Ido-0 THEN 50 
IF Ido-1 THEN 155 
IF I do- 2 THEN 185 
IF Ido-3 THEN 209 
IF Ido-4 THEN Z®9 
IF Ido-5 THEN Z5Z 
PRINT 



PRINT 

PRINT "SUMP TEff*ERATURE DEG C " 

PRINT 

OUTPUT 709 i "RST" 

Tot-0 

FOR J-1 TO 5 

OUTPUT 709 1 "CONFMEAS TEMPT. 51 1 .USE 600" 

ENTER 709 » A 

Tot-Tot+A 

NEXT J 

TavQ-Tot/5 

PRINT USIN6 "4X.MD0.D0'* iTbvq 

BEEP 

PRINT 

WAIT 6® 

60T0 50 



PRINT 

PRINT "LIQUID TEMPERATURE DEG C" 

PRINT 

OUTPUT 709 i "RST" 

FOR 1-1 TO 2 
Sun(I )-0 
TCI)-0 
NEXT I 

FOR J-1 TO 5 

OUTPUT 709 i "CONFMEAS TEfFT ,508-509 .USE 600" 
FOR 1-1 TO 2 
ENTER 709 1 T < I ) 

Sua(I)-Sua(I)4T(I) 

NEXT I 
NEXT J 

FOR 1-1 TO 2 
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174 Tenp(I)-SuR(I>/5 

175 NEXT I 

176 Tbvq-( T o«p( 1 > + Tonp( 2 ) )/Z 

18© PRINT USING M 4X,HD0.DD” i Tavg 
181 BEEP 
18Z WAIT 6© 

183 60TO 155 

1841 

185 PRINT 

186 PRINT “CHANNEL TEff’ERATURE DEG C M 

187 OUTPUT 7©9r'RST" 

188 FOR 1-1 TO 1Z 

189 Sur( I >-© 

19© T( I )•© 

19Z NEXT I 

193 FOR J-1 TO 5 

194 OUTPUT 7©9 1 "CONFHEAS TEff>T ,5©©-51 1 ,USE 6©© M 

195 FOR 1-1 TO 1Z 

196 ENTER 7©9»T( I ) 

197 Sur( I )-Sur( I )+T( I ) 

198 NEXT I 

2©© NEXT J 

Z©1 FOR 1-1 TO 12 

2©Z Tbrp< I )-Sur( I )/5 

Z©3 PRINT TAB(3)iIiTAB(l5)iTe«p(I> 

2 ©4 NEXT I 

2 ©5 BEEP 

2 ©6 WAIT 5 

2 ©7 GOTO 185 

Z©8! 

2 ©9 PRINT 

21© OUTPUT 7©9 i**RST*' 

211 FOR 1-1 TO 3 

212 Sur(I)-© 

215 NEXT I 

216 FOR J-1 TO 5 

218 OUTPUT 7©9 i "CONFMEAS DCV.512-514.USE 6©© M 

219 FOR 1-1 TO 3 

221 ENTER 7©9iE(I> 

222 Sur( I )-Sur( I HE( I ) 

223 NEXT I 

225 NEXT J 

226 Anp-© 

228 FOR 1-1 TO 3 

229 IF 1-1 THEN Volt-SuR( I )/5 
23© IF 1-2 AND Ido-3 THEN 

231 PRINT -MAKE SURE VOLTAGE BOX IS SET TO MAIN HEATERS 

232 Arp-Sur<I>/5 

233 END IF 

234 IF 1-3 AND Ido-4 THEN 

235 PRINT "MAKE SURE V0LTA6E BOX IS SET TO AUX HEATERS’* 

236 Arp-Sur(I>/5 

237 END IF 

238 NEXT I 

2 39 Anp-ABS( Arp* 1.9182) 

24© Volt-ABS(Volt*25> 

241 Pouer- Volt* Arp 

242 Resletence-Volt/Anp 

Z43 IF Arp-© THEN Reeletence-© 

245 PRINT 
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246 BEEP 

Z47 PRINT M VOLTAGE* V ) CURRENT* A ) RESI STENCE* ohne ) POUER*U>" 
Z48 PRINT 

Z49 PRINT USING M 1 X , 5* NDOOOOO. DOO , 3X > M \ Voi t , Anp .Resistance .Power 

250 UAIT 5 

Z51 GOTO 209 

25Z BEEP 

253 PRINT 

Z54 PRINT 

255 END 
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APPENDIX G: DATA REDUCTION PROGRAM 



The following is a listing of the complete data-reduct ion 
program, DRP8, written in Hewlett-Packard Basic 5.0 for both 
the Hewlett-Packard 9300 series computer and 9852A series data 
acquisition/control unit. 
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10001 FILE NAME: DRP8 
1 004 ! DATE: October 19, 1984 

1000! REVISED: Auguat 14,1991 by LT Dean Sugiyana 

1012 ! 

1016 COM /I dp/ I dp 
1020 PRINTER IS 1 
10Z4 CALL Select 

10Z8 INPUT "UANT TO SELECT ANOTHER OPTION < 1 -Y ,0-N )?" , I ael 

1 03Z IF Iael-1 THEN GOTO 1024 

1036 BEEP 

1040 BEEP 

1044 PRINTER IS 1 

1048 PRINT "DATA COLLECTION/REPROCESSING COMPLETED" 

1052 END 

1056 SUB Main 

1060 COM /I dp/ I dp 

1064 COM /Cc/ C(7) , leal 

1068 COM /Uil/ OZ.Di.Do.L.Lu.Kcu 

1 07Z DIM Enf < 1Z ) ,T< 1Z ) ,Dla( 1 3) ,DZa( 1 3) ,Oia( 1 3) ,Doa( 1 3> ,La( 1 3) .Lua( 13) ,Kcua( 13) , 
Et1(Z0).EtZ(20),Tnf(4)1 15) ,Sur< 20) ,Tt< 1 Z ) ,Ta< 1 Z ) . Vt ( Z ) , Va< Z ) ,E1 iq< 12),E( 12) 

1076 DATA 0. 10086091 ,25727.94369.-767345.8295,78025595.81 
1080 DATA -9247486589. 6. 97688E-H1 .-2.6619ZE-H 3,3. 94078E-H4 
1084 READ C<*> 

1088! DATA "SROoth" , "High Flux" , "TherRoexel-E" ,"TherRoexel-HE" 

1092 DATA SROoth.High Flux , Turbo-B, High Flux Mod,Turbo-B Mod 
1096 READ Tn«<0 
1100 PRINTER IS 701 
1104 BEEP 

1108 IF Idp-4 THEN PRINTER IS 1 
1 1 1 Z IF Idp-4 THEN 60T0 Z03Z 
1116 OUTPUT 709 i "RST" 

11Z0 OUTPUT 709 i "SET TIMEDATE" iTIMEDATE 
1 1Z4 OUTPUT 709 i "TIMEDATE" 

1128 ENTER 709 sDt 
1 1 32 PRINT 

1136 PRINT " Month, date and tiRe : " »DATEt( Dt ) ,TIMES( Dt ) 

1140 PRINT 

1144 PRINT USING M 10X,""NOTE: PrograR naRe : DRP8" " " 

1148 BEEP 

1152 INPUT "ENTER DISK NUMBER", Dn 

1156 PRINT USING "16X.""Diek nuRber - "".ZZ"»Dn 

1160 BEEP 

1164 INPUT "ENTER INPUT MODE ( 0-3852A, 1 -FILE )" . Ir 
1168 BEEP 

1172 INPUT "SELECT HEATING MODE <0-ELEC» 1 -WATER)" ,Ih* 

1176 BEEP 

1180 INPUT "ENTER THERMOCOUPLE TYPE (0-NEU, 1 -OLD)" .leal 
1184 IF Ir-0 THEN 
1188 BEEP 

1192 INPUT "GIVE A NAME FOR THE RAU DATA FILE" ,DZ_f ilet 
1196 PRINT USING "16X,""Neu file nane: "" . 14A" »02lf ilet 
1200 Sizel -20 

1204 CREATE BDAT OZ_f ilet, Sizel 
1208 ASSIGN 0FileZ TO OZ_filet 
12121 

12161 DUMMY FILE UNTIL Nrun KNOUN 
1 2Z0 01_filet-"DUMMY" 

1 ZZ4 CREATE BDAT Ol.fileS, Sizel 
1 2Z8 ASSIGN CFilel TO Ol.filet 
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1232 OUTPUT ©FileUDt 
1 2 36 IF Ih*-0 THEN 
1240 BEEP 

1244 INPUT "ENTER NUMBER OF DEFECTIVE TCS < O-DEFAULT ) ", I dtc 
1248 IF Idtc-0 THEN 
1252 Ldtcl "0 
1256 LdtcZ-0 

1260 PRINT USING "16X,""No defective TCe exist 

1264 END IF 

1268 IF Idtc-1 THEN 

1272 BEEP 

1276 INPUT "ENTER DEFECTIVE TC LOCATION" .Ldtcl 

1280 PRINT USING "16X,""TC is defective at location "".D0"iLdtc1 

1284 Ldtc2-0 

1288 END IF 

1292 IF Idtc-2 THEN 

1296 BEEP 

1300 INPUT "ENTER DEFECTIVE TC LOCATIONS" .Ldtcl ,Ldtc2 

1304 PRINT USING ,, 16X,""TC are defective at locations "" .DD.4X.D0" iLdtcl . Ldtc2 

1308 END IF 

1312 IF Idtc>2 THEN 

1316 BEEP 

1320 PRINTER IS 1 

1324 BEEP 

1328 PRINT "INVALID ENTRY" 

1332 PRINTER IS 701 
1336 60T0 1240 

1340 END IF 
1344 END IF 

1348 OUTPUT tiFllel (Ldtcl ,Ldtc2 
13521 Ir- 1 option 
1356 ELSE 
1360 BEEP 

1364 INPUT "61 VE THE NAME OF THE EXISTING DATA FILE" ,02_f ileS 
1368 PRINT USING "l6X f ""01d file nane: "" . 14A" »D2_f ile« 

1372 ASSIGN «File2 TO D2_fileS 
1376 ENTER ®File2iNrun 
1380 ENTER «File2;Dt 

1384 PRINT USING M 16X I ""This data set taken on : "" . 1 1 A" \ DATE$( Dt ) 

1388 ENTER ®File2 iLdtcl .Ldtc2 
1392 IF Ldtcl >0 OR Ldtc2>0 THEN 

1396 PRINT USING " 1 6X , ""Thermocouples were defective at locations: "" ,2( 3D, 4X>" » 

Ldtcl ,Ldtc2 

1400 END IF 

1404 ENTER OFlleZtltt 

1408 END IF 

1409 Idtc-0 

1411 IF Ldtcl >0 THEN Idtc-Idtc+1 

1412 IF Ldtc2>0 THEN Idtc-Idtc-H 
14131 IF Ir- 0 AND Ihn-1 THEN 1595 
1416 BEEP 

1420 INPUT "WANT TO CREATE A PLOT FILE? < 0-N. 1-Y )" .Iplot 
1424 IF Iplot-1 THEN 
1428 BEEP 

1432 INPUT "61 VE NAME FOR PLOT FILE" ,P_f ile* 

1436 CREATE BOAT P_file$,4 
1440 ASSI6N BPlot TO P_filet 
1444 END IF 
1448 IF Ihn-1 THEN 
1452 BEEP 
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145G INPUT “WANT TO CREATE Uo FILE? ( O-N . 1 -Y ) " . I uf 
1460 IF Iuf-1 THEN 
1464 BEEP 

1468 INPUT "ENTER Uo FILE NAME" ,Uf ileS 
1472 CREATE BOAT Ufile«,4 
1476 ASSI6N CUfile TO UfileS 
1480 END IF 
1484 BEEP 

1488 INPUT "WANT TO CREATE Re FILE? < O-N , 1 -Y > " , I re 
1492 IF Ire-1 THEN 
1496 BEEP 

1500 INPUT “ENTER Re FILE NAME" .Ref ilet 

1504 CREATE BOAT RefileS. 10 

1508 ASSI6N CRefile TO RefileS 

1512 END IF 

1516 END IF 

1520 PRINTER IS 1 

1524 IF In-0 THEN 



1528 


BEEP 








1532 


PRINT 


USIN6 


"4X, ""Select tube number 


1536 


IF Ihn-0 THEN 




1540 


PRINT 


USING 


“6X , " “0 


Smooth 4 inch Ref 


1544 


PRINT 


USIN6 


"6X , "“ 1 


Smooth 4 inch Cu (Press/Slide)' 


1548 


PRINT 


USING 


"6X , " "2 


Soft Solder 4 inch Cu""" 


1552 


PRINT 


USING 


"6X," "3 


Soft Solder 4 inch HIGH FLUX”" 


1556 


PRINT 


USIN6 


“6X , “ "4 


Uleland Hard 8 inch 


1560 


PRINT 


USING 


"6X,"“5 


HI6H FLUX B inch 


1564 


PRINT 


USIN6 


“6X,"“6 


GEUA-K 40 Fins/ in 


1568 


PRINT 


USIN6 


“6X,““7 


6EUA-K 26 Fins/ in 


1572 


PRINT 


USING 


“6X,““B 


GEUA-T 19 F ins/ in“" " 


1576 


PRINT 


USIN6 


"6X," "9 


GEUA-T OR 6EUA-TY 26 Fins/ in 


1580 


PRINT 


USIN6 


“6X , " " 10 


THERMOEXCEL-E 


1584 


PRINT 


USING 


" 6X , “ ” 1 1 


THERMOEXCEL- HE" " “ 


1585 


PRIN* USING 


"6X , "" 12 


TURBO- B 


1586 


PRINT 


USIN6 


"6X, ““ 1 3 


GEUA-K 19 Fins/in 


1588 


ELSE 








1592 


PRINT 


USING 


“6X,““0 


Smooth tube 


1596 


PRINT 


USING 


"6X, " " 1 


High Flux 


1600 


PRINT 


USING 


“ 6X , " " 2 


Turbo-B““" 


1604 


PRINT 


USIN6 


"6X , " "3 


High Flux Mod“"“ 


1608 


PRINT 


USING 


“6X , " "4 


Turbo-B Mod““" 


1612 


END IF 






1616 


INPUT 


Itt 







1620 OUTPUT OF lie 1 1 1 tt 

1624 END IF 

1628 PRINTER IS 701 

1632 IF Itt<10 THEN PRINT USIN6 “16X.""Tube Number: ““.D"iltt 

1636 IF I tt>9 THEN PRINT USIN6 " 16X , ""Tube Number: “".DO'iItt 

1640 IF Ihn-1 THEN PRINT USIN6 “16X,““Tube Type: ““ . 15A" iTn*( I tt ) 

1644 BEEP 

1648 INPUT “ENTER OUTPUT VERSION (O-L0N6, 1 -SHORT, 2 -NONE)" ,Iov 
1652 BEEP 

1656 INPUT “SELECT <0-LIQ, 1-VAP.2-<LIQ+VAP>/2 ) M .Ilqv 
16601 

16641 DIMENSIONS FO TESTED TUBES 

16681 ELECTRIC HEATED MODE 

16721 Dl-Dieneter at thermocouple positions 

1676 DATA .01 1 1 125, .01 1 1 125, .01 1 1 125, .0129540, .012446, ,0129540. .0100965 
1680 DATA .0100965, .01 157, .01 157, .01 157, .01 157, .01 157, .0100965 
1684 READ D1a<*) 
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1688 D1-01a(Itt> 

16921 

16961 DZ-Diameter of toot section to the baoe of fins 

1700 DATA .015875,. 01 5875.. 01 5875,. 01 5824, .015875, .015824, .01270 

1704 DATA .01 27,. 01 38, .01 38,. 01 38,. 01 38,. 01 38,. 01 27 

1708 READ DZa( O 

17121 

1716! Di-Inside diameter of unenhanced ends 

1720 DATA .01 27,. 01 27, .0127, .01 32,. 01 27, .0132 ,.01 1 1 125, .01 1 1 125 

1721 DATA .01 18, .01 18, .01 18. .01 18, .01 18, .01 1 1 125 
1724 READ 0ia(O 

1728! 

1732! Do-Outside diameter of unenhanced ends 

1736 DATA .015875, .015875, .015875. .015824, .015875, .015824, .01270, .01270 

1737 DATA .01 331 , .01 33 1 , .01 331 . .01 331 ,. 01 58, .01 27 
1740 READ Doa( * ) 

1744! 

1748! L-LenQth of enhanced surface 

1752 DATA . 101 6, . 101 6 ,. 1016 ,. 101 6, . 2032 ,. 2032 ,. 2032 , .2032 ,. 2032 ,. 2032 ,. 2032 ,. 20 

32.. 2032.. 2032 
1756 REAOLatO 
17601 

17641 Lu-Length of unenhanced surface at the ends 

1 768 DATA . 02 54 , . 02 54 , . 02 54 , . 02 54 , . 0762 , . 0762 , . 0762 , . 0762 , . 0762 ,. 0762 0762 ,. 07 

62.. 0762.. 0762 
1772 READ Lua( • ) 

1776! 

17801 Kcu-Thermal Conductivity of tube 

1784 DATA 398,344,344.45,344,45,344,344,398,398,398,398,398,344 
1788 READ Kcua( * > 

1792 IF Ihm-1 THEN 
1796! 

1800! Data statements for water heating mode 
18041 

1808 DATA 0.015875,0.015875,0.0169,0.0138,0.0169,0,0,0,0,0 
1812 READ DZa( * ) 

1816 DATA 0.0127,0.0127,0.0145,0.0127,0.0145,0,0,0,0,0 
1820 READ Dia(t) 

1824 DATA 0.015875,0.015875,0.0169,0.015875,0.0169,0,0,0,0,0 
1828 READ Doa< « ) 

1832 DATA 0.3048,0.3048,0.3048,0.3048,0.3048,0,0,0,0,0 
1836 READ La(«) 

1840 DATA 0.0254,0.0254,0.0254,0.0254,0.0254,0,0,0,0,0 
1844 READ Lua< « ) 

1848 DATA 398,45,398,45,398,0,0,0,0,0 
1852 READ Kcua(«) 

1856 END IF 
1860 D2-D2a( Iti > 

1864 Di-Dia(Itt) 

1868 Do-Doa(Itt) 

1872 L-La< Itt > 

1876 Lu-Lua(Itt) 

1880 Kcu-Kcua(Itt) 

1884 Xn- . 8 
1888 Fr- . 3 

1892 IF Itt-0 THEN Cf-1.70E+9 
1896 IF ltt>0 THEN Cf-3.7037E+ 10 
1900 A-PI«< Do*2-Di *2 )/4 
1904 P-PI «Do 
1908 IF Ihm-1 THEN 
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1912 BEEP 

1916 INPUT “TUBE INITIATION MODE. < 1-HOT WATER, 2-STEAM . 3-COLD WATER)*' , 1 1 in 

1920 IF Itin-1 THEN PRINT USING *' 16X ,* ,, 'Tube Initiate: Hot Water 

1924 IF Itin-2 THEN PRINT USING '* 1 6X , “"Tube Initiate: Steam 

1928 IF 1 1 in-3 THEN PRINT USING *' 1 GX . “ “Tube Initiate: Cold Water 

1932 INPUT “TEMP/VEL MODE: ( 0-T- CONST, V-DECi 1 -T-DEC ,V- CONST * 2-T-INC . V-CONST )" , 
Itv 

1936 IF I tv-0 THEN PRINT USING “ 16X, ““Tenp/Vel Mode: Constant/Decreas ing 

1940 IF Itv-1 THEN PRINT USING “ 16X,"" Tenp/Vel Mode: Decreasing/Conetant 

1944 IF I tv-2 THEN PRINT USING “ 1 6X . ““Temp/ Vel Mode: Increasing/Constant 

1948 INPUT “WANT TO RUN WILSON PLOT? < 1-Y ,0-N >“ , Iuil 

1952 IF Ihn-1 AND Iuil-0 THEN 

1956 IF Itt-0 THEN Ci-,032 

1960 IF Itt-1 OR Itt-3 THEN C1-.059 

1964 IF Itt-2 OR Itt-4 THEN Ci-.062 

1968 BEEP 

1972 INPUT “ENTER Cl < DEF : WH-.032 ,HF-.059 ,TB- .062 >" ,Ci 

1976 PRINT USING “ 16X . "“Sieder-Tate 

1980 PRINT USIN6 “1GX,“" Constant - ““ ,Z. 4D“ i Ci 
1984 END IF 
1988 END IF 

1992 IF Ihn-1 AND In-1 AND Iuil-1 THEN 
1996 IF Itt-0 THEN Ci-.032 
2000 IF Itt-1 OR Itt-3 THEN Ci-.059 
2004 IF Itt-2 OR Itt-4 THEN Ci-.062 
2008 ASSI6N «File2 TO • 

2012 CALL Wilson< Cf ,Ci ) 

2016 ASSIGN ®File2 TO 02_fileS 

2020 ENTER ®Fil e2 i Nrun . Dol dt , Ldtcl , Ldtc2 . It t 

2024 END IF 

2028 Neub-0 

2032 IF I dp- 4 THEN Ihn-1 
2036 IF Ihn-1 THEN Nsub-8 

2040 Ntc-6 

2041 IF Ihn-0 THEN Ntc-12 

2044 J-1 

2048 Sx-0 
2052 Sy-0 
2056 Sxe-0 
2060 Sxy-0 
2064 Repeat: I 
2068 IF In-0 THEN 
2072 Dtld-2. 22 
2076 Ido-2 

2080 ON KEY 0,15 RECOVER 2064 
2084 PRINTER IS 1 

2088 PRINT USING “4X, ““SELECT OPTION 

2092 PRINT USING “6X,““0-TAKE DATA*' “" 

2096 IF Ihn-0 THEN PRINT USING “6X, ““1-SET HEAT FLUX“““ 

2100 IF Ihn-1 THEN PRINT USING “6X,““1-SET WATER FLOW RATE"** “ 

2104 PRINT USING **GX , ““2-SET Teat 

2108 PRINT USING "4X.““N0TE: KEY 0 - ESCAPE 

2112 BEEP 
2116 I1PUT Ido 
2120 IF Ido>2 THEN Ido-2 

2124 IF Ido-0 THEN 2724 

21281 

2132! LOOP TO SET HEAT FLUX OR FLOWMETER SETTING 
2136 IF Ido-1 THEN 

2140 IF Ihn-0 THEN 
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2148 BEEP 

2152 INPUT "ENTER DESIRED Qdp",Dqdp 

2155 PRINT USING "4X , " "DESI RED Qdp ACTUAL Qdp 

2 16$ Err- 1000 
2154 FOR 1-1 TO 2 

2165 Sun<I>-0 

2166 Vt(I>-0 
216B NEXT I 

2169 OUTPUT 709 i " RST" 

2176 FOR Ji-1 TO 5 

2180 OUTPUT 709 t "CONFMEAS DCV ,51 2-51 3 .USE 600" 

2181 FOR 1-1 TO 2 

2182 ENTER 709|Vt<I) 

2184 Su«(I>-Su«<IHVt(I) 

2185 NEXT I 

2188 NEXT Ji 

2189 FOR 1-1 TO 2 

2192 IF 1-1 THEN VoH-Su«( I >/5 
2196 IF 1-2 THEN A*p-Su*<I)/5 
2200 NEXT I 

2 201 Anp- ABS( Amp* 1.9182) 

2202 VoH-ABS(Voli*25) 

2204 Aqdp-VoH*Amp/<PI*D2*L) 

2208 IF ABS< Aqdp-Oqdp)>Err THEN 

2212 IF Aqdp>Dqdp THEN 

2216 BEEP 4000, .2 

2220 BEEP 4000,. 2 

2224 BEEP 4000,. 2 

2228 ELSE 

2232 BEEP 250, .2 

2236 BEEP 250,. 2 

2240 BEEP 250,. 2 

2244 END IF 

2248 PRINT USING “4X,HZ.3DE,ZX.MZ.3D£"iOqdp,Aqdp 

2252 WAIT 2 

2256 GOTO 2164 

2260 ELSE 

2264 BEEP 

2268 PRINT USING ‘^X.MZ.SDE^X.HZ.SDE^Oqdp.Aqdp 

2272 Err-500 

2276 WAIT 2 

2280 GOTO 2164 

2284 END IF 

2288 ELSE 

2292 BEEP 

2296 INPUT "ENTER FLOWMETER SETTING" .Fwo 

2300 GOTO 2088 

2304 END IF 

2308 END IF 

2312! 

2316! LOOP TO SET Teat 
2320 IF Ido-2 THEN 
2324 IF Ikdt-1 THEN 2344 
2328 BEEP 

2332 INPUT “ENTER DESIRED Toaf’.Dtld 

2336! PRINT USIN6 "4X, M " DTsai ATsat Rate Tv 

2340 Ikdt-1 

2344 Old 1-0 

2348 Old2-0 

2352 Nn-1 



Rate""" 
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2356 Nre-Nn MOD 15 
2360 Nn-Nn+1 



2364 
2 365 


IF 

IF 


Nra-1 

Ihn-0 


THEN 

THEN 


PRINT 


USING "4X , "" 


Teat 


Tldl 


Tld2 


Tv 


T sun 


P 

2 368 


IF 


Ihn-1 


THEN 


PRINT 


USING “4X , " " 


Teat 


Tldl 


Tld2 Tv 


Teunp 


Tinle 



t Tpile Tout 

2372 ENO IF 

2373 FOR 1-1 TO 6 

2374 Sun(I)-0 

2375 NEXT I 

2376 OUTPUT 709 i "RST" 

2377 FOR Ji-1 TO 20 

2379 IF Ihn-0 THEN OUTPUT 709 i “CONFMEAS TEMPT, 508-51 1 .USE 600" 

2380 IF Ihn-1 THEN OUTPUT 709 1 ’’CONFMEftS TEMPT, 300- 305, USE 600" 

2384 FOR 1-1 TO 6 

2388 IF Ihn-0 AND I >4 THEN 2401 
2392 ENTER 709»Eliq(I) 

2396 Sun( I >-Sun< I )+El iq< I ) 

2400 NEXT I 

2401 NEXT Ji 

2402 FOR 1-1 TO 6 

2403 IF Ihn-0 AND I>4 THEN 2496 
2416 El lq( I )-Sun< I )/20 

2424 IF 1-1 THEN Tldl-Eliq(I) 

2428 IF 1-2 THEN Tld2-Eliq(I) 

2432 IF 1-3 THEN Tv-Eliq(I) 

2436 IF 1-4 THEN Teunp-El iq( I ) 

2440 IF 1-5 THEN T inlot-El iq( I ) 

2444 IF 1-6 THEN Tout-El iq< I) 

2448 NEXT I 

2452 IF Ihn-1 THEN 

2456 Sunl-0 

2460 OUTPUT 709 i M RST** 

2468 FOR Kk — 1 TO 20 

2469 OUTPUT 709 T’CONFMEAS TEtFT , 320.USE 600 M 
2472 ENTER 709»E1 

2476 Sun 1- Sun 1 +E 1 

2480 NEXT Kk 

2484 Enf ( 7 )-AB$( Sunl /20) 

2488 Tpile-Enf < 7)/3.96E-4 
2492 ENO IF 

2496 Atld-( Tldl +Tld2 )* . 5 

2500 IF ABS( Atld-Dt Id) >. 2 THEN 

2504 IF Atld>Dtld THEN 

2508 BEEP 4000,. 2 

2512 BEEP 4000,-2 

2516 BEEP 4000, .2 

2520 ELSE 

2524 BEEP 250,. 2 

2528 BEEP 250, .2 

2532 BEEP 250,. 2 

2536 END IF 

2540 Errl-Atld-Oldl 

2544 Oldl-At Id 

2548 Err2-Tv-01d2 

2552 01d2-Tv 

2553 IF Tldl > 100. THEN 2572 

2556 IF Ihn-0 THEN PRINT USING M 4X,5< MD00.DD,2X)” lOtld.Tldl ,Tld2 .Tv.Tounp 
2560 IF Ihn-1 ANO Idp-0 THEN PRINT USING "4X,7< MDD.D0,2X)* , iDtld,Tld1 ,Tld2 .Tv.Tb 
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uwp ,Tinl et , Tpi 1 e 

2564 IF Ih*-1 AND Idp-4 THEN PRINT USING '* 4X , 5< MDO .DD,2X) , 3< M30. DO , 2 X > i Dt 1 d , T 1 

dl ,Tld2 , Tv , Taunp.T inlet , Tpile.Tout 

2568 WAIT 2 

2572 GOTO 2356 

2576 ELSE 

258® IF ABS( At ld-Ot 1 d ) > . 1 THEN 

2584 IF At 1 d>Ot Id THEN 

2588 BEEP 3®®®.. 2 

2592 BEEP 3®®®,. 2 

2596 ELSE 

26®® BEEP 8®®,. 2 

26®4 BEEP 8®®,. 2 

26®8 END IF 

2612 Errl-Aild-Oldl 

2616 Oldl-At Id 

262® Err2-Tv-01d2 

2624 01d2-Tv 

2628 IF Ih*-® THEN PRINT USING **4X,5< MDOO. DO, 2X ) M i Otld.Tldl , Tld2 ,Tv , Tau«p 

2632 IF Ihn-1 THEN PRINT USING *’4X.5< HD0.D0.2X) , 3( M3D.DO, 1 X)" iDtld.Tldl ,Tld2 ,Tv 

,Tsu«p t Tlnlet .Tplle .Tout 

2636 WAIT 2 

264® GOTO 2356 

2644 ELSE 

2648 BEEP 

2652 Err 1 -At Id-Old 1 

2656 Oldl-Atld 

266® Err2-Tv-01d2 

2664 01d2-Tv 

2668 IF Ih«-® THEN PRINT USING "4X,5< MDOO. DO, 2 X >*' iDtld.Tldl ,Tld2 ,Tv ,Tau«p 

2672 IF Ih«-1 THEN PRINT USING "4X , 8< ADO. DD, 2X ) M iDild.Tldl ,Tld2,Tv, Taunp , Tinl et 

.Tplle .Tout 



2676 


WAIT 2 








268® 


60T0 2356 








2684 


END IF 








2688 


END IF 








2692 


END IF 








2696! 


ERROR TRAP FOR Ido OUT OF 


BOUNDS 






27®® 


IF Ido>2 THEN 








27®4 


BEEP 








27®8 


60T0 2®88 








2712 


END IF 








2716! 










272®! 


TAKE DATA IF I*-® LOOP 








2724 


IF Ikol-1 THEN 2737 








272 8 


BEEP 








2732 


INPUT "ENTER BULK OIL X*\ 


Bop 






2736 


Ikol-1 








2737 


FOR 1-1 TO Ntc 








2738 


E(I >-® 








2739 


Su«d>-® 








274® 


NEXT I 








2741 


OUTPUT 7®9 1 **RST*' 








2742 


FOR Jl-1 TO 2® 








2744 


IF Ih«-® THEN OUTPUT 7®9i 


"CONFMEAS 


TEMPT, 5®®-5l 1 .USE 


6®® M 


2745 


IF Ihw-1 THEN OUTPUT 7®9i 


"CONFMEAS 


TEMPT. 3®®- 3®5. USE 


6®®*' 


2748 


IF Ih«-® THEN Ntc- 12 








2752 


FOR 1-1 TO Ntc 








2756 


ENTER 7®9 i E< I ) 








276® 


Su«( I )-5uw< I )+E< I ) 
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Z776 IF I-< 9-Neub) THEN Et 1 ( Ji- 1 )-E< I > 

Z777 IF I-d0-Neub> THEN EtZ< Ji- 1 >-E< I ) 

Z778 NEXT I 
Z780 NEXT Ji 
Z7B4 Kdl-0 
Z7B5 FOR 1-1 TO Ntc 

Z788 IF I-( 9-Ncub) OR I-dO-Naub) THEN 
Z79Z Eave-Sum( I )/Z0 
Z797 Suw<I)-0 
Z798 END IF 

Z8O0 IF I-( 9-Neub ) THEN 

Z801 FOR Jk-0 TO 19 

Z804 IF ABS( Et 1 < Jk )-Eave X • 1 THEN 

2808 Sun< I )-Suw< I MEt 1 ( Jk ) 

Z81Z ELSE 
Z81B Kdl-Kdl+1 
Z817 END IF 
Z 8 1 9 NEXT Jk 
Z8ZO END IF 

Z8ZZ IF I-< lO-Nsub) THEN 

Z8Z3 FOR Jk-0 TO 19 

28ZS IF ABS< EtZ( Jk )-EaveK. 1 THEN 

Z8ZB Sun( I >-Suw< I >*EtZ< Jk ) 

Z8Z8 ELSE 
Z8Z9 Kdl-KdHI 
Z830 END IF 
Z83Z NEXT Jk 
Z833 END IF 

Z83B IF I-( 9-Neub) OR I-d0-Neub> THEN PRINT USING •*4X." , ‘Kdl - MM .DO"iKdl 
Z837 IF Kdl >10 THEN 
ZB38 BEEP 
Z840 BEEP 

Z844 PRINT USING "4X, M “Too nuch ecatterinQ in data - repeat data eet 

Z848 GOTO Z084 
Z85Z END IF 

2860 Ewf < I >-Suw< I )/( 20-Kdl > 

Z8B4 NEXT I 
Z8B8 IF Ih«-1 THEN 
Z87Z Sun< 1 )-0 

Z87B OUTPUT 709 » “CONFMEAS TEMPT , 3Z0.USE 600“ 

Z884 FOR Kk-1 TO ZO 
Z888 ENTER 709|E(Kk) 

Z89Z Suw< 1 >-Su*< 1 ME(Kk) 

Z89B NEXT Kk 

Z900 Enf(7)-ABS(Sun( 1 ))/Z0 

Z904 END IF 

Z908 IF Ihw-0 THEN 

Z909 FOR 1-1 TO Z 

Z910 Su«( I )-0 

Z911 NEXT I 

Z91Z OUTPUT 709i“RST“ 

Z 9 1 3 FOR Ji-1 TO 5 

Z915 OUTPUT 709 t “CONFMEAS DCV .51Z-S1 3 .USE BOO" 

Z91B FOR 1-1 TO Z 
Z9Z0 ENTER 709 \ E<I) 

Z93B Stm<I)-Sun<X)+E<I) 

Z937 NEXT I 

Z940 NEXT Ji 

Z94 1 FOR 1-1 TO Z 

Z944 IF 1-1 THEN Vr-Sun( I )/5 
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2948 IF 1-2 THEN Ir-Suw< I )• 1 .9 1 82/5 
2952 NEXT I 
2956 END IF 
2960 ELSE 

2964 IF Ihn-0 THEN ENTER tFile2iBop.Told.Emf ( • ) .Vr ,Ir 
2968 IF Ihm-1 THEN ENTER tFile2iBop.Told.Enf < * > ,Fna 
2972 END IF 
29761 

29801 CONVERT enf’S TO TEMP .VOLT .CURRENT 

2984 Tuo-0 

2988 FOR 1-1 TO Ntc 

2992 IF ldtc>0 THEN 

2996 IF I-Ldtcl OR I-Ldtc2 THEN 

3000 T( I )-0 

3004 60T0 3044 

3008 END IF 

3012 END IF 

3016 IF I tt <4 AND Ihn-0 THEN 

3020 IF I>4 AND I <9 THEN 

3024 T( I >-0 

3028 60T0 3044 

3032 END IF 

3036 END IF 

3040 T( I )-Enf ( I ) 

3044 NEXT I 

3048 IF Itt<4 THEN 

3052 FOR 1-1 TO 4 

3056 IF I-Ldtcl OR I-Ldtc2 THEN 

3060 Tua-Tua 

3064 ELSE 

3068 Tua-Tua+T( I ) 

3072 END IF 
3076 NEXT I 
3080 Tua-Tua/(4-Idtc ) 

3084 ELSE 

3088 IF I ha- 1 THEN 3128 

3092 FOR 1-1 TO 8 

3096 IF I-Ldtcl OR I-Ldtc2 THEN 

3100 Tua-Tua 

3104 ELSE 

3108 Tua-Tua+Td) 

3112 END IF 
3116 NEXT I 
3120 Tu-Tua/( 8-Idtc) 

3124 END IF 
3128 Tld-T< 9-Ntub ) 

3132 Tld2-T( 10-Neub) 

3136 Tlda-( Tld+Tld2 )* . 5 
3140 Tv-T( 1 1 -Neub) 

3158 Ttunp-T( 12-Neub) 

3160 IF Ihn-0 THEN 3176 
3168 Tinlet-T< 1 3-Neub) 

3172 Tout-T( 14-Neub) 

3176 IF Iha-0 THEN 
3180 Aap-ABS(Ir) 

3184 Volt-ABS( Vr )«25 
3188 Q-Volt*Aap 
3192 END IF 

3196 IF Itt-0 AND Iha-0 THEN 
3200 Kcu-FNKcu(Tu) 
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32 ©4 ELSE 

3208 Kcu-Kcua< I tt ) 

3212 END IF 
3216' 

3220! FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
3224 IF Ihn-0 THEN Tu-Tu-Q*LOG< D2/D1 )/( 2*PI*Kcu*L ) 

3228 IF Ilqv-0 THEN Teat-Tlda 

3232 IF Ilqv-1 THEN Tsat-< T1 da4Tv )• . 5 

3236 IF Ilqv-2 THEN Teat-Tv 

3244 IF Ih*-1 THEN 

3248 Tavg-Tinlet 

3252 Grad-37.98534. 104388* Tavg 

3256 Tdrop-ABS(ERf<7))»1.E4G/( 10*Grad> 

3260 Tavgc-Tinlet-Tdrop* .5 

3264 IF AB5( Tavg-Tavgc ) > . 01 THEN 

3268 Tavg-< Tavg4Tavgc )• . 5 

3272 GOTO 3252 

3276 END IF 

3280! 

32841 COMPUTE WATER PROPERTIES 
3288 IF Ihn-1 THEN 
3292 Ku-FNKu(Tavg) 

3296 Muua-FNMuu< Tavg) 

3300 Cpu-FNCpu< T bvq ) 

3304 Pru-FNPru(Tavg) 

3308 Rhou-FNRhou< Tavg ) 

3312 Twi-Tavg 
33161 

3320! Conpute MOOT 

3324 Mdot-3.9657E-34Fne*<3.61955E-3-Fne*<8.82006E-6-Fne*< 1 .23688E-7-F«e*4. 31897 

E- 10) ) ) 

33281 Mdot-Mdot*< 1 .0365-Tinl et*< 1 .96644E-3-Tinlet*5. 252E-6) >/1 .0037 
3332 Kdt-0 

3336 Q-Mdot*Cpu*Tdrop 

3340 L«td-Tdrop/LOG< < T Inlet -Tsai )/( Tinlet-Tdrop-Teat ) ) 

3344 Uo-Q/<PI*Do*L*Lntd) 

3348 Ru-Do*LOG< Do/Di )/< 2.*Kcu) 

3352 Tu*Teat4-Fr*Lntd 

3356 Vu-Mdot/< Rhou*PI *Di* 2/4 ) 

3360 Rew-Rhow*Vw*Dl/Muua 

3364 Hi-Ci*Ku/Di*Rew* . 8*Prw*< 1/3. )• < Muwa/FNMuw( Twi ) >* . 14 
3368 Tuic-Tavg-Q/< PI*Do*L*Hi ) 

3372 IF ABS< Tul-Twic)>.01 THEN 

3376 Twi-< Twl4Twic )• . 5 

3380 GOTO 3364 

3384 END IF 

3388 Twt-( Twl4Twlc )• . 5 

3392 Ho- 1 / < 1/Uo-Do/(Di*Hi)-Rw) 

3396 END IF 

3397 END IF 

3399 IF Ihw-1 THEN 

3400 Thetab-Q/< Ho*PI*Do*L ) 

3404 Tu-Tsat4Thet ab 

3405 ELSE 

3408 Thetab-Tu-Teat 
3411 END IF 
3413 IF Thetab<0 THEN 
3416 BEEP 

3420 INPUT "TUALL<TSAT (0-CONTINUE, 1-END) M ,Iev 
3424 IF Iev-0 THEN GOTO 2068 



* 
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3420 IF lev- 1 THEN 3832 

3432 END IF 

34401 

3444! COMPUTE VARIOUS PROPERTIES 
3452 Tfiln-(Tu*Teot)*.5 
3456 Rho-FNRho< Tf i I m ) 

3460 Mu-FNMutTfilw) 

3464 K-FNK(Tfiln) 

3468 Cp-FNCp(Tfilm) 

3472 Bet a-FNBe t a< T f i 1 n ) 

3476 Hf Q-FNHf q( T cot ) 

3480 Ni-Mu/Rho 
3404 Alpha-K/( Rho«Cp) 

3488 Pr-Ni/Alpho 
3492 Peat-FNPeattTeat > 

3496! 

35001 COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
35041 FOR UNENHANCED END(S) 

3508 Hbar-190 

3512 Fe-< Hbar«P/< K cu« A ) ) A . 5*Lu 
3516 Tonh-FNTonh(Fe) 

3520 Theta-Thetab«Tanh/Fe 

3524 Xx-(9.81»BetB«Thetab*Do*3«Tanh/(Fe»Nl»Alpha) )* . 166667 
3528 Yy-< 1+( .559/Pr)*( 9/ 16) r< 8/27) 

3532 Hbarc-K/Do*( .6+. 387«Xx/Yy r 2 

3536 IF ABS( ( Hbai — Hbar c ) / Hbar c ) > . 00 1 THEN 

3540 Hbar -( Hbar ♦ Hbar c > • . 5 

3544 GOTO 3512 

3548 END IF 

35521 

3556! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) 

3560 Ql-(Hbar*P*Kcu«A)* . 5«Thetab«Tanh 
3564 Qc-Q-2«Q1 

3568 Ab-PI»D2«L 
3572! 

3576! COMPUTE ACTUAL HEAT FLUX AND B0ILIN6 COEFFICIENT 

3580 Qdp-Qc/As 

3584 Htube-Qdp/Thetab 

3588 Cef •( Cp«Thetab/Hf q)/( Q dp/( Mu«Hfg)«( . 014/(9. 8MRhor.5K( 1/3. >«Pr*1 .7) 

3592! 

3596! RECORO TIME OF DATA TAKING 

3600 IF In-0 THEN 

3604 OUTPUT 709 1 " TI MEDATE“ 

3608 ENTER 709 t Told 

3612 END IF 

3616! 

3620! OUTPUT DATA TO PRINTER 
3624 PRINTER IS 701 
3628 IF Iov-0 THEN 
36 32 PRINT 

3636 PRINT USING M 10X. M ”Data Set Nunber - "" ,DDO,2X t M "Bulk Oil X - MM f DO.D t 5X,1 
4A” i J , Bop .TIME t( Told ) 

3640 IF Ihn-0 THEN 

3644 PRINT USING M 10X,""TC No: 1 2 3 4 5 6 7 

8“"" 

3640 PRINT USING " 10X, ""Tenp : — ,8< 1 X.MOO.OO)" i T( 1 ) , T( 2 ) ,T< 3) ,T< 4 ) , T< 5 ) ,T( 6) . T< 
7) ,T( 8) 

3652 PRINT USING , *10X. ,,M Tua Tliqd Tliqd2 Tvapr Peat Teunp 

3656 PRINT USIN6 M 10X,2( MDO.DO, IX) ,1X, MDO.DO. 1X t 2( IX, MOD. DO) ^X.MDO.D^Tu.Tld, T 
ld2 , Tv .Peat , Teunp 
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3660 PRINT USING "10X, MM Thetab Htube Odp 

3664 PRINT USING M 1 0X , HDD. 3D . 1 X , HZ. 3DE . I X , MZ . 3DE" « Thetab , Htube , Odp 
3668 ELSE 

367Z PRINT USING M 10X l " M Fne Vu Teat Tinl Tdrop Thetab q Uo 

Ho 

3676 PRINT USING " 1 0X, 4( ZD. DO, 1 X ) .2 . 3D , 1 X , DO. DO , 1 X , 3( H2. 3DE , 1 X) " ; Fne . Vu , Teat . T i 

nlet , Tdrop, Thet ab, Odp, Uo, Ho 

3680 END IF 

3684 END IF 

3688 IF Iov-1 THEN 

369Z IF J-1 THEN 

3696 PRINT 

3700 IF Ihn-0 THEN 



3704 

3708 


PRINT 

ELSE 


USING 


M 10X, MM 


RUN No 


Oil! 


Teat 


Htube 


Qdp 


Thetab' 


371 Z 


PRINT 


USING 


"10X, ,,M 


FMS 


OIL! 


TSAT 


HTUBE 


QDP 


THETAB' 



3716 END IF 
37Z0 END IF 
37Z4 IF Ihn-0 THEN 

37Z8 PRINT USING • , 1ZX,3D,4X,D0,ZX,HDD.D0,3< 1X.HZ.3DE)"i J.Bop.Teat .Htube ,Qdp, The 
tab 

373Z ELSE 

3736 PRINT USING M 1 ZX , 3D.4X ,DO,ZX ,MDO. DO, 3< 1 X ,MZ. 3DE >" iFne , Bop, Teat .Htube ,Qdp , T 
het ab 

3740 END IF 
3744 END IF 
3748 IF In-0 THEN 
375Z BEEP 

3756 INPUT “OK TO STORE THIS DATA SET ( 1 -Y t 0-N>?" ,0k 
3760 END IF 

3764 IF Ok- 1 OR In-1 THEN J-J+1 
3768 IF Ok- ! AND In-0 THEN 

377Z IF Ihn-0 THEN OUTPUT CFilel i Bop , Told , Enf ( • ) ,Vr , Ir 
3776 IF Ihn-1 THEN OUTPUT CFilel i Bop , Told , Enf ( • ) ,Fne 
3780 END IF 

3784 IF Iuf-1 THEN OUTPUT CUfileiVu.Uo 
3788 IF Ire-1 THEN OUTPUT CRef ileiFno ,Reu 

379 Z IF (In-1 OR Ok - 1 ) AND Iplot-1 THEN OUTPUT OPlot i Qdp .Thetab 
3796 IF In-0 THEN 
3800 BEEP 

3804 INPUT “WILL THERE BE ANOTHER RUN ( 1 -Y ,0-N )?" ,Go_on 
3808 Nrun-J 

381 Z IF Go_on-0 THEN 383Z 
3816 IF Go_on<>0 THEN Repeat 
36Z0 ELSE 

38Z4 IF J<Nrun+ 1 THEN Repeat 
38Z8 END IF 
383Z IF In-0 THEN 
3836 BEEP 

3840 PRINT USING - 1 0X , “ M NOTE : "" .22, MM data rune Mere stored in file M,, ,10A ,, iJ- 
1 ,DZ_f ile* 

3844 ASSIGN ®File1 TO • 

3848 OUTPUT CFileZ iNrun- 1 
385Z ASSIGN CFilel TO Dl_filef 
3856 ENTER CF il e 1 i Dt , Ldt cl . LdtcZ , lit 
3860 OUTPUT CFileZiDi.Ldtcl .LdtcZ.Itt 
3864 FOR I- I TO Nrun-1 
3868 IF Ihn-0 THEN 

387Z ENTER CFilel iBop.Told, Enf (•> .Vr.Ir 
3876 OUTPUT CFileZ iBop . Told, Enf (*> ,Vr . Ir 
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3880 ELSE 

3884 ENTER CFilel i Bop . Told ,E*f< •> .F rb 
3888 OUTPUT CF i 1 eZ i Bop . Tol d . E«H • ) ,Fms 
389 Z END IF 
3896 NEXT I 

3900 ASSIGN ©File 1 TO • 

3904 PUR6E "DUMMY" 

3908 END IF 
39 1 Z BEEP 
3916 PRINT 

39Z0 IF Iplot«1 THEN PRINT USING " 10X.""NOTE: "" ,ZZ."" X-Y pairs were stored in 
plot date file "" . 10A" i J- 1 ,P_f il e$ 

39Z4 ASSIGN CFileZ TO * 

39Z8 ASSIGN ©Plot TO • 

3932 IF Iuf-1 THEN ASSIGN ©Ufile TO * 

3936 IF Ire-1 THEN ASSI6N ©Refile TO • 

3940 CALL Stats 
3944 BEEP 

3948 INPUT "LIKE TO PLOT DATA < 1-Y.0-N)?*\Ok 
39 5Z IF Ok-1 THEN CALL Plot 
39S6 SUBEND 
39601 

3964! CURVE FITS OF PROPERTY FUNCTIONS 
3968 DEF FNKcu< T ) 

397Z I OFHC COPPER Z50 TO 300 K 

3976 Tk-T4Z73.15 I C TO K 

3980 K-434- • 1 1 Z • Tk 

3904 RETURN K 

3988 FNEND 

39 9 Z DEF FNflu(T) 

39961 170 TO 360 K CURVE FIT OF VISCOUSITY 

4000 Tk-T4Z73.15 1C TO K 

4004 Mu-EXPt -4.4636+( 101 1 ,47/Tk ) )• 1 .0E-3 

4008 RETURN Mu 

401 Z FNEND 

4016 DEF FNCp(T) 

40Z0I 180 TO 400 K CURVE FIT OF Cp 
40Z4 Tk-T4Z73.15 1C TO K 

40Z8 Cp-. 40 1884 1 . 65007E-3»Tk4 1 ,51494E-6*Tk *2-6. 67853E- 10«Tk A 3 

403Z Cp-Cp* 1000 

4036 RETURN Cp 

4040 FNEND 

4044 DEF FNRho< T ) 

4048 Tk-T4273.15 1C TO K 

405Z X-1-< 1.8«Tk/753.9S) IK TOR 

4056 Ro-36.3Z461 . 146414*X*< 1/3)4 16.41 801 5*X4 17. 476838*X* .541 . 1 19BZB*X*Z 

4060 Ro-Ro/ .06Z4Z8 

4064 RETURN Ro 

4068 FNEND 

407 Z OEF FNPr< T ) 

4076 Pr-FNCp< T ) *FNhu< T ) /FNK( T ) 

4080 RETURN Pr 
4004 FNEND 
4088 DEF FNK(T) 

409 Z 1 T<360 K WITH T IN C 
4096 K • . 07 1 - . 000Z 6 1 • T 
4100 RETURN K 
4104 FNEND 
4108 DEF FNTanh(X) 

4112 P-EXP< X ) 
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411 6 Q-1/P 

4120 Tanh-(P-Q)/<P+Q) 

4124 RETURN Tanh 
4128 FNEND 
4138! 

4184 DEF FNBata(T) 

4188 Rop-FNRho(T+. 1 ) 

4192 Row-FNRho( T- . 1 ) 

41 96 B«ta--2/( Rop+Ro* ) •( Rop-Rom )/ . 2 

4200 RETURN Bata 

4204 FNEND 

4208 DEF FNHf Q< T) 

4212 Hf Q-1 .3741 344E+5-T*< 3. 3094361 E+2+T*1 .2165143) 

4216 RETURN Hfg 

4220 FNEND 

4224 DEF FNPaat(Tc) 

42281 0 TO 80 deo F CURVE FIT OF Paat 
4232 Tf - 1 . 8*Tc+32 

4236 Pa - 5. 9 4552 5+ Tf • ( . 1 5352082+Tf • ( 1 .4B40963E-3*Tf *9 .6150671E-6) ) 

4240 Pg-Pa-14.7 

4244 IF Pq> 0 THEN I +-PSIG,--ln Hg 

4248 Paat-Pg 
4252 ELSE 

4256 Paat-Pg* 29. 92/14. 7 
4260 END IF 
4264 RETURN Past 
4268 FNEND 

4272 DEF FNHawooth< X .Bop , I aat ) 

4276 DIM A< 5 ) ,B( 5 ) ,C< 5 ) , D( 5 ) 

4280 DATA . 20526 B . 25322 ,. 3 1 9048 ,. 55322 , .79909 B 1 . 00258 
4284 DATA .74515, .72992 , .73189 . .71225. .68472. .64197 
4288 DATA .41092 ,. 17726, . 25142 , .54806, .819 16, 1 .0845 
4292 DATA . 7 1 403 , , 729 1 3 , . 72565 , . 69669 1 , . 665867 ,.61889 
4296 READ A< • ) ,B< • ) , C< • ) , D< • ) 

4300 IF Bop<6 THEN I-Bop 
4304 IF Bop-6 THEN 1-4 
4308 IF Bop-10 THEN 1-5 
4312 IF I aat- 1 THEN 
4316 He-EXP< A< I ) + B( I )*LOG( X) ) 

4320 ELSE 

4324 Ha-EXP(C(I ) + D( I )*L06(X) ) 

4328 END IF 
4332 RETURN Ha 
4336 FNEND 
4340 DEF FNPoly(X) 

4344 CON /Cply/ A< 1 0 , 1 0) ,C( 1 0) ,B( 5 ) , Nop , Iprnt , Opo , 1 1 oq , I f n , I j oin ,Nj oin 
4348 XI -X 
4352 Poly-B<0> 

4356 FOR 1-1 TO Nop 

4360 IF Ilog-1 THEN Xl-LOG(X) 

4364 Poly-Pol y+B( I )*X1 * I 
4368 NEXT I 

4372 IF Ilog-1 THEN Pol y-EXP( Pol y ) 

4376 RETURN Poly 
4380 FNEND 
4384 SUB Poly 

4388 DIM R( 10) ,S( 10) ,Sy( 12) ,Sx< 12) ,Xx( 100) ,Yy( 100) 

4392 COM /Cply/ A< 10, 10) ,C( 10) ,B( 5) ,N, Iprnt , Opo, I log, If n, Ijoin.Njoin 
4396 COM /Xxyy/ Xp( 25 ) , Yp( 25 ) 

4400 FOR 1-0 TO 4 
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4404 

4406 

441 Z 

4416 

4420 

44Z4 

44Z8 

4432 

4436 

4440 

4444 

4448 

445Z 

4456 

4460 

4464 

4468 

447Z 

4476 

4480 

4484 

4488 

44SZ 

4496 

4500 

4504 

4508 

4512 

4516 

45Z0 

45Z4 

45Z8 

453Z 

4536 

4540 

4544 

4548 

4552 

4556 

4560 

4564 

4568 

457Z 

4576 

4580 

4584 

4588 

459Z 

4596 

4600 

4604 

4608 

4612 

4616 

4620 

4624 

4628 

4632 

4636 

4640 



B(I)-0 
NEXT I 
BEEP 

INPUT “SELECT ( 0-FILE . 1 -KEYBOARD, 2 -PROGRAM )",In 

In-In4 1 

BEEP 

INPUT “ENTER NUMBER OF X-Y PAIRS" ,Np 

IF In-1 THEN 

BEEP 

INPUT “ENTER DATA FILE NAME" ,D_f ilet 
BEEP 

I**»UT “LIKE TO EXCLUDE DATA PAIRS ( 1 -Y ,0-N )?" , led 

IF Ied-1 THEN 

BEEP 

IW=UT “ENTER NUMBER OF PAIRS TO BE EXCLUDED" , Ipex 
END IF 

ASSIGN ®File TO D_file® 

ELSE 

BEEP 

ir#»UT “WANT TO CREATE A DATA FILE ( 1-Y,0-N>?“ .Yea 

IF Yea-1 THEN 

BEEP 

INPUT “6IVE A NAME FOR DATA FILE” ,D_f lie® 

CREATE BOAT 0_file$,5 
ASSIGN CFile TO Outlet 
ENO IF 
ENO IF 
BEEP 

INPUT “ENTER THE ORDER OF POLYNOMIAL" ,N 

FOR 1-0 TO N*Z 

Sy( I)-0 

Sx< I )-0 

NEXT I 

IF Ied-1 AND In-1 THEN 
FOR 1-1 TO Ipex 
ENTER CFileiX.Y 
NEXT I 
ENO IF 

FOR 1-1 TO Np 
IF In-1 THEN 

IF Opo-Z THEN ENTER ©FileiX.Y 
IF Opo<2 THEN ENTER CFileiY.X 
IF Opo-t THEN Y-Y/X 
IF IIoq-1 THEN 
IF Opo-Z THEN Xt-X/Y 
X-LOG(X) 

IF Opo-Z THEN Y-LOG(Xi) 

IF Opo<2 THEN Y-LOG(Y) 

END IF 
END IF 

IF In-2 THEN 
BEEP 

INPUT “ENTER NEXT X-Y PAIR".X,Y 
IF Yea-1 THEN OUTPUT «FileiX,Y 
END IF 

IF In<3 THEN 
Xx ( I )-X 
Yy( I )-Y 
ELSE 

X-Xp(I-1 ) 
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4644 Y-Ypd-1) ' 

4648 END IF 
4652 R<®)-Y 

4656 Sy<®>-Sy<0HY 
466® S< I )-X 
4664 Sx( 1 )-Sx< 1 HX 
4668 FOR J-1 TO N 
4672 R< J)-R< J-1 >«X 
4676 Sy( J)-Sy( J)+R< J) 

468® NEXT J 
4684 FOR J-2 TO N*2 
4688 S< J>-S< J-1 >«X 
4692 Sx< J)-Sx( J)4S< J) 

4696 NEXT J 
47®® NEXT I 

47®4 IF Yes-1 AND In-2 THEN 
47®8 BEEP 

4712 PRINT USING "12X.DO,"" X-Y pairs were stored in file **" . 1 ®A" ? Np . D_f i let 

4716 END IF 

472® Sx(®)-Np 

4724 FOR I-® TO N 

4728 C< I )-Sy< I ) 

4732 FOR J-® TO N 
4736 ft(I,J)-Sx(l4j) 

474® NEXT J 
4744 NEXT I 
4748 FOR I-® TO N- 1 
4752 CALL Divided) 

4756 CALL Subtract HI ) 

476® NEXT I 

4764 B< N)-C( N )/A< N.N ) 

4768 FOR I-® TO N- 1 
4772 B( N- 1-1 )-C( N- 1 - I ) 

4776 FOR J-® TO I 

478® B< N- 1 - I )-B( N- 1 -I )-A(N-1-I ,N-J)*B(N-J) 

4784 NEXT J 

4788 B< N- 1 -I )-B( N- 1 -I )/ A( N- 1 -I , N- 1 - 1 ) 

4792 NEXT I 

4796 {PRINTER IS 7®1 

48®® IPRINT B( • ) 

48®4 {PRINTER IS 7®5 
48®8 IF Iprnt-® THEN 

4812 PRINT USING ** 1 2X /“'EXPONENT COEFFICIENT*** 4 " 

4816 FOR I-® TO N 

482® PRINT USING '* 1 5X ,D0.5X ,HD.7DE" 1 1 ,B( I ) 

4824 NEXT I 
482 8 PRINT ** ** 

4832 PRINT USING ** 1 2X . ""DATA POINT X Y Y< CALCULATED ) DISCR 

EPANCY 

4836 FOR 1-1 TO Np 
484® Yc-B< ® ) 

4844 FOR J-1 TO N 
4848 Yc-Yc+Bt J)«Xxdr J 
4852 NEXT J 
4856 D-Yy( I )- Yc 

486® PRINT USING ** 1 5X , 3D ,4X . 4< HD. 5DE . 1 X )** » I , Xx ( I ) , Yy( I ) , Yc ,D 

4864 NEXT I 

4868 END IF 

4872 ASSIGN ©File TO * 

4876 SUBEND 
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4880 SUB Divide(M) 

4884 COM /Cply/ A( 10. 1 0 > ,C< 10) , B< 5 ) ,N, I prnt ,Opo , I Ioq , If n, I J oin.Nj oin 
4888 FOR I-M TO N 
489Z Ao-A<I,M) 

4896 FOR J-h TO N 

4900 A(I,J)-A(I # J)/Ao 

4904 NEXT J 

4908 C(I)-C(I)/Ao 

49 1 Z NEXT I 

49 IB SUBEND 

49Z0 SUB Subtract ( K ) 

49Z4 COM /Cply/ A< 10. 10) .C< 10) ,B< 5) ,N,Iprnt ,Opo. IloQ.Ifn.I join.Njoin 

49Z8 FOR I-K TO N 

493Z FOR J-K-l TO N 

4936 A(I.J)-A(K-1 ,J>-A(I,J> 

4940 NEXT J 

4944 C( I )*C< K- 1 )-C( I > 

4948 NEXT I 
49 5Z SUBEND 
4956 SUB PI in 

4960 COM /Cply/ A< 10,10) ,C( 10) ,B< 5> ,N, Iprnt ,Opo, Ilog, Ifn. I join. Njoin 
4964 COM /Xxyy/ Xx ( Z5 ) , Yy( 25 > 

4968 PRINTER IS 705 
497Z BEEP 

4976 INPUT "WANT TO PLOT Uo vs Vu7 < 1 -Y ,0-N) " . Iuo 
4980 IF Iuo-0 THEN 
4984 BEEP 

4988 INPUT "SELECT <0-h/h0X sane tube, 1 -h< HF >/h< en) " , Irt 
499Z BEEP 

4996 INPUT "SELECT h/h RATIO < 1-FILE ,0- COMPUTED)" , Ihrat 
5000 IF Ihrat *0 THEN 
5004 BEEP 

5008 INPUT "WHICH Tsat ( 1 -6.7,0— 2 . Z >" , Isat 

501 Z END IF 

5016 Xnin-0 

50Z0 Xnax-10 

50Z4 Xstep-Z 

50Z8 IF Irt-0 THEN 

503Z Ynin-0 

5036 Ynax-1.4 

5040 Ystep-.Z 

5044 ELSE 

5048 Ynin-0 

505Z Ynax-15 

5056 Ystep»5 

5060 END IF 

5064 ELSE 

5068 Opo«Z 

507Z Ynin-0 

5076 Ynax-IZ 

5080 Ystep*3 

5084 Xnin-0 

5088 Xnax*4 

509Z Xstep-1 

5096 END IF 

5100 IF I hr at -1 THEN 

5104 Ynin-0 

5108 Ynax-18 

5 1 1 Z Ystep-3 

5116 Xnin-0 
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5120 

5124 

5128 

5132 

5136 

5140 

5144 

5148 

5152 

5156 

5160 

5164 

5168 

5172 

5176 

5180 

5184 

5188 

5192 

5196 

5200 

5204 

5208 

52 12 

5216 

5220 

5224 

5228 

5232 

52 36 

5240 

5244 

5248 

5252 

5256 

5260 

5264 

5268 

5272 

5276 

5280 

5284 

5288 

5292 

5296 

5300 

5304 

5308 

5312 

5316 

5320 

5324 

5328 

5332 

5336 

5340 

5344 

5348 

5352 



X«ax-9 

Xtstep-2 
END IF 
BEEP 

PRINT "INiSPIUP 2 300,2200,8300,68001 M 
PRINT "SC 0, 100,0, 100?TL 2,0i" 

Sfx-100/(X«ax-X«in) 

Sfy-100/< Yaax-Yain) 

PRINT "PU 0,0 PD" 

FOR Xa-Xwin TO Xwax STEP Xstep 

X-< Xa-X«in)*Sfx 

PRINT "PR" i X , " ,0» XT i " 

NEXT Xa 

PRINT "PR 100,0*PU|" 

PRINT "PU PR 0,0 PD" 

FOR Ya-Yain TO Yaax STEP Yetep 
Y-< Ya-Y«in)*Sfy 
PRINT "PR 0, " i Y t " YT" 

NEXT Ya 

PRINT "PR 0,100 TL 0 2" 

FOR Xa-Xain TO Xaax STEP Xetep 

X-< Xa-XP»in>«Sfx 

PRINT "PR" \ X." ,100; XT" 

NEXT Xa 

PRINT "PR 100,100 PU PR 100,0 PD" 

FOR Ya-Yain TO Ynax STEP Yetep 

Y-< Ya-Y«in)*Sfy 

PRINT "PD PR 100," ,Y,"YT" 

NEXT Ya 

PRINT "PR 100,100 PU" 

PRINT "PR 0.-2 SR 1 .5,2" 

FOR Xa-Xain TO Xaax STEP Xetep 
X-< Xa-X*in)*Sfx 
PRINT "PR" i X, " ,0i " 

IF Iuo-0 THEN PRINT "CP -2,-1 iLB" » Xa* ”" 

IF Iuo-1 THEN PRINT "CP - 1 .5,-1 iLB" i Xai"" 

NEXT Xa 

PRINT "PU PR 0,0" 

FOR Ya-Y«in TO Ynax STEP Yetep 
IF RBS( Ya)<1 .E-5 THEN Ya-0 
Y-< Ya- Yain) «Sf y 
PRINT "PR 0 , " ; Y , " “ 

IF Iuo-0 THEN PRINT "CP -4 ,- . 25; LB" » Yai "" 

IF Iuo-1 THEN PRINT "CP -3.- . 25»LB" »Ya* 

NEXT Ya 

Xlabel$-"Oil Percent" 

IF Iuo-0 THEN 
IF Irt-0 THEN 
Ylabel»-"h/h0X" 

ELSE 

Ylabelt-"h/he«ooth M 
END IF 

PRINT "SR 1.5,2 |PU PR 50,-10 CP" i-LEN< XI abelt )/2 t“0»LB" ;Xlabel$i" M 
PRINT "PR -11,50 CP 0."t-LEN( Y1 abel • ) /2«5/6» "DI 0, 1 iLB" i Ylabeli; ““ 

PRINT "CP 0,0" 

ELSE 

PRINT " SP0i SP2" 

PRINT "SR 1 . 2 ,2 . 4 iPU PR -8,35|DI 0,1»LBUiPR 1,0.5»LBofPR -1,0.5«LB (kU/n 
PRINT "PR - 1 .0. 5t SR 1,1. 5t LB2 t SR l.5,2»PR .5,.5»LB.*PR .5,0|LBK>" 
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535G 

536® 

5364 

5368 

5372 

5376 

538® 

5384 

5388 

5392 

53361 

54®® 

54 ®4 

54®8 ! 

5412 

5416 

542® 

5424 

5428 

5432 

5436 

544® 

5444 

5448 

5452 

5456 

546® 

5464 

5468 

5472 

5476 

548® 

5484 

5488 

5492 

5496 

550® 

55®4 

5 5® 8 

5512 

5516 

552® 

5524 

552 8 

5532 

5536 

554® 

5544 

5548 

5552 

5556 

556® 

5564 

5568 

5572 

5576 

558® 

5584 

5588 

5592 



PRINT "PA 42 , - 1 ® t DI 1.®tLBV-,PR .4.-1,LBu;PR 1,.5tLB <m/a) 
PRINT " SP®; SP 1 " 

END IF 



Ipn-® 

BEEP 

INPUT "WANT TO PLOT DATA FROM A FILE < 1-Y.®-N)7" ,Okp 
Icn-® 

IF Okp-1 THEN 
BEEP 

INPUT "ENTER THE NAME OF THE DATA FILE" , D_f ile$ 

IF Iuo-® THEN 
BEEP 

INPUT "SELECT <®-LINEAR. 1 -L06< X , Y . I log 
END IF 

ASSIGN ©File TO D_file© 

BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER". Md 
BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" .Npaira 

IF Iuo-® AND Ihrat-® THEN 

BEEP 

INPUT "ENTER DESIRED HEAT FLUX" f Q 

END IF 

BEEP 

PRINTER IS 1 

PRINT USING "4X. ""Select a eymbol : 



PRINT USING "4X f "" 1 
PRINT USING "4X , 3 
PRINT USING "4X f ""5 
PRINT USING "4X , ""6 
PRINT USING "4X f ""7 
INPUT Syn 
PRINTER IS 7®5 
PRINT "PU DI" 

IF Syn-1 THEN PRINT 
IF Syn- 2 THEN PRINT 
IF Syn- 3 THEN PRINT 
Nn-4 

IF I log- 1 THEN Nn- 1 

IF Md> 1 THEN 

FOR 1-1 TO <Md-1) 

ENTER ©Fil et Xa, Ya 
NEXT I 
END IF 

IF Ihrat-® THEN 
Ql-Q 

IF I log- 1 THEN O-LOG(Q) 
END IF 



Plue eign'" 
Square""" 



Star 2 
Circle 4 
Ronbus 

Right-eide-up triangle" 
Up- aide- down triangle"" 



SM»" 

SH+" 

SMo" 



FOR 1-1 TO Npaira 
IF Iuo-® AND Ihrat-® THEN 
ENTER ©Filet Xa,B< O 
Ya-B(® > 

FOR K-1 TO Nn 
Ya-Ya+B< K ) • Q ^ K 
NEXT K 
END IF 

IF Iuo-1 OR Ihrat- 1 THEN 
ENTER ©Filet Xa, Ya 
IF Iuo-1 THEN Ya-Ya/1®®® 
END IF 
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5596 

5600 

5604 

5608 

5612 

5616 

5620 

5624 

562 8 

5632 

5636 

5640 

5644 

5648 

5652 

5656 

5660 

5664 

5668 

5672 

5676 

5680 

5684 

5688 

5692 

5696 

5700 

5704 

5708 

5712 

5716 

5720 

5724 

572 8 

5732 

5736 

5740 

5744 

5748 

5752 

5756 

5760 

5764 

5768 

5772 

5776 

5780 

5784 

5788 

5792 

5796 

5800 

5804 

5808 

5812 

5816 

5820 

5824 

5828 

5832 



IF Iuo-0 AND I hr at -0 THEN 
IF I loo- 1 THEN Ya-EXP(Ya) 
IF I loo-® THEN Ya-Ql/Ya 
IF Irt-0 THEN 
IF Xa-0 THEN 



Yo-Ya 

Ya-l 

ELSE 

Ya-Ya/Yo 
END IF 
ELSE 

HaR-FNHs«ooth< Q.Xa.Isat ) 

Ya-Ya/He« 

END IF 

END IF 

Xx( 1-1 >-Xa 

Yy(I-l )■ Ya 

X-(Xa-Xnin>*Sfx 

Y-< Ya-Y*in)*Sfy 

IF Syn>3 THEN PRINT "SH" 

IF Syw<4 THEN PRINT "SR 1.4,2. 4" 

PRINT *'PA” ,X ,Y t *'" 

IF Sy«>3 THEN PRINT "SR 1.2 ,1.6" 

IF Syn-4 THEN PRINT M UC2 , 4 , 99 ,0,-8 , -4 f 0 t 0 t 8 ,4 ,0, i" 

IF Syn-5 THEN PRINT "UC3,0,99 ,-3,-6,-3,6,3,6,3,-6r 
IF Syn-6 THEN PRINT "UC0, 5. 3 t 99 , 3,-8 f -6 f 0 , 3, 8| M 
IF Syn-7 THEN PRINT "UC0.-5. 3,99 ,-3,8 t 6 f 0 f -3 ,-8| * 

NEXT I 
BEEP 

ASSI6N CFil e TO • 

END IF 

PRINT M PU SM" 

BEEP 

INPUT "UANT TO PLOT A POLYNOMIAL ( 1-Y ,0-N >?" ,Okp 

IF Okp-1 THEN 

BEEP 

PRINTER IS 1 

PRINT USING "4X , ""Select lino types 1 — 

PRINT USIN6 M 6X , “ M 0 Solid lino MMM 

PRINT USING ”6X t " M 1 Dashed 

PRINT USING M 6X t “ M 2 , , ,5 Lonoer line - dash 

INPUT Ipn 
PRINTER IS 70S 
BEEP 

INPUT "SELECT (0-LINEAR, l-L06< X V Y) )" t IloQ 

Iprnt- l 

CALL Poly 

IF I u o" 1 THEN 

BEEP 

INPUT "DESIRE TO SET X Lower and Upper Linit <1-Y,0-N>?" 

IF Ixlin-0 THEN 

Xll-0 

Xul-7 

END IF 



IF Ixlin-1 THEN 
BEEP 

INPUT "ENTER X Lower LiRit",Xll 
BEEP 

INPUT "ENTER X Upper LiRit",Xul 
END IF 



Ixl in 
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5836 

5840 

5844 

5848 

5852 

5856 

5860 

5864 

5868 

587Z 

5876 

5880 

5884 

5888 

589Z 

5896 

5900 

5904 

5908 

59 1 Z 

5916 

59Z0 

59Z4 

59Z8 

59 3Z 

5936 

5940 

5944 

5948 

595Z 

5956 

5960 

5964 

5968 

5972 

5976 

5980 

5984 

5988 

599Z 

5996 

6000 

6004 

6008 

60! Z 

6016 

60Z0 

60Z4 

6028 

6032 

6036 

6040 

6044 

6048 

605Z 

6056 

6060 

6064 

6068 

607Z 



END IF 

FOR Xo-Xll TO Xul STEP Xetep/ZS 

Icn-Icn* I 

Yo-FNPoly(Xa) 

IF Iuo-1 THEN Ya-Ya/1000 

Y-< Ya-Y«in)*Sf y 

X-< Xa-Xnin )*Sf x 

IF Y<0 THEN Y-0 

IF Y> 100 THEN 60T0 5908 

Pu-0 

IF Ipn-1 THEN Idf-Icn MOD Z 
IF Ipn-Z THEN Idf-Icn MOD 4 
IF Ipn-3 THEN Idf-Icn MOD 8 

IF Ipn-4 THEN Idf-Icn MOD 16 

IF Ipn-S THEN Idf-Icn MOD 3Z 

IF Idf-1 THEN Pu-1 

IF Pu-0 THEN PRINT "PA" , X , Y , "PD" 

IF Pu-1 THEN PRINT "PA" .X. Y,"PU" 

NEXT Xa 
PRINT "PU" 

GOTO 537Z 
END IF 
BEEP 

INPUT “WANT TO QUIT ( 1 -Y ,0-N>?" , Iquit 
IF Iqult-1 THEN 5940 
60T0 537Z 
PRINT "PU SP0" 

SUBEND 

SUB State 

PRINTER IS 701 

J-0 

K-0 

BEEP 

IF Iplot-1 THEN ASSIGN «File TO P.flleS 
BEEP 

INPUT "LAST RUN No?< 0-QUIT >" ,Nn 

IF Nn-0 THEN 61 Z4 

Nn-Nn-J 

Sx-0 

Sy-0 

Sz-0 

Sxs-0 

Sys-0 

S 2 8-0 

FOR 1-1 TO Nn 
J-J+1 

ENTER CFiletO.T 

H-Q/T 

Sx-Sx+Q 

Sxe-Sxe+Q*2 

Sy-Sy+T 

Sys-Sye*T*2 

S2-SZ+H 

Szs-Sze-fH*Z 

NEXT I 

Oave-Sx/Nn 

Tavo-Sy/Nn 

Havo-Sz/Nn 

Sdovq- SOR( ABS< < Nn* Sx s - Sx * Z > / < Nn* < Nn- 1 > ) > > 
$devt-SQR< ABS( ( Nn*Sys-Sy* 2 >/ < Nn* < Nn- 1 > > > > 
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6076 Sdevh-SQR< ABS< ( Nn« Sz a- Sz * Z ) / ( Nn* ( Nn- 1 ) ) ) ) 

6060 Sh- 1 00* Sde vh/Have 
6084 Sq-lOO*Sdevq/Qave 
6088 St-100*Sdevt/Tave 
609 Z IF K-1 THEN 6116 
6096 PRINT 

6100 PRINT USING " 1 1 X , ""DATA FILE: , 14ft" t Filet 
6104 PRINT 

6108 PRINT USING " 1 1 X , “"RUN Htube SdevH Qdp SdevQ Thetab SdevT"" 

6112 K-1 

6116 PRINT USING M 1 1 X , DO, Z< ZX ,D. 3DE . 1 X . 3D. ZD) . ZX ,00. 3D. 1 X . 3D. ZD" i J , Have , Sh , Qave 

.Sq.Tave ,St 

61 ZO GOTO 597Z 

61Z4 ASSIGN CFllel TO • 

61Z8 PRINTER IS 1 
6132 SUBEND 
6136 SUB Coef 

6140 COM /Cply/ A( 10, 10) ,C( 10) ,B(5) ,N, Iprnt ,0po . I log , If n, IJoin.Njoin 
6144 BEEP 

6148 I*#>UT “61 VE A NAME FOR CROSS-PLOT FILE",Cpft 
61 SZ CREATE BOAT Cpft.Z 
6156 ASSIGN OFile TO CpfS 
6160 BEEP 

6164 INPUT "SELECT ( O-LI NEAR, 1-L0G< X . Y ))", Uog 
6168 BEEP 

6172 If#>UT “ENTER OIL PERCENT (- 1 -STOP > M .Bop 

6176 IF Bop<0 THEN 619Z 

6180 CALL Poly 

6184 OUTPUT GFilei Bop,B( • ) 

6188 GOTO 6168 

6192 ASSIGN CFile TO * 

6196 SUBEND 

6200 SUB Uileon( Cf ,Ci ) 

6204 COM /Uil/ DZ t Di .Do.L.Lu.Kcu 
6208 DIMEnf(IZ) 

62121 ULI SON PLOT SUBROUTINE DETERMINE CF AND Cl 
62 1 6 BEEP 

6220 If#>UT "ENTER DATA FILE NAME" .Filet 

6224 BEEP 

6228 PRINTER IS 1 



6232 


PRINT 


USING "4X , 


"“Select option:“"' 




6236 


PRINT 


USING "4X , 


"" 0 Vary Cf 


and 


C i " “ " 


6240 


PRINT 


USING "4X, 


1 Fix Cf 


Vary Ci“"" 


6244 


PRINT 


USING "4X , 


"" 2 Vary Cf 


Fix 


Ci 


6248 


IfPUT 


"ENTER OPTION", I cfix 







6252 PRINTER IS 701 

6256 IF IcflK-0 THEN 6Z7Z 

6Z60 IF Icf ix>0 THEN BEEP 

6264 IF Icfix-1 THEN IM^UT “ENTER Cf " .Cef 

6268 IF Icfix-Z THEN IftfHJT “ENTER CI M .Cl 

6272 PRINTER IS 1 

6276 INPUT "Want To Vary Coef f ?( 1 -Y ,0-N>" , Iccoef 
6280 IF Iccoef- 1 THEN ItfHJT “ENTER COEFF" ,R 
6284 PRINTER IS 701 

6288 IF Icfix-0 OR Icfix-Z THEN Cfa-.004 

6292 IF Icfix-1 THEN Cfa-Cef 

6296 Cia-Ci 

6300 Xn- . 8 

6304 Fr-.3 
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£308 Jj-0 
6312 Rr-3. 

6316 IF Iccoef ■ 1 THEN Rr-R 

6320 PRINTER IS 1 

6324 PRINT Do.Di.Kcu 

6328 ASSI6N CFile TO FileS 

6332 ENTER CFile»Nrun,Dt .Ldtcl , Ldtc2 , lit 

6336 Ru-Do*LOG< Do/Di )/<2«Kcu) 

6340 Sx-0 

6344 Sy-0 

6348 Sxy-0 

6352 Sx2-0 

6356 Sy2-0 

6360 FOR 1-1 TO Nrun 

6364 ENTER CFlleiBop ,TIMEt< Told) . Emf ( • ) ,Fms 
6368! CONVERT EMF * S TO TEMPERATURE 
6372 FOR J-1 TO 5 
6376 T< J )-FNTvsv< Emf < J ) > 

6380 NEXT J 

6384 Tsat-( T< 1 )+T<2))».5 
6388 Tavg-T(5) 

6392 6rad-37. 9853+ . 10438B«Tavg 

6396 Tdrop-Emf < 7 )• 1 . E+6/< 10.«6rad) 

6400 Tavgc-T< 5)-Tdrop* .5 

6404 IF ABS< T avg- T avgc ) > . 0 1 THEN 

6408 Tavg-< Tavg+Tavgc )• . 5 

6412 60T0 6392 

6416 END IF 

6420! 

6424! Compute properties of water 
6428 Kw-FNKw(Tavg) 

6432 Muwa-FNHuw< Tavg) 

6436 Cpw-FNCpw( Tavg) 

6440 Pru-FNPrw(Tavg) 

6444 Rhow-FNRhow( Tavg ) 

6448! 

6452! Compute properties of Freon- 114 

6456 Lmtd-Tdr op/L06( < T< 5)-Tsat )/< T< S)-Tdrop-Tsat ) ) 

6460 IF JJ-0 THEN 
6464 Tu-Tsat+Fr«Lmtd 
6468 Thetab-Tu-Tsat 
6472 Jj-1 
6476 ENO IF 
6480 Tf-< Tw+Teat )*.5 
6484 Rho-FNRho< Tf ) 

6488 Hu-FNMu< Tf ) 

6492 K-FNK(Tf) 

6496 Cp-FNCp(Tf) 

6500 Beta-FNBeta( Tf ) 

6504 Hfg-FNHfg(Tsat) 

6508 Ni-Mu/Rho 
6512 Alpha-K/< Rho*Cp ) 

6516 Pr-Ni/Alpha 
6520! 

6524! Analysis 
65281 COMPUTE MOOT 
6532 A-PI#<Do-2-Di-2>/4 

6536 P-PI#Do 

6540 Mdot-3.9657E-3+Fms*< 3.61955E-3-Fms«< 8. B2006E-6-Fns» < 1 .23688E-7-Fms«4. 31 897 

E- 10) ) ) 
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6544 Q-Hdot*Cpu*Tdrop 

65481 COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
6552 1 FOR UNENHANCED END(S) 

6556 Hbar-190 

6560 Fe-(Hbar*P/(Kcu*A> >“.5*Lu 
6564 Tanh-FNTanh<Fe> 

6568 Theta-Thetab*Tanh/Fe 

6572 Xx-( 9 . 81 *Beta*Thet eb*Do*3»Tanh/( Fe»Ni * Alpha ) > * . 1 66667 
6576 Yy-( 14 ( .SS9/Pr r < 9/ 1 6 > >* < 8/Z7) 

6580 Hbarc-K/Do*< .6+. 387*Xx/Yy)*2 

6584 IF ABS< ( Hbar-Hbarc )/Hbar >>.001 THEN 

6588 Hbar*-< Hbar+Hbarc) • . 5 

6592 60T0 6560 

6596 END IF 

6600' 

6604! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENDS 

6608 Q1 -( Hbar*P*Kcu* A ) * . 5*Thetab*Tanh 

6612 Qc-Q-2«Q1 

6616 Ae-PI»D2*L 

6620! COMPUTE ACTUAL HEAT FLUX 

6624 Odp-Qc/Ae 

6628 IF Icfix-0 OR Icfix>1 THEN Cef- 1 /Cf *( 1 . /Rr ) 

6632 Thetab-Cef/Cp*Hfg*<Qdp/<Hu*Hfg>*< .014/< 9.81 *Rho> >* .5>*< 1/Rr)*Pr*1 .7 

6636 Ho-Qdp/Thetab 

6640 O*ega-Ho/Cf 

6644 Uo«Q/(Pl*Do*L»L*td> 

6648 Vu-Mdot/(Rhou*PI»Di*2/4> 

6652 Rew-Rhou*Vu*Di/Muua 
6656 Tui-Tu+Q*Ru/<PI*Do*L) 

6660 6ana-Ku/Di*Reu* .8*Pru*( 1/3. )•( Muua/FNMuu< Tui ) ) A . 14 

66641 PRINTER IS 1 

6668 Yu-( 1 . /Uo-Ru )*0«ega 

6672 Xu-Q«ega*Do/( 6a«a*Di > 

6676 Sx-Sx+Xu 

6680 Sy-Sy+Yu 

6684 Sxy-Sxy4Yu*Xu 

6688 Sx2-Sx2+Xu*Xu 

6692 Sy2-Sy24Yu»Yu 

6696 NEXT I 

6700 ASSI6N CFile TO • 

6704 M-< Sx*Sy-Nr*un»Sxy )/( Sx*Sx-Nrun* Sx2 > 

6708 C-(Sy-Sx*M)/Nrun 

6712 IF Icfix-0 OR Icfix-3 OR Icfix-4 THEN 

6716 Cic-1/h 

6720 Cfc-1/C 

6724 END IF 

6728 IF Icfix-1 THEN 

6732 Cic-1/H 

6736 Cfc-Cf 

6740 END IF 

6744 IF Icfix-2 THEN 

6748 Cic-Ci 

6752 Cfc-1/C 

6756 END IF 

6760 IF ABS((Ci-Cic)/Cic)>.001 OR ABS( ( Cf-Cf c >/Cf c) > .001 THEN 
6764 Ci-( Ci+Cic)* .5 
6768 Cf-(Cf+Cfc>*.5 
6772 PRINTER IS 1 

6776 PRINT USING "2X, MM Caf - " M .M2.3DE t 2X t M ** Ci - t MZ . 3DE" iCsf ,Ci 

6780 PRINTER IS 701 
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6784 

6788 

679Z 

6796 

680® 

6804 

6808 

6812 

6816 

6820 

6824 

6828 

6832 

6836 

6840 

6844 

6848 

6852 

6856 

6860 

6864 

6868 

6872 

6876 

6880 

6884 

6888 

6892 

6896 

6900 

6904 

6908 

6912 

6916 

6920 

6924 

6928 

6932 

6936 

6940 

6944 

6948 

6952 

6956 

6960 

6964 

6968 

6972 

6976 

6980 

6984 

6988 

6992 

6996 

7000 

7004 

7008 

7012 

7016 

7020 



SOTO 6328 
END IF 
PRINT 

PRINTER IS 701 

PRINT USING "23X t " " Cf Ci 

PRINT USING "8X , " "ASSUMED "" ,M2 . 3DE , 3X , MZ . 3DE” i Cf a ,Cla 

PRINT USIN6 "SX , " "CALCULATED "" ,MZ . 3DE , 3X ,MZ . 3DE" ; Csf ,Ci 
PRINT 

Sun2-Sy2-2*M* Sxy-2«C*Sy+M* 2» Sx2+2*M*C*Sx+Nrun*C A 2 
PRINT USING " 10X , " "Sun of Squares - "" ,2. 3DE" i Sum2 
PRINT USING "10X. ""Coefficient - "" ,D.DDD M »Rr 
SUBEND 

DEF FNMuu(T) 

6-247. 8/<T+133. IS) 



Mu-2. 4E-5* 10-6 
RETURN Mu 
FNEND 

DEF FNCpu(T) 

Cpu-4.2n20858-T*(2.26826E-3-TM4.4Z361E-S+2.71428E-7*T)> 

RETURN Cpw* 1000 
FNEND 

DEF FNRhou(T) 

Ro-999 . 52946+T* { .01 269-T*< 5.48251 3E-3-T* 1 . 2341 47E-5) ) 

RETURN Ro 
FNEND 

DEF FNPru(T) 

Pru-FNCpu< T )*FNMuu( T )/FNKu< T ) 

RETURN Pru 
FNEND 

DEF FNKu(T) 

X-( T+273. 15)/273. 15 

Ku— .92247+XM2.8395-XM 1 .8007-XM .52577- .07344«X)> ) 

RETURN Ku 
FNEND 
SUB Plot 

COM /Cply/ 6( 10,10) f C( 10) f B( 5) ,Nop, Iprnt .Opo.Ilog.Ifn.I join.Njoin 
DIM Bs( 3 ) 

INTEGER li 
PRINTER IS 1 



Idv-0 



BEEP 

INPUT "LIKE DEF6ULT V6LUES FOR PLOT ( 1-Y ,0-N)?" , Idv 
Opo-0 



BEEP 

PRINT USING "4X , " 
PRINT USING "6X, " 
PRINT USING "6X f " 
PRINT USING "6X," 
INPUT Opo 
BEEP 

INPUT "SELECT UNITS 
PRINTER IS 705 
IF IdvOl THEN 
BEEP 
INPUT 
BEEP 
INPUT 
BEEP 
INPUT 
BEEP 



h vereufi q""" 

( O-SI , 1 "ENGLISH)" ,Iun 

OF CYCLES FOR X-6XIS" 
OF CYCLES FOR Y-6XIS" 



"Select Option:”"" 

"0 q versus delto-T"”" 
"1 h versus delta- T 

"2 



"ENTER NUMBER 
"ENTER NUMBER 
"ENTER MIN X-V6LUE (MULTIPLE OF 10) 



Cx 

Cy 

,Xnin 
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7024 

7028 

7032 

7036 

7040 

7044 

7048 

7052 

7056 

7060 

7064 

7068 

7072 

7076 

7080 

7081 

7084 

7088 

7092 

7096 

7097 

7098 

7099 

7100 

7101 

7102 

7104 

7105 

7108 

71 12 

7116 

7120 

7124 

7128 

7132 

7136 

7140 

7144 

7148 

7152 

7156 

7160 

7164 

7168 

7172 

7176 

7180 

7184 

7188 

7192 

7196 

7200 

7204 

7208 

7212 

7216 

7220 

7224 

7228 

72 32 



INPUT “ENTER MIN Y-VALUE (MULTIPLE OF 10)" ,Y*in 
ELSE 

IF Opo-0 THEN 

Cy-3 

Cx-3 

Xmn- . 1 

Y*in-100 

END IF 

IF Opo-1 THEN 

Cy-2 

Cx-2 

Xmn- . 1 

Ymn- 1 00 

END IF 

IF Opo-2 THEN 
IF Iun-0 THEN 
Cy-3 
Cx-2 

Xmin-1000 
Ymn- 1 00 
ELSE 
Cy-3 
Cx-3 

X* in- 1 00 
Y«in- 10 
END IF 
END IF 
END IF 
BEEP 

PRINT "IN;SP1sIP 2300, 2200.8300, 6800*“ 

PRINT "SC 0, 100.0, 100sTL 2,0s" 

Sfx-100/Cx 
Sf y- 1 00/Cy 
BEEP 

INPUT "WANT TO BY-PASS CAGE? ( 1 -Y ,0-N>" , Ibyp 
IF Ibyp- 1 THEN 7632 
PRINT "PU 0,0 PD" 

Nn-9 

FOR 1-1 TO Cx*1 

Xot-X«in*10"(I-n 

IF I-Cx+1 THEN Nn-1 

FOR J-1 TO Nn 

IF J-1 THEN PRINT "TL 2 0" 

IF J-2 THEN PRINT "TL 1 0" 

Xo-Xot* J 

X-LGTC Xa/X«in )*Sf x 
PRINT "PA" » X, ",0s XT s M 
NEXT J 
NEXT I 

PRINT "PA 100,0s PU i M 
PRINT "PU PA 0,0 PD" 

Nn-9 

FOR 1-1 TO Cy+1 
Yot-Y«in*10*<I-1 > 

IF I-Cy+1 THEN Nn-1 

FOR J-1 TO Nn 

IF J-1 THEN PRINT "TL 2 0" 

IF J-2 THEN PRINT "TL 1 0" 

Yo-Yot * J 

Y-L6T( Yo/Ynin )»Sf y 



198 



72 36 

72 40 

7244 

7248 

7252 

725B 

7260 

7284 

7268 

7272 

7276 

7280 

7284 

7288 

7292 

7296 

7300 

7304 

7308 

7312 

7316 

7320 

7324 

7328 

7332 

7336 

7340 

7344 

7348 

7352 

7356 

7360 

7364 

7368 

7372 

7376 

7380 

7384 

7388 

7392 

7396 

7400 

7404 

7408 

7412 

7416 

7420 

7424 

7428 

7432 

7436 

7440 

7444 

7448 

7452 

7456 

7460 

7464 

7468 

7472 



PRINT “PA 0 , " i Y , " YT" 

NEXT J 
NEXT I 

PRINT "PA 0,100 TL 0 2" 

Nn-9 

FOR 1-1 TO Cx*1 
Xat-Xwin* 1 0* < I - 1 ) 

IF I-Cx+1 THEN Nn-1 

FOR J-1 TO Nn 

IF J-1 THEN PRINT “TL 0 2“ 

IF J>1 THEN PRINT “TL 0 1“ 

Xa-Xat* J 

X-L6T( Xa/X«in)»5f x 
PRINT "PA“»X, “,100s XT M 
NEXT J 
NEXT I 

PRINT “PA 100,100 PU PA 100,0 PD" 

Nn-9 

FOR 1-1 TO Cy+ 1 
Yat-Ynin* 10*< I- 1 > 

IF I-Cy+1 THEN Nn-1 

FOR J-1 TO Nn 

IF J-1 THEN PRINT “TL 0 2" 

IF J>1 THEN PRINT “TL 0 1 M 
Ya-Yat * J 

Y-L6T< Ya/Ynin )*Sf y 
PRINT “PD PA 100, " ,Y ,“YT M 
NEXT J 
NEXT I 

PRINT “PA 100,100 PU" 

PRINT “PA 0,-2 SR 1.5,2“ 

I i-L6T< Xnin) 

FOR 1-1 TO Cx-H 
Xa-Xnin*10"(I-1 > 

X-L6T< Xa/Xnin) *Sf x 
PRINT “PA" i X," ,0i“ 

IF Ii>-0 THEN PRINT "CP -2 ,-2iLB10iPR -2 ,2 |LB“ i I ii"“ 
IF li<0 THEN PRINT “CP -2 .-2|LB10|PR 0.2 i LB" i I i i 
Ii-Ii+1 
NEXT I 

PRINT “PU PA 0,0“ 

Ii-L6T< Ynin) 

Y10-10 

FOR 1-1 TO Cy+1 
Ya-Ynin«10 A < 1-1 > 

Y-L6T< Ya/Y«in>*Sfy 
PRINT "PA 0," i Y ,"“ 

PRINT "CP -4,-.25iLB10iPR -2,2iLB"»Iii"“ 

Ii-Ii+1 
NEXT I 
BEEP 

INPUT "WANT USE DEFAULT LABELS < 1-Y ,0-N>?" , Idl 

IF IdlOl THEN 

BEEP 

INPUT “ENTER X-LABEL" .Xlabelt 
BEEP 

INPUT "ENTER Y-LABEL" , Ylabelt 
END IF 

IF Opo<2 THEN 

PRINT "SR 1 ,2 iPU PA 40,-14i“ 
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7476 PRINT *'LB( T ;PR -1.6,3 PD PR I.Z,0 PUiPR . 5 , -4; LBuoiPR .S.li" 

7480 PRINT "LB- T i PR . 5 , - 1 iLBaat i PR .5,1;" 

7484 IF Iun-0 THEN 
7488 PRINT "LB) <K)" 

7492 ELSE 

7496 PRINT "LB) <F)" 

7500 END IF 
7504 END IF 
7508 IF Opo-2 THEN 
7512 IF Iun-0 THEN 

7516 PRINT "SR 1.5.2|PU PA 40,-14;LBq (U/niSR 1,t.5iPR 0.5,1 iLB2iSR 1.5.2«PR 
0. 5,-1 i LB)" 

7520 ELSE 

7524 PRINT "SR 1.5,2iPU PA 34,-14iLBq (Btu/hnPR .5,.5iLB.|PR .5,-.5i" 

7528 PRINT "LBf t ;PR .5,1iSR 1,1.5iLB2|SR 1.5,2iPR .5,-liLB);" 

7532 END IF 
7536 END IF 
7540 IF Opo-0 THEN 
7544 IF Iun-0 THEN 

7548 PRINT "SR 1.5,2iPU PA -12.40|DI 0,hLBq <U/mPR -1,0.5 iSR 1,1.5|LB2;SR 1 
.5 ,2|PR 1 , .5»LB) M 
7552 ELSE 

7556 PRINT "SR 1.5 t 2;PU PA -12.32?DI 0,1iLBq (Btu/hnPR - .5, .5iLB. iPR .5,.5i" 
7560 PRINT "LBf i i SR 1,1.5|PR -1,.5|LB2|PR 1,.5iSR 1.5,2iLB)" 

7564 END IF 
7568 END IF 
7572 IF Opo>0 THEN 
7576 IF Iun-0 THEN 

7580! PRINT "SR 1.5,2iPU PA -12,38«DI 0,1iLBh (U/«iPR -1,.5|SR 1,1.5 iLB2iSR 1. 
5 ,2 i PR .5, .5i" 

7584 PRINT "SR 1.2,2. 4, PU PA -12.37|DI 0,1iLBh»PR 1 ,0.5iLBoiPR -1,0.5iLB (U/n 

7588 PRINT "PR -1,.5|SR 1,1.5|LB2|SR 1.5,2|PR .5,.5iLB.iPR . 5,0|LBK )" 

7592 ELSE 

7596 PRINT "SR 1.5.2iPU PA -12,28iDI 0,1iLBh (Btu/hnPR - .5, .5|LB. |PR .5,.5i" 

7600 PRINT "LBftiPR -1,.5iSR 1,1.5 iLB2iSR 1.5,2iPR .5,.5iLB.iPR .5,.5jLBF)" 

7604 END IF 

7608 END IF 

7612 IF Idl-0 THEN 

7616 PRINT "SR 1.5,2|PU PA 50,-16 CP" i -LEN< XI abel« ) /2 \ "01 LB" i XI abel *\ "" 

7620 PRINT "PA -14,50 CP 0 , " i -LEN( Ylabel t )/2«5/6j "DI 0 , 1 i LB" i Y1 abel • i "" 

7624 PRINT "CP 0,0 DI" 

7628 END IF 
7632 Ipn-0 
7636 X11-1.E+6 

7640 Xul — 1.E+6 
7644 Icn-0 
7648 Ifn-0 
7652 Ijoin-1 
7656 BEEP 

7660 INPUT "WANT TO PLOT DATA FROM A FILE < 1 -Y ,0-N)7" , Ok 
7664 IF Ok-1 THEN 
7668 BEEP 

7672 INPUT "ENTER THE NAME OF THE DATA FILE" ,D_f ile« 

7676 ASSIGN «File TO D.fileS 
7680 BEEP 
7684 BEEP 

7688 INPUT "ENTER THE BE6INNING RUN NUMBER" ,Md 
7692 BEEP 

7696 INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" .Npoira 
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77001 

7704! 

7708 

7712 

7716 

7720 

7724 

7728 

7732 

7736 

7740 

7744 

7748 

7752 

7756 

7760 

7764 

7768 

7772 

7776 

7780 

7784 

7788 

7792 

7796 

7800 

7804 

7808 

7812 

7816 

7820 

7824 

7828 

7832 

7836 

7840 

7844 

7848 

7852 

7856 

7860 

7864 

7868 

7872 

7876 

7880 

7884 

7888 

7892 

7896 

7900 

7904 

7908 

7912 

7916 

7920 

7924 

7928 

7932 

7936 



BEEP 

INPUT “CONNECT DATA UITH LINE < 1-Y ,0-N)?" . Icl 
BEEP 

PRINTER IS 1 

PRINT USING "4X # ""Select a ©ynbol : * 

PRINT USING "6X , " " 1 Star 2 Plus siQn 

PRINT USING "6X , "”3 Circle 4 Square 

PRINT USING "6X.""5 Ronbus 

PRINT USING "6X,""6 Right-side-up triangle 

PRINT USING * , 6X t ""7 Up-side-doun triangle 

INPUT Syn 
PRINTER IS 705 
PRINT "PU DI" 

IF Syw-1 THEN PRINT "Sh* M 
IF Syw-2 THEN PRINT "SMV 
IF Syn-3 THEN PRINT "Sho" 

IF Md> 1 THEN 
FOR 1-1 TO < Hd- 1 ) 



ENTER OFileiYa.Xa 
NEXT I 
END IF 

FOR 1-1 TO Npair® 

ENTER OFiletYa.Xa 

IF 1-1 THEN Ql-Ya 

IF I-Npairs THEN Q2-Ya 

IF Opo-t THEN Ya-Ya/Xe 

IF Opo-2 THEN 

Q-Ya 

Ya-Ya/Xa 



Xa-Q 
END IF 

IF Xa<Xll THEN Xll-Xa 
IF Xa>Xul THEN Xul-Xa 
IF Iun-1 THEN 



IF Opo<2 THEN Xa-Xa*1.8 
IF Opo>0 THEN Ya-Ya» . 1761 
IF Opo-0 THEN Ya-Ya» .317 
IF Opo-2 THEN Xa-Xe*.317 
END IF 



X-L6T ( Xa/X*in)*Sf x 



Y-LGT< Ya/Ynin)»Sf y 
KJ-0 

CALL SyabtX.Y.Syw.Icl ,Kj > 

60T0 7924 

IF Syn>3 THEN PRINT 
IF Syfi<4 THEN PRINT 
IF Icl-0 THEN 
PRINT "PA" ,X ,Y, "" 

ELSE 

PRINT "PA" t X, Y. "PD" 

END IF 

IF Syw>3 THEN PRINT 
IF Syw-4 THEN PRINT 
IF Syn-5 THEN PRINT 
IF Syn-6 THEN PRINT 
IF Syn-7 THEN PRINT 
NEXT I 
PRINT "PU" 

BEEP 

INPUT "WANT TO LABEL? ( 1 -Y ,0-N>" , Ilab 



"SH" 

"SR 1. 4,2.4" 



"SR 1.2 ,1.6" 

"UC2 ,4, 99,0, -8, -4,0,0,8,4,0i" 
N UC3 a 0 a 9S , -3,-6, -3,6, 3.6.3,-Br 
"000,5.3,99,3,-8,-6,0,3,81" 
"UC0.-5. 3, 99, -3, 8, 6,0, -3, -8i" 



\ 
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7940 


IF Uab-l THEN 




7944 


PRINT "SP0* SPZ” 




7948 


BEEP 




795Z 


IF Kleb-0 THEN 




7956 


Xlab-5 




7960 


Ylab-85 




7964 


INPUT “ENTER INITIAL X,Y LOCATIONS" .Xlab.Yli 


796B 


Xtt-Xlab-5 




7972 


Ytt-Ylab+8 




7976 


PRINT "SR 1 . t.5" 




7980 


PRINT "SMiPft" ,Xtt , Ytt ,"LB X Heat F 


7984 


Ytt-Ytt-3 




7988 


PRINT “PA" , Xtt , Yt t , "LB Oil Flux 


Name 1 


799Z 


IF Syn-1 THEN PRINT "SM«" 




7996 


IF Sy*-Z THEN PRINT “SM+" 




8000 


IF Syw-3 THEN PRINT "SHo" 




8004 


K lab- 1 




8008 


END IF 




8012 


Kj-1 




8016 


CALL Synb( XI ab , Ylab , Syn, Icl ,K j ) 




8020 


PRINT "SR 1 , I .SiSM” 




8024 


IF Sy*<4 THEN PRINT “PR 2,0“ 




8028 


BEEP 




8032 


INPUT “ENTER BOP". Bop 




8036 


IF Bop< 10 THEN PRINT “PR 2 ,0iLB“ \ Bop; “ 




8040 


IF Bop>9 THEN PRINT “PR . 5 ,0aB“ * Bop i “ 


" 


8044 


Ihf-0 




8048 


IF Q1>QZ THEN Ihf-1 




8052 


IF Ihf-0 THEN PRINT “PR 4.0*LBInc" 




8056 


IF Ihf-1 THEN PRINT “PR 4,0»LBDec“ 




8060 


PRINT “PR Z,0*LB"»D_f ile«»"“ 




8064 


PRINT "SP0; SP 1 i SR 1.5,2" 




8068 


Ylab-Ylab-5 




8072 


END IF 




8076 


BEEP 




8080 


ASSI6N fFile TO • 




8084 


X11-X11/1.Z 




8088 


Xul-Xul* 1 . 2 




8092! 


60T0 8040 




8096 


END IF 




8100 


PRINT "PU SM M 




8104 


BEEP 




8108 


INPUT “WANT TO PLOT A POLYNOMIAL ( 1-Y, 


0-N)?' 


81 12 


IF 6o_on-1 THEN 




8116 


BEEP 




8120 


PRINTER IS 1 




8124 


PRINT USING “4X, "“Select line type: 




8128 


PRINT USING “6X,“ “0 Solid line“““ 




8132 


PRINT USIN6 “ 6X , “ “ 1 Daehed 




8136 


PRINT USING “6X,““Z , , ,5 Longer line - 


dash’" 


8140 


INPUT Ipn 




8144 


PRINTER IS 705 




8148 


BEEP 




8152 


INPUT "SELECT ( 0-LIN, 1 -LOG< X . Y ) )“ , I log 


8156 


Iprnt-1 




8160 


CALL Poly 




8164 


IF Ifn-0 THEN 




8168 


BEEP 




8172 


INPUT “ENTER NUMBER OF FILES TO JOIN?" 


.Njoii 


8176 


END IF 
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818 ® 

8184 

8188 

8 1 SZ 

8196 

8200 

8204 

8208 

8212 

8216 

822 ® 

8224 

8228 

82 32 

8236 

824® 

8244 

8248 

8252 

8256 

826® 

8264 

8268 

8272 

8276 

828® 

8284 

8288 

8292 

8296 

83®® 

8304 

8 3® 8 

8312 

8316 

832® 

8324 

8328 

8332 

8336 

834® 

8344 

8348 

8352 

8356 

836® 

8364 

8368 

8372 

8376 

838® 

8384 

8388 

8392 

8396 

84®0 

84 ®4 

84® 8 

8412 

8416 



I ) oin-0 

IF Ifn<Njovn THEN Ijoin-1 
IF Ifn-0 OR Ijoin-1 THEN 
FOR Ij-0 TO 3 
Ba( I j >-Bs< I j )+B( I j ) 

NEXT Ij 
If n-If n* 1 
END IF 

IF Njoin-Ifn THEN 
FOR Ij-® TO 3 
B( I j )-Bc< Ij )/Nj oin 
B*< Ij >-® 

NEXT Ij 
Ifn-® 

ELSE 

GOTO 7656 
END IF 
BEEP 

INPUT “ENTER Y LOUER AND UPPER LIMITS" , Y1 1 , Yul 
FOR Xx-0 TO Cx STEP Cx/2®® 

Xa-X«in* 1®*Xx 

IF Xa<Xl 1 OR Xa>Xul THEN 8372 

Icn-Icn4l 

Pu-® 

IF Ipn-1 THEN Idf-Icn MOD 2 

IF Ipn-Z THEN Idf-Icn MOD 4 

IF Ipn-3 THEN Idf-Icn MOD 8 

IF Ipn-4 THEN Idf-Icn MOD 16 

IF Ipn-5 THEN Idf-Icn MOD 28 

IF Idf-1 THEN Pu— 1 

IF Opo-® THEN Ya-FNPoly( Xa) 

IF Opo- 2 AND I loo*® THEN Ya-Xa/FNPoly< Xa) 

IF Opo- 2 AND Ilog-1 THEN Ya-FNPoly< Xa) 

IF Opo* 1 THEN Ya-FNPoly(Xa) 

IF Ya<Y«in THEN 8372 
IF Ya<Yl 1 OR Ya>Yul THEN 8372 
IF Iun-1 THEN 
IF Opo<2 THEN Xa-Xa* 1 . 8 
IF Opo>® THEN Ya-Ya* . 1761 
IF Opo-® THEN Ya-Ya* .317 
IF Opo-2 THEN Xa-Xa*. 317 
END IF 

Y-L6T< Ya/Ynin>*Sfy 
X-LGT(Xa/Xnin)*Sfx 
IF Y<® THEN Y-® 

IF Y>1 ®® THEN GOTO 8372 

IF Pu-® THEN PRINT "PA" , X,Y ,-PD M 

IF Pu— 1 THEN PRINT "PA" ,X. Y,"PU" 

NEXT Xx 
PRINT "PU" 

GOTO 7656 
END IF 
BEEP 

INPUT "WANT TO PLOT REILLY* S DATA? < 1-Y.®-N>" .Irly 

IF Opo-® OR Opo- 1 THEN 

Xll-3 

Xul-2® 

END IF 

IF Opo-2 THEN 
Xll-1®®®® 
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84Z0 

84Z4 

84Z8 

843Z 

8436 

8440 

8444 

8448 

845Z 

8456 

8460 

8464 

8468 

847Z 

8476 

8480 

8484 

8488 

849Z 

8496 

8500 

8504 

8508 

851 Z 

8516 

85Z0 

85Z4 

85Z8 

853Z 

8536 

8540 

8544 

8548 

855Z 

8556 

8560 

8564 

8568 

857Z 

8576 

8580 

8584 

8588 

859Z 

8596 

8600 

8604 

8608 

861 Z 

8616 

86Z0 

86Z4 

86Z8 

863Z 

8636 

8640 

8644 

8648 

865Z 

8656 



Xul-100000 
END IF 

IF Irly-1 THEN 

Y11-Z0 

Yul-70 

BEEP 

INPUT “ENTER LOUER AND UPPER Y-LIMITS FOR PLOTTING" ,Y11 ,Yul 

FOR Xx-0 TO Cx STEP Cx/ZOO 

Xa-Xnin»10*Xx 

IF Xa<Xl 1 OR Xa>Xul THEN 8508 
Xl-LOG(Xa) 

IF Opo-O THEN Yl — Z . 540Z837E-1 4X1 •< 4 .97Z01 5 1 -XI *Z .51 34787E- 1 ) 

IF Opo-1 THEN Yl — Z . 540Z837E- 1 4X1 *< 3. 97Z01 5 1 - XI *Z . 5 1 34787E- 1 ) 

IF Opo-Z THEN Yl — 3.7073801E- 1 4X1»< 8.7ZS91 90E- 1 -Xl»6.88Z684ZE-3) 
Ya-EXP(Yl) 

Y-L6T< Ya/Y«in>*Sfy 

X- L6T<Xa/Xnin)«Sfx 
Ipu-0 

IF Y<YI 1 THEN Ipu-1 

IF Y>Yul THEN 60T0 8508 

IF Ipu-0 THEN PRINT "PA" . X . Y , “PO" 

IF Ipu-1 THEN PRINT “PA" , X , Y . “PU" 

NEXT Xx 
PRINT "PU" 

END IF 
BEEP 

INPUT “WANT TO PLOT ROHSENOU CORRELATION? < 1-Y ,0-N)" .Iroha 

IF Iroha-1 THEN 

YI1-15 

Yul-80 

BEEP 

INPUT "ENTER Taat ( Daa CT.Taat 

Caf- .0040 

BEEP 

INPUT “ENTER Caf ( DEF-0.004)" ,Caf 

Tf -Taat 42 

FOR Xx-0 TO Cx STEP Cx/ZOO 
Xa-X«in»10“Xx 

IF Xa<Xl 1 OR Xa>Xul THEN 8684 

XI- LOG(Xa) 

IF Opo<Z THEN Tf -Taat 4Xa/Z 
Rho-FNRho< Tf ) 

K-FNK< Tf ) 

Hu-FNHu< T f ) 

Cp-FNCp< Tf ) 

HfQ-FNHfQ(Taat) 

Ni-Hu/Rho 

Pr-Cp«Hu/K 

0«eoa-Caf«Hfo/Cp*( ( . 014/(9. 81*Rho> P .5/<Hu«Hf 0 > P< 1 ./3>«Pr*1 .7 

IF Opo-O THEN Ya-( Xa/Onoga) *“ 3 

IF Opo-1 THEN Ya-( Xa/ Onega ) * 3/Xa 

IF Opo-Z THEN Ya-Xa*( 2 ./3) /Omega 

IF Opo-Z THEN 

Tfc-Taat4Xa/Ya«.S 

IF ABS< Tf-Tf c)>.01 THEN 

Tf-< Tf 4Tf c )• .5 

GOTO 8584 

END IF 

END IF 

Y-L6T(Ya/Y«in)«Sf y 
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8660 

r Jb64 

8668 

867Z 

8676 

8680 

8684 

8688 

869 Z 

8696 

8700 

8704 

8708 

871 Z 

8716 

87Z0 

87Z4 

87ZB 

873Z 

8736 

8740 

8744 

8748 

875Z 

8756 

8760 

8764 

8768 

8772 

8776 

8780 

8784 

8788 

879 Z 

8796 

8800 

8804 

8808 

881 Z 

8816 

88Z0 

88Z4 

8828 

8832 

8836 

8840 

8844 

8848 

8852 

8856 

8860 

8864 

8868 

8872 

8876 

8880 

8884 

8888 

8900 

8904 



X-LGT( Xo/Xwin) *Sf x 
Ipu-0 

IF Y<Y1 1 THEN Ipu-1 

IF Y>Yul THEN 8684 

IF Ipu-0 THEN PRINT "PA" ,X , Y . "PO" 

IF Ipu-1 THEN PRINT "PA" . X , Y . "PIT 
NEXT Xx 
PRINT "PU" 

END IF 
BEEP 

INPUT "WANT TO QUIT ( 1 -Y.O-N) " , Iqt 
IF Iqt-1 THEN 8712 
GOTO 7636 

PRINT “PU PA 0.0 SPO" 

SUBEND 

SUB Sy«b( X , Y , Syn , Icl ,K j ) 

IF Syw>3 THEN PRINT “SM" 

IF 5y«<4 THEN PRINT “SR 1.4. 2. 4“ 

Yad-0 

IF Kj-1 THEN Yad-.8 
IF Icl-0 THEN 
PRINT “PA" ,X, Y+Yad,“" 

ELSE 

PRINT “PA" , X , Y+Yad . "PD" 

END IF 

IF Sy«>3 THEN PRINT “SR 1.2. 1.6“ 

IF Syw-4 THEN PRINT “UCZ ,4 ,99 .0.-8.-4 ,0,0 .8,4 ,0i " 

IF Sy«-5 THEN PRINT "UC3.0.99 .- 3.-6 .-3,6 , 3 , 6.3. -6i " 

IF Syw-6 THEN PRINT “UCO.5.3,99 ,3.-8.-6.0,3.8»“ 

IF Syw-7 THEN PRINT “UC0,-5.3,99,-3,8,6.0,-3,-8i“ 

IF Kj-1 THEN PRINT "SMiPR 0.-.8" 

SUBEND 
SUB Purg 
BEEP 

I1PUT "ENTER FILE NAME TO BE DELETED" .Filet 

PURGE Filet 

GOTO 8792 

SUBEND 

SUB Tdcn 

COM /Cc/ C< 7 ) , leal 
DIM Ewf< 1 ) 

DATA 0. 10086091 .25727.94369.-767345.8295,78025595.81 
DATA -9247486589 , 6. 9768BE+ 11,-2. 661 9ZE+ 1 3 , 3. 94078E+ 1 4 
READ C<*> 

BEEP 

IffHJT “GIVE A NAME FOR FILE TO BE CREATED" .Filet 
BEEP 

irtHJT “SELECT TUBE <0-UH, 1-HF.Z-UT)" ,Ht 
BEEP 

INPUT "SELECT THERMOCOUPLE TYPE ( O-NEU. 1 -OLD>“ , leal 

IF Itt <2 THEN Di-.OI 27 

CREATE BDAT Filet, 4 

ASSI6N CFile TO Filet 

OUTPUT tFileiltt 

J-0 

BEEP 

OUTPUT 709 i “TIME DATE " 

ENTER 709 i Di 
PRINTER IS 1 
PRINT 
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8908 PRINT " Month, date and tine: ” ; DATE$( Dt ) , TI ME$( Dt ) 

8912 PRINT 

8916 PRINT USING ,, 10X. ,,M Fme Tin Tev Vy*2 Tdrop' 

8920 IF K-0 THEN 
8924 PRINTER IS 701 
8928 PRINT 

8932 PRINT M Month, date and ti*e: “ * DATEK Dt ) ,TI ME*( Dt ) 



89 36 


IF Itt-0 THEN 


PRINT 


USING 


M 10X .’’“Tube 


Type: 




Uielend Snooth*' 


8940 


IF Itt-1 THEN 


PRINT 


USING 


M 10X t M,, Tube 


Type: 




High Flux 


8944 


IF Itt-2 THEN 


PRINT 


USING 


" 10X, M,, Tube 


Type: 




Turbo-B"" " 


8948 


PRINT 














8952 


PRINT USING " 


10X, ,,M 


Frs 


Tin 


Tev 


Vy* 2 


Tdrop 


8956 


PRINTER IS 1 















8960 K-1 

8964 ENO IF 
8968 BEEP 

8972 INPUT “ENTER FLOUMETER READING” ,F«a 
8976 OUTPUT 709; M AR AF0 AL4 VR1 M 
8980 FOR L-0 TO 4 
8984 OUTPUT 709 r AS SA“ 

8988 IF L>0 AND L<4 THEN 9020 
8992 S-0 

8996 FOR 1-0 TO 9 
9000 ENTER 709 i E 
9004 S-S+E 
9008 NEXT I 

9012 IF L-0 THEN Enf < 0)-ABS( S/ 10) 

9016 IF L-4 THEN E«f ( 1 )-ABS( S/ 10) 

9020 NEXT L 

9024 OUTPUT 709 t " AR AF20 AL20 VR1" 

9028 OUTPUT 709 ; “AS SA n 

9032 Etp-0 

9036 FOR 1-0 TO 9 

9040 ENTER 709 ;Et 

9044 Etp-Etp+Et 

9048 NEXT I 

9052 Etp-Etp/10 

9056 Tin-FNTvsv( Enf ( 1 ) > 

9060 Tov-FNTvav( Enf ( 0) ) 

9064 6rad-37. 9853+ . 104 388* Tin 

9068 Mdot-3.9657E-3+Fna*( 3.61955E-3-F«a*( 8. 82006E~6“Fna»( 1 ,23688E-7-F*a*4. 31897 

E - 10) ) ) 

9072 Vu-Mdot/( 1 000»PI »Di * 2 ) *4 
9076 Tdrop- Etp» 1 . E + 6/( 10*Grod) 

9080 PRINT USING M 1 0X , 3( DO. D0.4X > , 1 X , Z . DO ,4X ,MZ. 40" i F ro , T in ,Tev , Vu* 2 , Tdrop 
9084 BEEP 

9088 INPUT “WANT TO ACCEPT THIS DATA SET7 < 1 -Y ,0-N >" ,0k 

9092 J-J+1 

9096 IF Ok -0 THEN 

9100 J-J-1 

9104 60T0 8968 

9108 ELSE 

9112 OUTPUT ®F il c*Frb ,E nf ( * ) ,Etp 
9116 PRINTER IS 701 

9120 PRINT USING “ 10X , 3( DO. D0.4X ) , 1 X , Z . DO ,4X ,MZ . 4D“ i Fna , Tin ,Tev , Vy*2 ,Tdrop 
9124 PRINTER IS 1 
9128 BEEP 

9132 INPUT "WILL THERE BE ANOTHER DATA SET? < 1 - Y ,0-N ) " ,6o_on 
9136 IF 6o_on- 1 THEN 8968 
9140 END IF 



206 



9144 ASSIGN fiFile TO * 

9148 PRINTER IS 701 
915Z PRINT 

9156 PRINT USIN6 "10X. M "NOTE: "" .ZZ,"" data sets ore stored in file 



llet 

9160 PRINTER IS 1 
9164 SUBEND 
9168 SUB Uoprt 
917Z PRINTER IS 1 
9176 BEEP 

9180 INPUT “Enter Uo File None" , Filet 
9184 BEEP 

9188 INPUT "Number of Data Rune M f Nrun 
9 1 9Z INPUT "Do You Want a Plot File?( 1-Y.0-N)" .Iplot 
9196 BEEP 

9Z00 IF Iplot-1 THEN 

9Z04 INPUT "6ive Plot File Name" ,P_f ilet 

9Z08 CREATE BOAT P_filet,4 

9 Z 1 Z ASSIGN ©Plot TO P_filet 

9Z16 END IF 

9ZZ0 PRINTER IS 701 

9ZZ4 PRINT 

9ZZ8 PRINT 

9Z3Z PRINT USING M 10X f MM Water Vel Uo 

9Z36 ASSI6N ®File TO FileS 

9Z40 IF Iplot-1 THEN ASSIGN tFilel TO P_filet 
9Z44 FOR 1-1 TO Nrun 
9Z48 ENTER CFileiVu.Uo 

9Z5Z IF Iplot-1 THEN OUTPUT «File1|Vu,Uo 
9Z5B PRINT USING " ISX.D.OD.BX.HZ.SDE-iVu.Uo 
9Z60 NEXT I 
9Z64 ASSIGN tFile TO • 

9Z68 ASSI6N CFllel TO • 

9Z7Z PRINT USING "10X, MM NOTE: " M ,ZZ f M " data sets are stored in file 
n.Filet 

9 Z76 IF Iplot-1 THEN 

9Z80 PRINT USING "10X. M "NOTE: M, \ZZ f "" X-Y Paire are etored in file 



n,P_f ilet 



9Z84 
9Z88 
9Z9Z 
9Z96 
9300 
9304 
9 308 
931 Z 
9316 
93Z0 
93Z4 
93Z8 
933Z 
9336 
9 340 
9344 
9 345 

9346 

9347 

9348 
9 35Z 
935G 



END IF 

PRINTER IS 1 
SUBEND 
SUB Select 



COM /I dp/ I dp 
BEEP 

PRINTER IS 1 

PRINT USING "4X .""Select option: 

PRINT USING "6X, M " 0 Taking data or re-processing previous data 1 

PRINT USING " 6X , " " 1 Plotting data on Log-Log 

PRINT USING "6X, "" Z Plotting data on Linear 

PRINT USING "6X. 3 Hake cross-plot coefft file 

PRINT USING " 6X , " " 4 Re-circulate water 

PRINT USING " 6X , " " 5 Purge 

PRINT USING "6X . " " G T-Drop correction 

PRINT USING " 6X , " " 7 Print Uo File * 

PRINT USING "6X , 8 Hodify X-Y file 

PRINT USING "BX."" 9 Hove 

PRINT USING "6X, "" 10 Conb/Flxup 

INPUT I dp 

IF I dp-0 THEN CALL Hain 
IF Idp-1 THEN CALL Plot 



, 15A"i J.F 



. 1 5A" i Nru 



, 1 5A" ; Nru 
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9360 IF I dp-2 THEN CALL PI in 
9364 IF I dp- 3 THEN CALL Coef 
9368 IF I dp-4 THEN CALL Main 
9372 IF I dp-5 THEN CALL Purg 
9376 IF I dp-6 THEN CALL Tdcn 
9380 IF I dp-7 THEN CALL Uoprt 
938! IF I dp-8 THEN CALL Xynod 

9382 IF I dp-9 THEN CALL Move 

9383 IF I dp- 10 THEN CALL Conb 
9385 SUBEND 

9394 SUB Xynod 
9434 PRINTER IS 1 
9444 BEEP 

9454 INPUT "ENTER FILE NAME" .Filet 
9464 ASSI6N (File! TO Filet 
9474 BEEP 

9484 INPUT "ENTER NUMBER OF X-Y PAIRS" .Np 
9494 BEEP 

9504 INPUT "ENTER NEU FILE NAME" ,File2t 
9514 CREATE BOAT File2t.5 
9524 ASSIGN «File2 TO File2t 
9534 BEEP 

9544 INPUT "ENTER NUMBER OF X-Y PAIRS TO BE DELETED" .Ndel 
9554 IF Ndel-0 THEN 9594 
9564 FOR 1-1 TO Ndel 
9574 BEEP 

9584 INPUT "ENTER DATA SET NUMBER TO BE DELETED" t Nd< I ) 

9594 NEXT I 
9604 FOR J-l TO Np 
9614 ENTER tFilel iX.Y 
9624 FOR 1-1 TO Ndel 
9634 IF Nd( I )- J THEN 9674 
9644 NEXT I 
9654 OUTPUT CFile2iX.Y 
9664 PRINT J,X t Y 

9674 NEXT J 

9675 PRINTER IS 701 
9684 ASSIGN tFilel TO • 

9694 ASSI6N ®File2 TO • 

9704 SUBEND 

9714 SUB Move 

97181 FILE NAME: MOVE 

97221 

9726 DIM Bop( 66) t A( 66) .B< 66) ,C( 66) ,D< 66) f E< 66) ,F< 66) ,6< 66) ,H< 66) , J< 66) ,X(66),L< 
66) ,M( 66) 

9730 DIM Told( 66) ,N( 66) ,Vr( 66) ,Ir( 66) 

9734 BEEP 

9738 INPUT "OLD FILE TO MOVE" ,02_fi let 

9742 ASSI6N fFile2 TO 02_filet 

9746 ENTER ®File2 iNrun.Dt .Ldtcl ,Ldtc2 .lit 

9750 FOR 1-1 TO Nrun 

9754 ENTER «File2 iBop< I ) f Told< I ) 

9758 ENTER «F-ile2 1 A< I ) t B< I ) ,C< I ) ,D( I > .E< I ) ,F( I ) t G< I ) t H< I ) . J< I ) ,K< I ) ,L( I ) ,M( I ) # N 
(I) 

9762 ENTER ®File2 I Vr( I ) , Ir( I ) 

9766 NEXT I 

9770 ASSI6N «File2 TO • 

9774 BEEP 

9778 INPUT "SHIFT DISK AND HIT CONTINUE" ,0k 
9782 BEEP 
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9786 INPUT “INPUT BDAT SIZE** , Size 
9790 CREATE BOAT DZ_f i le$ . Size 
9794 ASSIGN CFilel TO DZ_file« 

9798 OUTPUT CF i leltNrun ,Dt .Ldtcl .LdtcZ , lit 

980Z FOR 1-1 TO Nrun 

9806 OUTPUT CFilel ;Bop( I ) ,Told( I ) 

9810 OUTPUT CFile1;A<I).B(I) i C(I),0(I),E(I>,F(I),G(I> a H(I).J(I),K(I),L(I),M(I), 
N( I ) 

9814 OUTPUT CFilel t Vr( I ) , Ir( I ) 

9818 NEXT I 

98ZZ ASSIGN CFilel TO * 

98Z6! RENAME “TEST** TO DZ.filet 

9830 BEEP Z000..Z 

9834 BEEP 4000,. Z 

9838 BEEP 4000.. Z 

984Z PRINT "DATA FILE MOVED** 

9846 SUBEND 
9874 SUB Comb 
98781 FILE NAME: COMB 
988Z1 

9886 DIM Emf ( 1 Z ) 

9890 BEEP 

9894 INPUT "OLD FILE TO FIXUP" ,DZ_f ileS 
9898 ASSI6N CFileZ TO DZ.fileS 
990Z D1_f i let-" TEST" 

9906 CREATE BOAT D1_fileS.Z0 
9910 ASSI6N CFilel TO D1_filet 

9914 ENTER CFileZ iNrun.Dt .Ldtcl .LdtcZ .Itt 

9915 Nrunm-Z0 

9918 IF K-0 THEN OUTPUT CFilel »Nrunn.Dt .Ldtcl .LdtcZ . Itt 
99ZZ FOR 1-1 TO Nrun 

99Z6 ENTER CFileZ;Bop.Told.Emf< • ) .Vr.Ir 
9930 OUTPUT CFilel; Bop .Told. Emf ( • ) a Vr , I r 
9934 NEXT I 

9938 ASSI6N CFileZ TO • 

994ZI RENAME "TEST" TO DZ_filet 
9946 BEEP Z000..Z 
9950 BEEP 4000.. Z 
9954 BEEP 4000.. Z 
9958 BEEP 

996Z IfPUT "WANT TO ADO ANOTHER FILE ( 1 -Y .0-N>?" .Oko 
9966 IF Ok a-1 THEN 
9970 K-1 

9974 BEEP 

9978 INPUT "GIVE NEU FILE NAME" .Nf lie* 

998Z ASSI6N CFileZ TO NfileS 
9986 60T0 9914 

9990 END IF 

9994 ASSI6N CFileZ TO • 

9998 SUBEND 
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